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ABSTRACT 


This handbook provides estimates of average available solar insola- 
tion to fixed, flat-plate, south-facing collector surfaces at various 
array tilt angles at numerous sites in the United States. This first 
volume contains average daily, total insolation estimates, by month, and 
annual totals for 235 locations. The second and third volumes contain 
the daily profiles by hour nsed to compute the daily totals for the 235 
locations (at selected array tilt angle.s). 

A model that estimates the direct, diffuse, and reflected components 
of total insolation on an hourly, daily, and monthly basis is presented. 

A shadow loss model and a reflector augmentation model providing esti- 
mates of the losses and gains associated with various fixed array geom- 
etries are also described. These models can be used with the insolation 
model provided or with other recorded data. 

A FORTRAN computer program with user's gride is presented. The pro- 
gram can be used to generate additional handbook values or to examine the 
effects of array shadowing and fixed reflector augmentation effects on a 
daily, monthly, or annual basis. Array shadowing depends on location, 
array size, array tilt, array separation, and time. The program can be 
used to examine trade-offs between array spacing and insolation losses 
due to shadowing. The reflector augmentation program can be used to 
examine trade-offs among array size and tilt, separation, and reflector 
tilt to determine the combination of design values that optimize the 
economic objectives or technical criteria of the system. 
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GLOSSARY 


latitude angle 

declination angle (tilt of Earth with respect to plane 
through equator of sun) 

array tilt angle from horizontal 

reflector tilt angle from horizontal 

sun hour angle, 0° at solar noon, negative in the morning 
and positive in the afternoon; the Earth rotates through an 
angle of Aw* 15° each hour 

solar altitude angle 

incident angle between direct beam of sun and normal to 
surface of array 

azimuth angle (0° for north, measured clockwise) 
zenith angle between direct beam of sun and vertical 
solar constant (1.377 kW/tn^) 
sunrise hour angle 

the total daily radiation on a horizontal surface (measured) 
at ground level divided by the total daily extraterrestrial 
radiation on a horizontal surface (computed) (i.e., the 
fraction of the total daily insolation penetrating the 
atmosphere); computed for each day in t’ e month and condensed 
into an average that represents a "typical" day for the morth 
(H/H 0 ) 

the fraction of total daily extraterrestrial radiation on a 
horizontal surface (computed) reaching the Earth's surface as 
diffuse (scattered) radiation on a horizontal surface 
(measured); computed for each day of the month and condensed 
into an average that represents a typical day for the month 
(H d /H 0 ) 

the average daily total radiation reaching a horizontal sur- 
face on the ground (measured) during a given month 

daily total extraterrestrial radiation on a horizontal sur- 
face for a given month 

monthly adjustment to the solar constant to account for 
Earth-sun distance variation 
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m 


B atmospheric absorption coefficient 

Hj average daily total terrestrial diffuse radiation on a 

horisontal surface for a given month 

r<j fraction of total daily extraterrestrial radiation reaching 

the ground as diffuse radiation during a given hour u on a 1 

horisontal surface . 

?dh average hourly terrestrial diffuse radiation on a horisontal 

surface for a given month 

Ton average hourly terrestrial direct normal radiation for a 
given month 

A adjusted average daily terrestrial direct normal radiation 

for a given month 

l D h average hourly terrestrial direct radiation on a horizontal 

surface for a given month 

Ifh average hourly terrestrial total ’’’diation on a horisontal 

surface for a given month 

I 

_ I 

lot average hourly terrestrial direct radiation on a fixed tilted 

surface for a given month 

I<j t average hourly terrestrial diffuse radiation on a fixed 

tilted surface for a given month 

I rC average hourly terrestrial diffuse radiation from ground 

reflectance for a given month 

Lpt average hourly total terrestrial radiation on a fixed tilted 

surface for a given month 

p ground reflectance cf surrounding area 

p' spectral reflectance of rack-mounted arrays 

d array spacing 

a array surface length 

F g (u>) shadow loss fraction during hour angle u> 

R reflector length on reflector augmented racks 

Note: The phrase "average hour ly/daily.. .for a given month" indicates 

that the value represents an average hour or day during the given 
month (not the total number of average hours or days in the month). 
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SECTION t 


INTRODUCTION 


This document is designed for engineers, designers, and home- 
owners. For the homeowner this handbook provides a quick reference 
for obtaining estimates of available solar insolation for numerous 
locations and array tilt angles. For the engineers and designers it 
provides a model and computer program that considers the effects of 
array shadowing and reflector augmentation as design variables. 

The first objective of this publication is to provide a handbook 
reference document for estimating available solar insolation on south- 
facing fixed surfaces. The user interested in such data is referred 
directly to the Appendices. Hourly and daily insolation estimates are 
provided for 235 U.S. sites at a variety of array tilt angles. In 
the event the tilt angle is to be adjusted monthly or fixed at some 
intermediate angle other than the handbook angles, a computer program 
is provided to aid the user with more detailed requirements. The 
computer program is based on a parametric solar insolation model that 
takes into account the site specific solar climate of a particular 
location. (This document does not consider tracking surfaces although 
the computer program in Appendix A could be modified to do so.) The 
computer program enables the engineer/designer to conduct extensive 
sensitivity analyses on the orientation and location of solar 
collectors . 

This publication is directed toward the engineer/designer who is 
interested in average insolation estimates that represent the typical 
solar climate of a particular site. In this sense, the model pre- 
sented here is an "aggregate" model and does not attempt to model the 
atmospheric physics in a detailed fashion. In particular, hourly 
trends within a day are not modeled; that is, a symmetric daily pro- 
file is presented. The latest available input data required by this 
model are provided in Appendix B of this document. 

The second objective of this publication is to address insolation 
problems specific to residential/commercial applications. Specifi- 
cally, it is concerned with the placement of solar collectors on resi- 
dential housing units and commercial properties (shopping centers, 
parking lots, etc.) where the amourt of space available for collectors 
is limited by shadowing problems, e ergy output requirements, and 
costs. The array spacing problem is also relevant to central station 
solar power plants although spacing may not be as critical due to the 
larger land areas involved. The document addresses the relationships 
among available insolation, spacing, and shadowing for south-facing 
surfaces — the costs associated with various spacing arrangements are 
not addressed. The reflector augmentation problem is related to the 
shadowing model. Given a tandem rack array design (sec page 4-2), a 
reflective surface built into the support structure of the arrays 
vields a higher level of insolation at the collector surface. The 
insolation available to such a surface depends not only upon the 
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amount of reflector augmentation, but alao on array apacing and aun 
elevation angles at which shadowing occurs. 

A multipurpose package of programs has been written based on the 
models presented in this report. Details of this package may be found 
in Section VI and Appendix A. For the user interested in south-facing 
fixed-surface data, Appendices C and D contain data for 235 sites. 
Shadowing and reflector augmentation data are not presented due to the 
many possible combinations of input parameters. In order to obtain 
insolation estimates that consider shadowing and reflector 
augmentation, Sections VI and Appendix A should be consulted for 
directions on how to run the program. 

It should be noted that the modeling of the spectral response of 
photovoltaic cells by time of day and type of material has not been 
considered here due to the complexities associated with modeling light 
spectrum changes. This is discussed more in Section II. 

Sections II, III, IV, and V present the models used to generate 
the handbook data and computer program package as well as a discussion 
of the assumptions and accuracv of the techniques. A careful examina- 
tion of these Sections is advised prior to altering the computer 
program. 
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SECTION II 


THE INSOLATION MODEL 


ngg 

There are basically two ways that solar radiation data can be used 
for design purposes in the photovoltair.s area. One way is to use past 
hourly recorded data for a particular site or similar locality and esti- 
mate through computer simulation what the performance of the system and 
its components would have been under these historical conditions. 

Future performance can then be predicted based on this historical anal- 
ysis. The second method is to use data that correspond to specific 
types of days (clear* average* "typical*" etc.) to calculate the per- 
formance of the system and its components. 

The first method can account for the transient nature of the pro- 
cess. However, the use of this technique is limited by the cost and 
complexity associated with the use of detailed models and computer 
programs and the almost complete absence of good hourly data for solar 
radiation on tilted surfaces. The primary source of solar radiation 
data in the United States is the National Weather Service and its 
network of stations. Data collection began in 1902 with the first 
measurements of the direct component of the solar radiation at normal 
incidence being made in Asheville* N.C. 

Until 1977 there were 86 stations covering the United States and the 
West Indies. However, due to errors in data collection and instrument 
calibration, the accuracy of these data was not guaranteed. Reference 1 
gives a complete history of the network. In 1977 a new network with 
anproximately 30-35 stations was established with the aid of the U.S. 
Department of Energy. This network collects only horizontal surface 
(direct, diffuse, and total) data by hour. 

With the inc.'*nsing emphasis on the development of solar energy use, 
much of the past solar radiation data has been rehabilitated. Major 
work has been done in an effor*. to smooth the data record using models 
to fill in gaps due to missing data or improperly calibrated instru- 
ments. The acronym used to describe this new data base is SOLMET 
(Reference 2), which consists of solar and meteorological data for 
27 sites throughout the United States. 

The second method — using solar radiation data for specific types 
of days and orientations — has the advantage that computing insolation 
estimates for surfaces of arbitrary orientation is a simple computa- 
tional procedure However, the input data required are aggregated solar 
insolation measurements. The approach adopted here was to use the 
latest available SOLMET derived input data to calibrate an average or 
typical day insolation model and thus combine the advantages of both 
methods . 
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A. NOMENCLATURE AND ASSUMPTIONS 


The model presented here yields insolation levels for "typical 
days" — a term with specific meaning. The method of procedure may be 
outlined as follows. The amount of solar radiation intensity outside 
the atmosphere is relatively constant , so the "typical" amount selected 
for a given month is the quantity around the midpoint (15th or 16th). 

The proportion of this radiation that actually reaches the Earth's sur- 
face varies from day to day (hour to hour, in fact) depending upon the 
weather characteristics of the location. The amount of solar radiation 
reaching the Earth's surface during a given month, divided by the total 
radiation outside the atmosphere during the given month, is used as an 
estimate of the average or "typical" fraction of radiation reaching the 
gri md. The direct, diffuse, and reflected components of the total 
radiation are derived from this fraction and are referred to as typical 
estimates since they represent the available insolation one could expect 
on any day in the particular month. Thus, all estimates generated by 
the model are estimates of the long-run average insolation. 

A number of assumptions are associated with the use of this model. 
First, it is assumed that most of the diffuse radiation comes from an 
origin near the sun. In other words, the scattering of the solar radia- 
tion is mostly forward scattering. Second, it is assumed that the dif- 
fuse component is uniformly distributed over the sky. In general, for 
clear days there is little variation in diffuse solar radiation across 
the United States.* However, while this may not be true for particular 
partly cloudy days, it is probably sufficient for the "typical" esti- 
mate. It is further assumed that the effects on the average total inso- 
lation of tilting the surface (Reference 3) are accounted for by the 
atmospheric attenuation adjustment contained in the model discussed 
later. 


B. PREVIOUS WORK 

Liu and Jordan (References 4-6) conducted extensive analyses during 
the early 1960s on the available solar radiation data and developed an 
approach that can be used to estimate the available daily solar radia- 
tion for each month of the year for numerous locations in the United 
States and Canada. Using this approach, it is possible to take the 
daily total radiation on a horizontal surface, divide the daily total 
into its direct and diffuse (scattered radiation) components, convert 
each component into hourly values, and then compute the hourly value of 
either component on a surface of any orientation desired. In addition, 
with the availability of the rehabilitated SOLMET data, the calibration 
of this model has been extended and improved. 


♦This can be seen in the relationship between Kf and Kj (page 2-7) 
where the percentage of the total daily extraterrestrial radiation that 
reaches the ground as diffuse ranges from 12% to 19%. 
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This Section presents s version of the Liu and Jordan model modified 
by Kusuda and Ishii at the National Bureau of Standards (Reference 7) to 
account for early morning and late afternoon errors. These errors are 
caused by increased atmospheric absorption of the direct radiation beam 
at extreme sun-hour angles. (During early morning and late afternoon, 
the direct radiation must pass through a longer length of atmosphere.) 
The model was adapted to allow for varying array orientations and the 
algorithm for calculating the total radiation on a tilted surface 
follows. Furthermore, the model has been updated to include the most 
recent insolation input data and the revised solar constant 
(1.377 kW/m^). The model yields hourly solar radiation values for a 
given site location (latitude), month (declination), tilt angle, array 
orientation, and reflectivity. A number of examples are presented in 
Appendix A of this report. 


C. CALCULATION PROCEDURE 

The calculation procedure is a modification of a method developed by 
Kusuda and Ishii, which is based on the historical model of Liu and 
Jordan. For brevity, the approaches are presented in combined form. 
Details may be found in Reference 7. 

The relevant relationships derived by Liu and Jordan (References 4-6) 

are: 

(1) The monthly average daily diffuse radiation on a horizontal 
surface as a function of the monthly average daily total radia- 
tion on a horizontal surface; 

(2) The ratio ,/f the hourly diffuse radiation to the daily diffuse 
radiation (both on a horizontal surface) as a function of the 
number of daylight hours; 

(3) The ratio of the hourly total radiation to the daily total 
radiation (both on a horizontal surface) as a function of the 
number of daylight hours. 

The relationships were based on available solar radiation data from 
Blue Hill, Mass.; Nice, France; Helsingfors, Finland; and Kew 
Observatory, London, England. More recent data have led to minor 
adjustments in the diffuse component relationships (Reference 8). 

However, comparisons of these adjustments are inconclusive due to the 
small differences involved and the limited data. The approach used here 
retains the original Liu and Jordan results for the diffuse component. 

Liu and Jordan (References 6, 9) also compiled for each of 12 months, 
at 80 different locations in the United States and Canada, the "average" 
daily total radiation on a horizontal surface (H) and the ratio Kj ■ 

H/H 0 , where H 0 is the daily total extraterrestrial radiation on a 
horizontal surface calculated for the 15th or 16th of each month. It 
should be noted that the choice of the midpoint of the month is somewhat 
arbitrary since the monthly variation in daily extraterrestrial 


radiation is only +1.5%. Using the rehabilitated SOLMET data base, the 
National Climatic Center has updated the data used for computing the 
fractions (Reference 10). Using regression models, the data base 
has been extended to cover 248 locations. These values were used to 
generate the 235 selected handbook sites in Appendices C and D. Since 
the Kf values are site specific and reflect the cloudiness and weather 
conditions of the site, they more accurately represent the expected 
insolation levels than do the so-called "clear day" estimates of the 
earlier literature in this area. The K-j ratio represents the fraction 
of the total average daily extraterrestrial radiation that is trans- 
mitted through the atmosphere and thus provides a useful index for the 
solar climate of a particular location. 

The step-by-step technique for using the above relationships to com- 
pute the incident solar radiation on a horizontal or tilted surface 
follows. Figures 2-la and 2-lb illustrate the relevant angles involved 
in the procedure. 

Step 1 

The average daily total extraterrestrial radiation (H 0 ) on a horizon- 
tal surface for a given site (latitude) and month (declination angle) 
must be calculated first. This calculation will give the total avail- 
able radiation on a horizontal surface before the radiation beam passes 
through the atmosphere and is broken into its direct and diffuse compo- 
nents. 

The expression is: 

H 0 = — X Y X I c X (cos <j> X cos 6 X sin W g + W s X sin <|> X sin 5 ) 
where 

H 0 = daily total extraterrestrial radiation on a horizontal 
surface (kWh/m^/day) 

Y = ratio of solar radiation intensity at normal incidence 
outside the Earth's atmosphere to the solar constant 
(monthly average values are given in Table 2-1) 

I c ■ solar constant, 1.377 kW/m^ 

(j> = latitude, radians* 

6 * solar declination angle, radiar.s; monthly values are given 
in Table 2-1; the declination is the angular position of 
the sun at solar noon with respect to the plane of the 
equator (north positive) 

W s = sunrise hour angle (radians) = cos - ^ (-tan 4* tan 6 ) 


♦Degrees are converted to radians by multiplying by tt/180. All the 
trigonometric functions used here require radian arguments. 
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Figure 2-lb. Angles Used in Solar Insolation Calculation 
(Modified from Reference 7) 



Table 2-1. Solar Declination (6), the Ratio of Solar Radiation 
Intensity at Normal Incidence Outside Earth's 
Atmosphere to the Solar Constant (Y), and 
Atmospheric Absorption Coefficient (B) (Reference 7) 


Month 

Declination 
(Degrees . Minutes) 

6 

Ratio 

Y 

Atmospheric 

Absorption 

Coefficient 

B 

1 January 

-19.51 

1.0300 

0.142 

2 February 

-10.28 

1.0207 

0.144 

3 Mi»rch 

0.20 

1.0057 

0.156 

4 April 

11.56 

0.9875 

0.180 

5 Msy 

20.14 

0.9727 

0.196 

6 June 

23.27 

0.9670 

0.205 

7 July 

20.26 

0.9692 

0.207 

8 August 

12.03 

0.9785 

0.201 

9 September 

0.37 

0.9945 

0.177 

10 October 

-10.47 

1.0133 

0.160 

11 November 

-19.58 

1.0267 

0.149 

12 December 

-23.27 

1.0327 

0.142 


It should be noted that all times are expressed as angles of the sun's 
position relative to a vertical plane oriented north to south. 

Step 2 

The next step is to obtain the monthly average daily total radiation un 
a horizontal surface, H, from Reference 10 or other available sources 
and calculate the ratio: 


K T - H/H 0 

The value ftp represents the average fraction of daily total extra- 
terrestrial radiation that penetrated the atmosphere. If the value of 
Kj is available, this calculation may be omitted. Appendix B contains 
a table of values for 235 locations in the United States (Reference 10) 
derived from the rehabilitated SOLMET data base. 

Step 3 


Liu and Jordan determined the relationship between the fractions of 
extraterrestrial radiation that were daily diffuse (K<j) and daily 
total radiation (K-p) (on a monthly average basis) for any given 
month. This empirical relationship is given by: 


Kj 

Kd 

0.3 

0.179 

0.4 

0.133 

0.5 

0.188 

0.6 

0.174 

0.7 

0.149 

0.75 

0.125 


For example, if 60% of the average daily extraterrestrial radiation is 
reaching the Earth at a particular location, 17.4% of the average daily 
extraterrestrial radiation which reaches the ground is diffuse 
(scattered by the atmosphere). 

Thu K-p, K(j relationship is used_to determine the daily diffuse radi- 
ation on a horizontal surface, H<j, by: 

H d - Kd x H 0 


Step 4 


The next step is to take the daily vclue above and convert it to an 
hourly value. This is done by calculating the value of a factor, r^, 
from the following equation that gives the f~action, r<j, of average 
daily extraterrestrial radiation that reaches the ground as diffuse 
radiation during a given hour: 
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^ cos u> - cos W # 

r d 24 X sin W - W x cos W 
s s s 

where w * the hour angle corresponding to a given hour of the day 
(solar time), radians. 

Note that the use of this relationship "smooths" the average daily 
diffuse profile into a symmetric function as mentioned earlier. Also 
note that solar time and standard time are different. Solar time is a 
local phenomenon, determined so that solar noon occurs when the sun is 
most nearly overhead. Thus, two observers displaced east and west of 
each other would measure different solar times at the same Pacific 
Standard Time (PST). The following table illustrates the relationship 
between solar time and hour angle: 

Solar Time Hour Angle to (Degrees) 



10 

• 

A.M. 

• 

-30° 


11 

A.M. 

-15° 

Solar Neon: 

12 

Noon 

0 


1 

P.M. 

15° 


2 

P.M. 

30° 


Using rj, the average hourly terrestrial diffuse radiation on a 
horizontal surface (I^h) is the average daily diffuse radiation 
multiplied by the fraction available during hour u*: 

Idh = r d x ^d 


Step 5 


This step involves the calculation of the direct normal component of the 
hourly radiation adjusted to more accurately account for the larger 
amount of atmosphere the beam must pass through near the sunrise and 
sunset hours. The details of this procedure may be found in Reference 7. 


*A small d subscript indicates a diffuse component and a large D 
denotes a direct component. The second letter of the subscript 
denotes whether the surface is horizontal (h), tilted (t), or normal 
to the direct beam (N). A large T denotes total radiation amounts. 


The value of the factor m, which representa a correction for the 
increased atmospheric absorption of the direct normal component , can be 
determined from Table 2-2.* So, the hourly direct normal radiation is 
calculated using: 


Idn "AX exp (-B/cos 9 t ) (with T ■ 0) 

where TpN * average hourly terrestrial direct normal 

radiation (kWh/m 2 ) 

A - (K T - Kd) X f x H 0 (kWh/m 2 ); the 

correction factor m ■ f X H 0 where f is a 
function of latitude and time of year 

B * an atmospheric absorption coefficient 

(see Table 2-1) that describes the varying 
beam attenuation by the atmosphere during 
the year 

■ the angle of incidence between the sun's 

direct beam and a normal to the array surface. 

It is this variable that takes into account 
the azimuthal orientation, a, and tilt angle, t, 
of the surface. 

The equation for 0^ is derived from Reference 12: 


(sin <p x cos T + cos <p x sin T x cos a) x sin 6 

+ (cos <j> X cos t - sin <#> X sin t X cos a) x cos 6 

cos u) - sin T X sin a x cos 6 X sin u 

if cos 0 t > 0 

0 if cos 6 t SO 



*The computer program presented in Section VI of this document generates 
these m factors internally. 
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Table 2 - 2 a . The Value of the o Factor ( kWh / m *) 


LAT . MONTH 


DEG . 


1 

< 

> 


3 


A 


5 


6 


7 

i 

5 


9 

10 

11 

12 

0 . 

1 

.77 

1 

. 7 A 

1 

. 7 A 

1 

.77 

1 

.81 

1 

.83 

1 

.83 

1 

.81 

1 

.77 

1 . 

77 

1 . 

78 

1 . 

78 

1 . 

1 

.77 

1 

. 7 A 

1 

. 7 A 

1 

.77 

1 

.80 

1 

.83 

1 

.83 

1 

.81 

1 

.77 

1 . 

77 

1 . 

78 

1 . 

78 

2 . 

1 

. 77 

1 

. 7 A 

1 

. 7 A 

1 

.77 

1 

.80 

1 

.82 

1 

.82 

1 

.81 

1 

.77 

1 . 

77 

1 . 

78 

1 . 

79 

3 . 

1 

.77 

1 

. 7 A 

1 

. 7 A 

1 

.77 

1 

.80 

1 

.82 

1 

.82 

1 

.80 

1 

.77 

1 . 

77 

1 . 

79 

1 . 

79 

A . 

1 

.78 

1 

. 7 A 

1 

. 7 A 

1 

.77 

1 

.80 

1 

.82 

1 

.82 

1 

.80 

1 

.77 

J . 

77 

1 . 

79 

1 . 

79 

5 . 

1 

.78 

1 

.75 

1 

. 7 A 

1 

.77 

1 

.79 

1 

.81 

1 

.82 

1 

.80 

1 

.77 

1 . 

77 

1 . 

79 

1 . 

80 

6 . 

1 

.78 

1 

.75 

1 

. 7 A 

1 

.77 

1 

.79 

1 

.81 

1 

.81 

1 

.80 

1 

.77 

1 . 

78 

1 . 

80 

1 . 

80 

7 . 

1 

.79 

1 

.75 

1 

. 7 A 

1 

.77 

1 

.79 

1 

.81 

1 

.81 

1 

.80 

1 

.77 

1 . 

78 

1 . 

80 

1 . 

80 

8 . 

1 

.79 

1 

.75 

1 

. 7 A 

1 

.76 

1 

.79 

1 

.81 

1 

.81 

1 

.80 

1 

.77 

1 . 

78 

1 . 

80 

1 . 

81 

9 . 

1 

.79 

1 

.75 

1 

. 7 A 

1 

.76 

1 

.79 

1 

.80 

1 

.81 

1 

.80 

1 

.77 

1 . 

78 

1 . 

81 

1 . 

81 

10 . 

1 

.80 

1 

.76 

1 

. 7 A 

1 

.76 

1 

.78 

1 

.80 

1 

.81 

1 

.80 

1 

.77 

1 . 

79 

1 . 

81 

1 . 

82 

11 . 

1 

.80 

1 

.76 

1 

. 7 A 

1 

.76 

1 

.78 

1 

.80 

1 

.80 

1 

.80 

1 

.78 

1 . 

79 

1 . 

82 

1 . 

82 

12 . 

1 

.81 

1 

.76 

1 

. 7 A 

1 

.76 

1 

.78 

1 

.80 

1 

.80 

1 

.80 

1 

.78 

1 . 

79 

1 . 

82 

1 . 

83 

13 . 

1 

.81 

1 

.76 

1 

.75 

1 

.76 

1 

.78 

1 

.80 

1 

.80 

1 

.80 

1 

.78 

1 . 

80 

1 . 

83 

1 . 

83 

1 A . 

1 

.81 

1 . 

.77 

1 

.75 

1 

.76 

1 

.78 

1 

.80 

1 

.80 

1 

.80 

1 

.78 

1 . 

80 

1 . 

83 

1 . 

5 A 

15 . 

1 

.82 

1 . 

.77 

1 

.75 

1 

.76 

1 

.78 

1 

.79 

1 

.80 

1 

.80 

1 

.78 

1 . 

80 

1 . 

8 A 

1 . 

8 A 

16 . 

1 

.82 

1 , 

.77 

1 

.75 

1 

.77 

1 

.78 

1 

.79 

1 

.80 

1 

.80 

1 

.78 

1 . 

81 

1 . 

8 A 

1 . 

85 

17 . 

1 

.83 

1 , 

. 78 

1 

.75 

1 

.77 

1 

.78 

1 

.79 

1 

.80 

1 

.80 

1 

.79 

1 . 

81 

1 . 

85 

1 . 

86 

18 . 

1 

. 8 A 

1 , 

.78 

1 

.76 

1 

. 77 

1 

.78 

1 

.79 

1 

.80 

1 

.80 

1 

.79 

1 . 

81 

1 . 

85 

1 . 

86 

19 . 

1 

. 8 A 

1 , 

.78 

1 

.76 

1 

.77 

1 

.78 

1 

.79 

1 

.80 

1 

.80 

1 

.79 

1 . 

82 

1 . 

86 

1 . 

87 

20 . 

1 

.85 

1 , 

.79 

1 

.76 

1 

.77 

1 

.78 

1 

.79 

1 

.80 

1 

.80 

1 

.79 

1 . 

82 

1 . 

87 

1 . 

88 

21 . 

1 

.86 

1 

. 79 

1 

.76 

1 

.77 

1 

.78 

1 

.79 

1 

.80 

1 

.81 

1 

.80 

1 . 

83 

1 . 

87 

1 . 

89 

■»** 

(. fm . 

1 

.86 

1 , 

.80 

1 

.77 

1 

. 77 

1 

.78 

1 

.79 

1 

.80 

1 

.81 

1 

.80 

1 . 

83 

1 . 

88 

1 . 

90 

23 . 

1 

.87 

1 , 

.80 

1 

.77 

1 

.77 

1 

.78 

1 

.79 

1 

.80 

1 

.81 

1 

.80 

1 . 

8 A 

1 . 

89 

1 . 

91 

2 A . 

1 

.88 

1 

.81 

1 

.77 

1 

. 77 

1 

.78 

1 

.79 

1 

.80 

1 

.81 

1 

.81 

1 . 

8 A 

1 . 

90 

1 . 

92 

25 . 

1 

.89 

1 . 

.81 

1 

.78 

1 

.78 

1 

.78 

1 

.79 

1 

.80 

1 

.81 

1 

.81 

1 . 

85 

1 . 

91 

1 . 

93 

26 . 

1 

.90 

1 

.82 

1 

.78 

1 

.78 

1 

.78 

1 

.79 

1 

.80 

1 

.82 

1 

.81 

1 . 

86 

1 . 

92 

1 . 

9 A 

27 . 

1 

.91 

1 . 

.82 

1 

.78 

1 

.78 

1 

.78 

1 

.79 

1 

.80 

1 

.82 

1 

.82 

1 . 

86 

1 . 

93 

1 . 

95 

28 . 

1 

.92 

1 . 

.83 

1 

.79 

1 

.78 

1 

.78 

1 

.79 

1 

.80 

1 

.82 

1 

.82 

1 . 

87 

1 . 

9 A 

1 . 

96 

29 . 

1 

.93 

1 

. 8 A 

1 

.79 

1 

.79 

1 

.79 

1 

.79 

1 

.81 

1 

.82 

1 

.83 

1 . 

88 

1 . 

95 

1 . 

98 

30 . 

1 

. 9 A 

1 

. 8 A 

1 

.79 

1 

.79 

1 

.79 

1 

.80 

1 

.81 

1 

.83 

1 

.83 

A . 

89 

1 . 

96 

1 . 

99 

31 . 

1 

.95 

1 

.85 

1 

.80 

1 

.79 

1 

.79 

1 

.80 

1 

.81 

1 

.83 

1 

• 8 A 

1 . 

89 

1 . 

97 

C • 

01 

32 . 

1 

.96 

1 . 

.86 

1 

.80 

1 

.79 

1 

.79 

1 

.80 

1 

.81 

1 

.83 

1 

.SA 

1 . 

90 

1 . 

99 

2 . 

02 

33 . 

1 

.98 

1 

.87 

1 

.81 

1 

.80 

1 

.79 

1 

.80 

1 

.81 

1 

. 8 A 

1 

.85 

1 . 

91 

2 . 

00 

2 . 

0 A 

3 A . 

1 

. 99 

1 . 

.87 

1 

.81 

1 

.80 

1 

.80 

1 

.80 

1 

.82 

1 

.SA 

1 

.85 

1 . 

92 

2 . 

02 

C • 

06 

35 . 

2 

.01 

1 

.88 

1 

.82 

1 

.80 

1 

.80 

1 

.80 

1 

.82 

1 

. 8 A 

1 

.86 

1 . 

93 

2 . 

03 

2 . 

08 

36 . 

2 

.02 

1 

.89 

1 

.83 

1 

.81 

1 

.80 

1 

.81 

1 

.82 

1 

.85 

1 

.87 

1 . 

9 A 

2 . 

05 

2 . 

10 

37 . 

•> 

C 

. 0 A 

1 . 

. 90 

1 

.83 

1 

.81 

1 

.80 

1 

.81 

1 

.83 

1 

.85 

1 

.87 

1 . 

95 

2 . 

07 

2 . 

12 

38 . 

2 

.06 

1 , 

.91 

1 

. 8 A 

1 

.82 

1 

.81 

1 

.81 

1 

.83 

1 

.86 

1 

.88 

1 . 

97 

2 . 

09 

2 . 

15 

39 . 

2 

.08 

1 

.92 

1 

.SA 

1 

.82 

1 

.81 

1 

.82 

1 

.83 

1 

.86 

1 

.89 

1 . 

98 

2 . 

11 

2 . 

18 

AO . 

2 

.10 

1 , 

. 9 A 

1 

.85 

1 

.83 

1 

.81 

1 

.82 

1 

. 8 A 

1 

.87 

1 

.90 

1 . 

99 

2 . 

1 A 

2 . 

21 

Al . 

2 

.13 

1 . 

. 95 

1 

.86 

1 

.83 

1 

.82 

1 

.82 

1 

. 8 A 

1 

.87 

1 

.91 

2 . 

01 

2 . 

16 

2 . 

2 A 

A 2 . 

2 

.15 

1 . 

. 96 

1 

.87 

1 

. 8 A 

1 

.82 

1 

.83 

1 

.SA 

1 

.88 

1 

.92 

2 . 

02 

C . 

19 

2 . 

28 

A 3 . 

2 

.18 

1 

.98 

1 

.88 

1 

. 8 A 

1 

.83 

1 

.83 

1 

.85 

1 

.89 

1 

.93 

2 . 

0 A 

2 . 

22 

2 . 

32 

AA . 

2 

.21 

1 

. 99 

1 

.89 

1 

.85 

1 

.83 

1 

.83 

1 

.85 

1 

.89 

1 

. 9 A 

2 . 

06 

2 . 

26 

2 . 

36 

A 5 . 

2 

.25 

2 

.01 

1 

.89 

1 

.85 

1 

.83 

1 

.SA 

1 

.86 

1 

.90 

1 

.95 

2 . 

08 

2 . 

29 

2 . 

Al 

A 6 . 

2 

.29 

2 

.03 

1 

.90 

1 

.86 

1 

. 8 A 

1 

. 8 A 

1 

.86 

1 

.91 

1 

.96 

2 . 

10 

2 . 

3 A 

2 , 

, A 7 

A 7 . 

2 

.33 

2 

.05 

1 

.92 

1 

.87 

1 

. 8 A 

1 

.85 

1 

.87 

1 

.92 

1 

. 97 

2 . 

12 

2 . 

38 

2 , 

.53 

A 8 . 

2 

.37 

2 

. 07 

1 

.93 

1 

.88 

1 

.85 

1 

.85 

1 

.87 

1 

.92 

1 

.98 

2 . 

1 A 

2 . 

A 3 

2 . 

.60 

A 9 . 

2 

• A 3 

2 

.09 

1 

. 9 A 

1 

.88 

1 

.86 

1 

.86 

1 

.88 

1 

.93 

2 

.00 

2 . 

17 

2 . 

. A 9 

2 . 

.68 

50 . 

2 

. A 8 

2 

.11 

1 

.95 

1 

.89 

1 

.86 

1 

.86 

1 

.89 

1 

. 9 A 

2 

.01 

2 . 

20 

2 . 

55 

2 

.78 
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Table 2-2b 


The Value of the m Factor (BTU/ft 2 ) 


UT. MONTH 


DEG. 

1 

2 

3 

9 

5 

6 

7 

8 

9 

10 

11 

12 

0. 

561. 

552. 

551. 

563. 

573. 

581 . 

581. 

579. 

561. 

561. 

565. 

565. 

1. 

561. 

552. 

551. 

562. 

572. 

580. 

580. 

579. 

561. 

561 . 

566. 

566. 

2. 

562. 

553. 

551. 

562. 

571 . 

579. 

579. 

573. 

561 . 

562. 

566. 

567. 

3. 

563. 

553. 

551. 

561. 

571. 

578. 

578. 

573. 

561. 

562. 

567. 

568. 

9, 

569. 

559. 

551. 

561. 

570. 

577. 

577 . 

572«. 

561. 

563. 

568. 

56 9. 

5. 

565. 

559. 

552. 

561. 

569. 

576. 

576. 

572. 

562. 

563. 

569. 

570. 

6. 

566. 

555. 

552. 

560 . 

568. 

375. 

575. 

572. 

562. 

569. 

570. 

571. 

7. 

567. 

555. 

552. 

560. 

568. 

579. 

575. 

571. 

562. 

565. 

571. 

573. 

8. 

568. 

556. 

552. 

560 . 

567. 

573. 

579. 

571. 

562. 

565. 

572. 

579 . 

9. 

569. 

557. 

553. 

560. 

567 . 

573. 

573. 

571. 

563. 

566 . 

579. 

575. 

10. 

570. 

558. 

553. 

560. 

566. 

572. 

573. 

571 . 

563. 

567. 

575. 

577. 

11. 

572. 

558. 

553. 

560. 

566. 

571. 

573. 

571. 

563. 

568. 

576. 

578. 

12. 

573. 

559. 

559. 

560 . 

565. 

571 . 

572. 

571. 

569. 

569. 

578. 

580. 

13. 

579. 

560. 

559. 

560 . 

565. 

570. 

572. 

571. 

569. 

570 . 

579. 

582. 

19. 

576. 

561 . 

555. 

560. 

565. 

570. 

572. 

571 . 

565. 

571 . 

581. 

589. 

15. 

577 . 

562. 

555. 

560. 

565. 

570 . 

571 . 

571 . 

566. 

572. 

583. 

585. 

16. 

579. 

563. 

556. 

560 . 

565. 

569. 

571 . 

571 . 

566. 

573. 

589. 

587. 

17. 

581. 

569. 

557. 

560 . 

569. 

569. 

571 . 

572. 

567. 

579. 

586 . 

590. 

18. 

583. 

565. 

557. 

561 . 

5o9 . 

569. 

571 . 

572. 

568. 

576. 

588. 

592. 

19. 

585. 

566 . 

558. 

561 . 

569 . 

569. 

571 . 

572. 

569. 

577 . 

590. 

599. 

20. 

587. 

568. 

559. 

561 . 

569. 

568. 

571 . 

573. 

569. 

579. 

592. 

597. 

21. 

589. 

569. 

560. 

562. 

569. 

568. 

571 . 

573. 

570. 

580. 

595. 

599. 

22. 

591 . 

570. 

560. 

562 . 

565. 

568. 

57 3 . 

579. 

571. 

582. 

597. 

602. 

23. 

599. 

572. 

561 . 

563. 

565. 

568. 

57 3 . 

579 . 

572. 

583. 

600 . 

605. 

29. 

596. 

579 . 

562. 

563. 

565. 

568. 

57 t . 

575. 

579. 

585. 

602. 

608. 

25. 

599. 

575. 

563. 

569. 

565. 

569. 

572. 

575 . 

575. 

587. 

605. 

612. 

26 . 

602. 

577 . 

569. 

565. 

565. 

569. 

572. 

576 . 

576 . 

589. 

608. 

615. 

27 . 

605. 

579. 

566. 

565. 

566. 

569. 

572. 

577 . 

577. 

591 . 

612. 

619. 

28. 

608. 

581. 

567. 

566 . 

566 . 

569. 

573. 

578. 

579. 

599. 

615. 

623. 

29. 

611 . 

583. 

568. 

567. 

567 . 

570 . 

573. 

579. 

580. 

596. 

619. 

627. 

30. 

615. 

585. 

569. 

568. 

567. 

570 . 

579. 

580. 

582. 

598. 

622. 

631. 

31. 

619. 

587. 

571. 

568. 

568. 

570. 

579. 

581. 

553. 

601 . 

627. 

636. 

32. 

623. 

590 . 

572. 

569. 

568. 

571 . 

575. 

582. 

585. 

609. 

631. 

691 . 

33. 

627. 

592. 

579. 

570 . 

569. 

572. 

576. 

583. 

587. 

607 . 

635. 

697. 

39. 

632. 

595. 

576. 

572. 

570 . 

572. 

577 . 

589. 

589. 

610. 

690 . 

653. 

35. 

637. 

598. 

577. 

573. 

571 . 

573. 

577. 

585. 

591 . 

613. 

696. 

659. 

36 . 

692. 

601 . 

579. 

579. 

571. 

579. 

578. 

587 . 

593 

617 . 

651 . 

666 . 

37. 

698. 

609 . 

581. 

575. 

572. 

579. 

579. 

588. 

595. 

620 . 

657 . 

679. 

38. 

659. 

607. 

583. 

577 . 

573. 

575. 

580 . 

590 . 

597. 

629 . 

669 . 

682. 

39. 

660. 

611 . 

586 . 

578. 

579. 

576. 

581 . 

591 . 

600. 

628 . 

671 . 

691 . 

90. 

668. 

615. 

588. 

580 . 

575. 

577 . 

582. 

593. 

602. 

633. 

679. 

700 . 

91 . 

675. 

619. 

590 . 

581 . 

577. 

578. 

589. 

595. 

605. 

637 . 

687. 

711. 

92. 

689. 

623. 

593. 

583. 

578. 

579. 

585. 

597. 

608. 

692. 

696. 

722. 

93. 

693. 

628. 

595. 

585. 

579. 

580. 

586. 

599. 

611 . 

698. 

706 . 

735. 

99. 

703. 

633. 

598. 

587. 

581 . 

582. 

588. 

601 . 

619. 

653. 

716. 

799. 

95. 

719. 

638. 

601. 

589. 

582. 

583. 

589. 

603. 

618. 

659. 

728. 

765. 

96. 

726. 

693. 

609. 

591 . 

589. 

585. 

591 . 

605. 

621. 

666. 

791 . 

783. 

97. 

739. 

699. 

608. 

593. 

585. 

586. 

593. 

608. 

625. 

673. 

755. 

802. 

98. 

759. 

656. 

611. 

595. 

587. 

588. 

595. 

611 . 

629. 

680. 

771. 

825. 

99. 

770. 

663. 

615. 

598. 

589. 

589. 

597. 

613. 

639. 

689. 

789. 

851 . 

50. 

788. 

671 . 

619. 

600 . 

591 . 

591. 

599. 

616. 

638. 

698. 

809. 

881 . 
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The azimuth angle, a , is measured clockwise by a north-facing observer 

where I 0° for north orientation, 

a ■ < 90° for east orientation, 

I 180° for south orientation, etc* 
and T - the array surface tilt angle from the horizontal. 

Note that for a south-facing surface at any tilt angle 

cos 0 t ■ sin ( ♦ - T ) X sin 6 ♦ cos ( $ - T ) X coo 6 X cot u> 

Step 6 

The average hourly terrestrial direct radiation on a horizontal surface 
is given by: 

T Dh " T DN X CO< 6 t (with T “ 0) 

Step 7 

The average hourly terrestrial total radiation on a horizontal surface 
is given by the sum of the hourly direct and diffuse components: 

*Th " *Dh + *dh 


Step 8 


The average hourly terrestrial direct radiation on a tilted surface i 
given by: 

I Dt ■ I DN x cos 6 t (with t ■ given array tilt angle) 

Step 9 


The average hourly terrestrial diffuse radiation on a tilted surface is given 
by: 


I 


dt 


1 + COS T - 

2 cih 


Step 10 


Since the array is tilted up at some angle, T , it receives reflected 
radiation from the surrounding area that it faces. This reflected 
radiation will vary depending on whether there is high reflectivity 
(e.g, snow) or low reflectivity (e.g., forest). Tables 2-3 and 2-4 
show average reflectivities for various surroundings. Hunn and 
Calafell (Reference 13) have developed a method for ascertaining these 
average reflectivities. They conclude: 
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1. Snow (freshly fallen or with ice filn) 0.75 | 

2. Water surfaces (relatively large incidence 0.07 | 

angles) 1 

3. Soils (clay, loam, etc.) 0.14 

4. Earth roads 0.04 

5. Coniferoi. ; forest (winter) 0.07 

6. Forests in autumn, ripe field crops, plants 0.26 

7. Weathered blacktop 0.10 

8. Weathered concrete 0.22 

9. Dead leaves 0.30 

10. Dry grass 0.20 

11. Green grass 0.26 

12. Bituminous and gravel roof 0.13 | 

13. Crushed rock surface 0.20 

14. Building surfaces, dark (red brick, dark 0.27 

paints, etc.) 

15. Building surfaces, light (light brick, 0.60 

light paints, etc.) 


IS. Building surfaces, light (light brick, 0.60 

1 ight paints, etc.) 



Table 2-4. Average Reflectivity Valuea for 12 Representative 
Winter Landscapes (Reference 13) 


Surface Description Average 

Reflectivity 


Rural Areas 

Fields with Snow Cover 


1. 

Field with wooded area in background 

0.66,0.73 

2. 

Open field (soil and dry grass) near road 

0.61,0.70 

3. 

Trees dispersed in field 

0.62 

Wooded Areas 


4. 

Conifer forest (with heavy snow covar) 

0.61 

5. 

Deciduous forest (with heavy snow cover) 

0.72 

Water 


6. 

Open water 

0.16 

7. 

Ice/snow-covered water 

0.68 

8. 

Partially open waterway (trees and houses 
in background) 

0.43,0.66 

Urban Areas 


9. 

Urban areas (commercial , institutional) 

0.16,0.38 

10. 

Residential areas (dwellings and roadway) 

0.21,0.35,0.45 

11. 

Educational institution 

0.36,0.42 

12. 

Recreational area (park) 

0.49 


9 
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(1) An average ground reflectivity of 0.6 Co 0.7, similar Co that 
uatd by Liu «nd Jordan (Reference 14), la accuraca for aoat 
rural landaeapaa in winter and whara anow covar ia predomi- 
nanc. Tha exception Co Chia ia for locations adjacanc Co 
opan, or parCially opan, bodiaa of water, whara a conaidarably 
lowar valua of reflectivity ia obcainad. 

(2) For urban araaa, no characCariaCic ground raflaccivicy for 
wincar aay ba apacifiad dua Co cha ride variacion in landacapa 
dacaila poaaibla. Howavar, cha ranga for ground raflaceivicy 
for urban araaa in wincar ia 0.16 eo 0.49, conaidarably lowar 
than eh«c uaad by Liu and Jordan. 

(3) For all landaeapaa conaidarad, ground raflactiviciaa era quite 
aenaicive to cha fraction of field of view in anow covar. 


Thua, the ground reflectivity ia an input paraneter denoted hare by p . 
* 

The average hourly terrestrial reflected radiation from the ground ia 
given by: 



1 - coa T 


x x i, 


Th 


where p ia the average ground reflectance. 


Step 11 


The final atep ia the calculation of the average hourly terreatrial 
total radiation on a tilted surface. Thia total ia the aum of the 
direct, diffuae, and reflected components calculated in Stepa 8-10. 


*Tt * lot ♦ Idt ♦ Irt 

Since the model dividea daily values into hourly intenaitiea (kilowatt* 
per aquare meter), the intenaity over one hour can be integrated to get 
kilowatt-houra per aquare meter (if it ia aaaumed intenaity ia rela- 
tively constant over the hour). So at thia point, the unite of I-j-t 
are kilowatts during hour u per aquare meter. 


D. SPECTRAL CONSIDERATIONS 

As mentioned earlier, the spectral response of solar cells is 
limited in comparison to the spectrum of light. The standard curve for 
the spectral distribution of light was in part established by the JPL 
Mariner Mars probe (References 15 and 16) and is illustrated in 
Figure 2-2. The area under this standard curve represents the solar 
constant I c , defined earlier. Since the establishment of this con- 
stant, an adjustment has been made due to Frolich (Reference 17) 
updating the constant to 1.377 kW/m?. This new value is used 
throughout this document. Superimposed on the same curve is the range 
of a tvpical solar cell relative response curve. As shown, the region 
that the solar cell responds to is smaller than the overall spectrum. 
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IRRADIANCE W 


In fact, the extraterrestrial radiation available to a solar cell can 
be obtained by multiplying the cell spectral response curve by the 
standard curve of Figure 2-2 and calculating the area under the 
resulting curve. 

In order to model the transmission of light through the atmosphere, 
a Rayleigh type equation (Reference 11) could be used to generate 
transmission coefficients for each wavelength. However, the difficulty 
in using a transmission equation is determining values of the 
atmospheric constants in the equations. Furthermore, data that 
represent profiles over t^me of these constants are not generally 
available. 

Thus the use of a non-spectrum oriented model will over-estimate 
the actual radiation to which a solar cell will respond. However, 
since this model estimates available insolation, it is implicitly 
assumed that, ir. the case of photovoltalcs, conversion efficiencies 
must account for this difference when calculating power output from 
available insolation. 



Figure 2-2. Solar Spectral Irradiance, Standard Curve, 
bolar Constant (1.377 kW/m^) 
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PHOTOVOLTAIC CELL RESPONSE FRACTION 



SECTION III 


ARRAY SHADOWING 


A. THE PROBLEM 

The problem of array shadowing arises when arrays are placed too 
close together and the lower part of an array is caught in the shadow 
of the array in front of it. The model presented here is based on a 
model that assumes fixed flat-plate photovoltaic arrays facing south at 
some tilt angle (see Figure 3-1) and sufficiently close to cause shad- 
owing at some point during the year (References 18 and 19).* It is 
further assumed that the array length in the east-west direction is 
large relative to the array height and that cells within the arrays are 
connected in horizontal series strings so that shadowing does not cause 
an open circuit condition within modules. 

The solar insolation loss caused by shadowing is the product of the 
available insolation to an unshadowed array times the fraction of the 
array in the shadow of the adjacent array to the south. So, available 
insolation energy loss is the time integral of the insolation intensity 
loss. 


B . THE MODEL 

To calculate the available insolation intensity (or energy) lost 
due to shadowing, the following parameters are required (see Figure 
3-1): the length of the array surface, the array tilt angle, the sun 
elevation angle, and the array spacing. 



n ARRAY SURFACE LENGTH d ARRAY SPACING 

T ARRAY TILT ANGLE < SOLAR ALTITUDE ANGLE 



*The front line of arrays will not receive any shadowing and must be 
considered separately. 


For a given distance, d, between arrays, it can be shown that the 
fraction of available energy lost during hour u is given by: 


F (w) 
s 


1 

1 - (cos 
0 


t ♦ d/a) 


sin g 

sin (f ♦ t ) 


if €<0 

if 0 < * < 6 (Equation 3-1) 
if f >6 


where F 8 (w) * shadow loss fraction during hour angle at 

T ■ array tilt angle (radians) 

d ■ array spacing (same units as a) 

a =• array surface length (same units as d) 

sin e “ sine of sun elevation angle (a function of u) 

■ cos <j> cos 6 cos co + sin $ sin 6 where <P f <*>, and 5 
are as defined in Section II (pages 2-4 and 2-8). 


Since the above values are hourly, in order to calculate the energy 
losses over time, hourly power values (assuming the power during the hour 
is constant) can be summed using: 


Available energy 
lost due to 
shadowing 


sunset 


CO 


sunrise 


I (w) 
Tt 


F (w) 
s 


since the power values from Section II are discrete.* This sum represents 
a daily insolation estimate. To compute energy losses on a monthly or 
annual basis, the daily value must be multiplied by the number of days in 
the particular month and the sum taken over months in the year. 


The program discussed in Section VI can be used to trade-off the 
energy losses due to shadowing with the total amount of space required 
for the system. The effects of shadowing become more pronounced in the 
higher latitudes so different spacings are required, depending on the 
economic parameters of the design. 


It should be noted that this model is not restricted to the insolation 
model of Section II and recorded data could be used. 


♦Similarly, if the available energy is to be calculated, (1 - F s (<*0) 
must be substituted for F s (<0 in the above expression. 
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SECTION IV 

REFLECTOR AUGMENTATION 


A. THE PROBLEM 

The question of reflector augmentation is concerned with a specific 
array design known as the tandem rack design (Reference 20). The geom- 
etry of this design is illustrated in Figure 4-1 and the derivation of 
the model for solar noon may be found in References 18 and 19. Again it 
is assumed that the arrays are south facing, but in this case, the north 
side of the array support structure is covered with a highly reflective 
surface (such as Mylar) facing the adjacent array. The insolation 
available to the array surface is determined by calculating an augmenta- 
tion factor. This augmentation factor represents the additional insola- 
tion from the reflective surface. 

This model is especially useful in conducting trade-offc among array 
spacing, tilt angles, and reflector angles. 


B. THE MODEL 

Examining Figure 4-1, it can be seen that when the sun elevation 
angle, € , is less than the reflector angle, 8 , there is no reflector 
augmentation, and in fact there are losses due to shadowing. In this 
case, the expression for hourly insolation is given by: 

I Tt (co) X (1 - F s (u>)) 

where I<pt(<*>) is the total insolation available to the surface during 
hour u , and F s (w) is the shadow loss fraction during hour ui (see 
Section III). 

However, when the sun elevation angle is greater than the reflector 
tilt, reflection occurs. If the augmentation of array illumination by 
reflected sunlight is given by A I/I, where AI is the additional insola- 
tion and I is the available insolation without reflectors, the multi- 
plier that yields the total insolation at the collector surface is given 
by 1 + A I/I. 

For a specular reflectance, p', the expression for the augmentation 
multiplier is given by: 


1 + AI/I = 1 + p' X sin ( e - 8 ) X (for 6 < <r < 2 6) 

where: 6 is the reflector tilt angle (= tan - ^ [a X (sin t)/d]); a is 
the array length in the meridian plane; R is the reflector length in the 
meridian plane; and L is the portion of the reflector length for which 
reflected insolation is not intercepted by the adjacent array for solar 
elevation angles greater than 28. 
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PARTIAL 

REFLECTION 

0<e<20 



FULL 

REFLECTION 


0<<< 20 



FULL 

REFLECTION 


€ >2 0 



Figure 4-1. Reflector Augmentation Geometry 



Since L * 0 for e<28, a simpler form can be used for sun elevation 
angles less than twice the reflector inclination angle: 


1 + AI/I » 1 + p' x sin (e - 8)x^£~- 

sin 6 


(for 8< e<28) 
(Equation 4-1) 


And for sun elevations greater than twice the reflector inclina- 
tion, the multiplier is corrected for reflector losses by: 

1 + AI/I - 

1 + p 1 X sintxsin (e - 8 ) X 


(for e>28) 
(Equation 4-2) 

So, for a given hourly insolation estimate from Section II or recorded 
data, multiplication by (1 + Al/l) yields the total available plus reflected 
insolation on the array surface. That is, 


CSC 8 - 


cot 8 + cot T 


sinS 


tan (e - 28) 


+ cos8 


I T1 >) X 

(1 - (Equation 3-1*)) 

if e< 8 

T xt >) x 

(Equation 4-1) 

if 8 < e<2 8 

T Tt (w) X 

(Equation 4-2) 

if e > 2 8 


Using the energy equation in Section III, this model can be used 
to determine the optimal spacing and reflector tilt to maximize energy 
output by multiplying the insolation level by either the shadow factor 
or the augmentation factor during hour to , depending on whether € is 
less than or greater than the reflector tilt, 8 . The computer program 
discussed in Section VI can be used to determine the optimal reflector 
angle and array separation distance for a particular site. 

Again, recorded data may be used with this model if desired. 


*Equation 3-1 is on page 3-2. 
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SECTION V 


DISCUSSION 


A number of models have been presented that enable the solar array 
system engineer/designer to examine insolation parameters related to the 
system design. These models are concerned solely with the available 
insolation that can be expected for a particular location and array 
spacing, with or without reflecting surfaces. The estimates are average 
values intended for use in engineering applications. The advantages of 
the approach in this document are that sensitivity analyses can easily 
be performed and, given the appropriate cost function for array spacing, 
an optimal array spacing can be calculated. 


A. INSOLATION MODEL SENSITIVITY 

Certain limitations should be considered in using this_model. The 
accuracy of the insolation model depends primarily on the values. 

If an hourly ftp profile for a typical day were available^ it would 
capture the shape of the profile more accurately. When ftp values 
greater than 0.75 occur, the tabulated correlations between ftp and 
K,j in Section II can introduce extrapolation error. Sunny locations 
can yield ftp values greater than 0.8, which requires extrapolation 
into the range where no correlated values exist. In the model and 
£rogram presented here, extrapolation is limited to the upper bound of 
ftp * 0 • 80 • 

In applications where load matching within the day is important, 
the insolation model presented here may be deficient. (Again, this 
depends on the randomness of the solar climate for the location.) 
Climates with highly variable cloud patterns will cause additional 
error; however, since the ftp values are derived from monthly average 
values, daily and monthly totals should be reliable. 

A sensitivity analysis was conducted in order_to determine the 
effect of changes in the solar climate parameter lCp on the total inso- 
lation I-p t . It was found that for a horizontal surface, there is 
essentially a one-to-one correspondence between percentage changes in 
ftp and the resulting changes in Lp t . However, as the array is 
tilted upward, the ftp value is more sensitive in winter months and 
less sensitive in the summer months; that is, the Kf parameter sensi- 
tivity depends on the relationship between the orientation of the array 
surface and the sun elevation angle. In the higher north latitudes, the 
magnitudes of the sensitivities increase due to the more dramatic 
effects of changes in the declination and sun elevation angles. 


B. INSOLATION MODEL VERIFICATION 

A weakness of the current model is the difficulty in evaluating its 
accuracy. With errors in data collection, comparisons are difficult to 
evaluate. The lack of detailed hourly insolation data for various tilt 
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angles at numerous sites is a further problem. Some work has been done 
in Reference 12 (pages 263-265); however, since clear-day profiles were 
used, this error analysis has limited applicability. 

Comparing the insolation model to three months of data collected at 
JPL (for a south-facing surface inclined at 15 degrees) yielded daily 
total estimates within 3% of the data values. These results are limited 
by the short record, but, as mentioned earlier, the model takes the 
daily amount of radiation and distributes it in a symmetric fashion 
about solar noon, thus introducing variation at the hourly level. The 
sum of the hourly absolute deviations of the model estimates from the 
data was much larger (6% to 8Z); therefore, this model is to be used 
c&utiously for applications where the profile within the day is impor- 
tant since it does not capture hourly trends. However, uncertainty 
about the trends in hourly data may exceed the error in using the model, 
depending on the micro-climate of the location. Improvement could be 
achieved using hourly Kf values if the data were available. 

The primary source of error is th£ data (total radiation on a hori- 
zontal surface) used to estimate the Kf parameters and which limit the 
accuracy of the model. Since the recorded data may have measurement 
inaccuracies, these errors translate into model inaccuracy. Because of 
this, small differences cannot be considered very significant (Reference 
21). Similarly, the extrapolation of this model to a site without 
recorded data must be done cautiously, since the model does not consider 
solar energy availability variations on a micro-climate scale and local 
topography might create local differences in solar insolation patterns. 

Two components are influenced by the sunrise hour and array orienta- 
tion. As can be seen from the definition of cos 0 t in Section II, the 
direct normal component is defined to be zero when the sun is behind the 
array. Thus, this component is unaffected. However, the diffuse and 
reflected components depend on the sunrise hour angle via the rj 
term. Since the diffuse and reflected components are small during the 
sunrise/sunset hours, the error introduced by using the horizontal sur- 
face sunrise angle is not believed to be significant. However, the 
effect of spectrum changes on morning and afternoon insolation values 
may complicate this problem. The actual sunrise time that the surface 
"sees" may be calculated by setting the expression for cos 0 t equal to 
zero and solving iteratively for w ■ W s . The author knows of no 
detailed study of the errors introduced by non-south-facing tilted sur- 
faces. Thus, the user of the SOLINS subroutine provided with this docu- 
ment should be aware of the above problems with non-south-facing tilted 
arrays. 

Finally, it should be noted that the model produces estimates of 
total available insolation and that photovoltaic cells respond only to 
certain wavelengths of the spectrum — a spectrum that changes during 
the day. Thus, actual available insolation will be less for solar cells 
since they respond only to a certain region of the spectrum (see Section 
II-D above). This problem is in a sense a conversion efficiency problem 
specific to the technology of photovoltaic conversion and was not 
addressed in this document. 
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C. SHADOW AND REFLECTOR MODELS 


The shadow loss and reflector augmentation models provide estimates 
of the insolation energy losses and gains associated with these 
problems. Obviously, there are a number of shortcomings. The accuracy 
of these models depends directly 'bn the accuracy of the insolation data 
input. The assumption that there is a one-to-one relationship between 
shadowing and insolation losses may be an optimistic estimate in the 
case where partial shadowing severely reduces module output. Another 
drawback of this model is that it is limited to rectangular-shaped 
collectors — non-uniform shadowing is not considered. The reflector 
augmentation model shares this problem. In addition, there may be 
double counting of reflected insolation in areas of low ground reflec- 
tance. The specular and ground reflectance must be chosen carefully; 
e.g., if there are no reflective objects (mountains, trees, etc.) facing 
the arrays, it may be sufficient to set the ground reflectance to zero 
and adjust the specular reflectance. This is an area requiring further 
examination. 

The remainder of this document presents a SOLINS computer program 
user's guide; a listing of SOLINS and related computer programs; and 
insolation profile data for 235 sites in the United States for various 
south-facing tilt angles. The computer program also computes shadow 
losses and reflector augmentation estimates on a monthly and annual 
basis. The values computed were verified by sampling 3% of the sites 
and reproducing all of the calculations using an APL computer program 
first and then by hand calculation. 
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SECTION VI 


SOLINS COMPUTER PROGRAM USER'S GUIDE 


The model presented here and the resulting computer program (called 
SOLINS) listed in Appendix A were written to support work being 
performed for the Low-Cost Solar Array Project at the Jet Propulsion 
Laboratory. In addition to calculating hourly and daily insolation 
estimates by month, the program calculates monthly and annual total 
available energy with array shadowing and/or reflector augmentation 
(that is, a reflecting surface facing the array surface). A copy of the 
program is available from the Computer Software Management and 
Information Center (COSMIC) operated for NASA by the Information 
Services Division of the University of Georgia Computer Center 
(Reference 22). 

Eight subprograms are contained in the SOLINS package. In addition, 
the user must provide a main program that will call SOLINS. Appendix A 
illustrates a sample main program and each of the eight subprograms. 

The example main program runs each of three options available with 
SOLINS and the resulting outputs are shown in Appendix A. Also included 
in Appendix A is the program that generated the m factor tables in 
Section II. 


A. SUBPROGRAM INS 

The program is based upon a modified version of the National Bureau 
of Standards model (see Section II) and calculates the direct, diffuse, 
and reflected components of total insolation on a tilted surface with a 
given azimuthal orientation (south facing, east facing, etc.). The 
subroutine that produces these estimates (called INS) generates a 
monthly, average day profile by hour of the insolation energy. The 
model uses a parameter that captures the monthly solar climate, so the 
estimates include clouding factors and other weather patterns — they 
are not clear day estimates. The user of the program specifies the 
latitude of the site, the array tilt, the azimuth, the month, the solar 
climate parameter, the ground reflectance, and the printout 
specification. 


B. SUBPROGRAMS SHADOW AND REEL 

In order to generate the shadowing and reflector estimates, the user 
must also specify the array spacing, the array surface length, and the 
spectral reflectance. Given the user's design data, the information 
required crrun the program is the site latitude and solar climate 
parameter. This will be discussed in further detail later in this 
Section. 
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C. COMPUTER CONFIGURATION REQUIRED — MEMORY AND TIMING 


The program was written in FORTRAN IV using the IBM 370/3032 at 
the California Institute of Technology. The program was designed to 
run in a batch mode and requires a maximum of 128 kilobytes of 
storage. The execution time is between 3 and 20 CPU seconds, 
depending on the specific user requirements. 


D. ACCURACY OF RESULTS 

Due to the uncertainties associated with insolation data (see 
Section V, page 5-2), all calculations ara performed in single 
precision. Thus, output tables have accuracy no greater than three 
significant figures. The program output was verified by hand 
calculation of a sampling of sites. 


E. CONTROL CARD REQUIREMENTS 

The deck setup shown in Figure 6-1 is required to run the 
program. First are the control cards, which consist of a job card 
with accounting information specific to the site installation and 
compiler instructions (cards for an IBM installation are given in 
Figure 6-1) which are instructions to the operating system. Next, the 
main program and the SOLINS package are followed by a control card 
(installation specific) and then any data that are to be read in by 
the main program. At the end of the deck is a card that indicates the 
end of the job. This card is also installation specific. The program 
given here has no data included with it. 
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F. USER INSTRUCTIONS 


The ueer hat three options available with the solar insolation pro- 
gram package (SOLINS): 

Option 1 — Calculate the insolation tables only. That is, for 
one or all months, calculate the average hourly and 
daily insolation for a tilted surface at a given 
aaimuthal orientation. 


Option 2 — Calculate the insolation at five different array 
spacings (one tilt angle only) and on a monthly or 
annual basis. The unobstructed insolation is also 
calculated and a fraction of total monthly and 
annual energy lost due to shadowing is computed. 

Option 3 — Calculate the reflector augmentation at a particular 
array spacing (one tilt angle only) and the fraction 
energy gain, by month or annually. 


Each of these options has a distinct set of requirements and a gen- 
eral set of input rules that apply to all three. 

The general input date requirements are that the variables ZKT (the 
12 solar climate values), XLAT (the latitude), TOWN (the location), and 
RH0(the ground reflectance) must be specified when SOLINS is called. 

XLAT is represented in the form degrees and minutes (e.g., 32.51 is 
equivalent to 32° 51'). 

In addition to the general input data requirements, the parameters 
of the subroutine call must be specified. The CALL statement is of the 
following form: 

CALL SOLINS (ISWTCH, MONTH, IPRINT) 

where ISWTCH ■ 1 represents Option 1 

■ 2 represents Option 2 

■ 3 represents Option 3 

MONTH ■ The calendar month (1-12) to be run or, if equal to 0, 
all 12 months are run automatically. This enables the 
user to change the parameters on a monthly basis to 
account for monthly or seasonal adjustment of the 
array tilt angle or monthly variation in ground 
reflectance to account for snow cover. 

IPRINT » 0 prints nothing but returns all daily values in the 
following variables: 

SUM(l) — daily total (unobstructed) insolation at 
tilt angle I (59). 

SUMS(1,J) — daily total with shadowing (at array 
spacing J) and one (1) tilt angle. 
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1 


SUMRE(l) — daily total with raflaetor augmentation 
at ona (1) tilt angla. 

■ 1 prints tha hourly and daily tablas onlyt 

■ 2 prints only tha shadow subrouting results; 

■ 3 prints the hourly and daily tablas 

associated with the INS subroutine (which 
is called by SHADOW) as well as the output 
from SHADOW; 

■ 4 prints the results of subroutine REFL only; 

■ 3 prints the hourly and daily tables 

associated with the INS subroutine (which 
is called by RBFL) as well as the output results 
from REFL. 

The tabulated output is in kilowatt-hours per square sseter. If 
other units are desired, the user oust redefine the variable RST in the 
INS subroutine from 1.377 kW/m^ to the equivalent value in the new 
units (e.g., 437 BTU/Oir-ft^)). Therefore, depending on the needs of 
the user, a variety of outputs are possible. 

To aid the user, an optional subroutine DATCK (ISHTCH) has been 
included. This subprogram should be called prior to calling SOLINS and 
it will perform a cursory examination of the input data by checking 
ranges of values. It should be noted that DATCK does not have to b« 
run for SOLINS to work; in fact, it may be removed from the package if 
desired. 

The only program restriction is that the calling program to SOLINS 
must contain the following labeled common statements 

COMMON/ SOL/ ZKT( 12), XLAT, RHO, TOWN(5), TILT(9), NO, NA, X(5), 
IUNIT, SREFL, BETA, ALNGTH, SUM(9), SUM(9,5), SUMRE(9), 
ZIT(24,9), SHD(24,1,5), REF(24,1) 

where these variable names are reserved for use by SOLINS. In addi- 
tion, the output unit must be specified using the integer variable 
IUNIT. This concludes the general requirements of the SOLINS package. 

Option 1 

The Option 1 (ISWTCH ■ 1) of the SOLINS subprogram computes the 
hourly and daily insolation estimates by calling the INS subprogram. 

The requirements for using this option ares 

1. Up to but not more than nine tilt angles may be run at one 
time. The tilt angles (in decimal degrees) are stored in 
TILT(I) and a zero tilt angle may be used. 

2. The number of tilt angles must be specified with the integer 
variable NA (< 9). 
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specified with the integer 

9 

6 ■ Northeast 

7 ■ Bast 

8 ■ Southeast 

9 ■ Horieontal 


orientation angle somewhere 
between the above choices, the variable WAZewould have to be 
redefined with that angle (0° ■ north, azimuth aeasured 
clockwise) . 

4. The ground reflectance oust be specified using the real 
variable RHO (a value of 0.2 is typically used). 

5. The general input requirements mist be satisfied. 

Option 2 

The Option 2 (ISWTCH ■ 2) of the SOLZNS subprogram computes the monthly 
total insolation estimates including shadowing by calling the sub- 
program, SHADOW, which in turn calls the INS subprogram. If the MONTH 
parameter is set equal to aero, the annual totals are also computed. 

The requirements for using this option aret 

1. All of the Option 1 restrictions must be satisfied except the 
number of tilt angles, NA, must equal one. 

2. The real variable ALNGTH, the array surface length (same units 
as X(O), must be specified (see Figure 3-1). 

3. The real vector X(«) of five array spacing values must be 
specified (same units as ALNCTH). 

4. The azimuthal orientation must be south facing (NO “ 1). 

5. The tilt angle should not equal zero (redundant to purpose of 
Option 2). 

Option 3 

The Option 3 (ISWTCH ■ 3) of the SOLINS subprogram computes the 
total insolation estimates including shadowing and reflector augmenta- 
tion by calling the subprogram RE7L, which in turn calls the INS sub- 
program. If MONTH ■ 0, the annual totals are also computed. The 
requirements for using this option are: 

1. All the Option 2 restrictions must be satisfied — except the 
arrav spacing value must be stored in X(l) (array spacing is 
not allowed to vary within a single call to REFL). 


3. The azimuthal orientation must be 
variable NO (1 & NO 19), where 

1 ■ South 

2 ■ Southwest 

3 ■ West 

4 ■ Northwest 

5 ■ North 

If the user requires an azimuthal 
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2. The real variable SREFL must contain the spectral reflectance 
value of the reflective surface. If no reflective surface is 
present, SREFL * 0.0. 

A sample main program with the corresponding output may be found in 
Appendix A. While this program package does not allow for every possi- 
ble user need, a listing of the programs is presented in Appendix A to 
enable the user to adapt the program to individual requirements. 
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APPENDIX A 


LISTING OF SOLINS AND RELATED COMPUTER PROGRAMS 


An example main program is given followed by the SOLINS, DATCK, 
SHADOW, REFL, and INS subroutines. Subroutine INS is fol lowed by its 
FUNCTION subprograms, ZD, SHDLOS, and AUG. The sample main program 
runs each of the SOLINS three options. 

The first option call in the main program generates the monthly 
tables of insolation values (kilowatt-hours/square meter) as shown on 
pages A-23 through A-26. 

The second option call in the main program generates the monthly 
insolation estimates and annual totals with and without array shad- 
owing for five array spacings. In addition, the fraction of energy 
lost due to shadowing is printed (page A-27). The monthly losses are 
computed by multiplying the daily amounts (with and without shadowing) 
by the number of days in the en month and computing a percentage 
loss. 


The third option call of the main program makes several passes 
through the reflector augmentation subroutine while varying the array 
spacing (pages A-28 through A-30). It can be seen that a maximum 
occurs at a spacing of 7 meters; however, a more detailed search 
(finer increments) would be required to determine the global optimum. 
Again, the daily amounts (before multiplication by the number of days 
in the month) are printed. 

Finally, the program used to generate Tables 2-2a and 2-2b on 
pages 2-10 and 2-11 is given on page A-31. 


A-l 



c 

C SAMPLE MAIN PKOGRAM 
C 

C THIS EXAMPLE RUNS THE INSOLATION PROGRAM INS, THE SHATOW PROGRAM, SHADOW, 

C AND THE REFLECTOR AUGMENTATION PROGRAM REFL . 

C 

COMMON /SOL/ZKT(12)»XLAT,RHO , TOWN (5),TILT(9),N0.NA,X(5)»IUNIT, 

1 SREFL»BETA,ALNGTH,SUM(9), SUMS ( 9 , 5 ) , SUMREC 9),ZIi!2A,9),REF(2A,l), 

2 SHD(2A ,1,5) 

C THE VARIABLE N RUNS OVER THE NUMBER OF SITES. 

N = 1 

DO 1 K=l,N 

C READ THE SITE LOCATION AND THE SOLAR CLIMATE VECTOR. 

READ( 5, 2 ) ( TOWNC I),I=1,5)»XLAT 
READ( 5 , A ) (ZKT(L),L=I,12) 

C SPECIFY THE NUMBER OF TILT ANGLES 
NA-9 

C SET THE TILT ANGLES. 

TILTC1 )=0.0 
TILT(2)=I0.0 
TILT! 3)=20 . 0 
TILT(A>=?5.0 
TILT(5)=30.0 
TILT(6)=35.0 
TILT(7)=A0 .9 
TILT(8>=50.0 
TILT (9) =6 0.0 

C SPECIFY THE AZIMUTHAL ORIENTATION. 

NO’S 

C SPECIFY THE OUTPUT UNIT 
IUNIT=6 

C PRINT SWITCH! SEE SUBROUTINE SOLINS). 

IPRINT=1 

C GROUND REFLECTANCE! MOUNTAINS, TREES, ETC. FACING THE ARRAY). 

RHO = . 2 

C ISWTCH DETERMINES WHICH OPTION IS TO BE RUN. 

ISWTCH=1 

C SET THE MONTH TO BE RUN. 

MONTH=0 

C CALL THE SOLINS ROUTINE FOR OPTION ONE. 

CALL SOL INS (ISWTCH, MONTH, IPRINT) 

35 CONTINUE 
C SET OPTION TWO. 

ISUTCH=2 
NO = l 

C SPECIFY THE LENGTH OF THE ARRAY FROM THE GROUND EDGE ALONG THE SURFACE 
C TO THE TOP EDGE(METERS). 

ALNGTH=A . 0 

C SPECIFY THE VARIOUS ARRAY SPACINGS (METERS ) --HORIZONTAL DISTANCE FROM TOP OF 
C FRONT ARRAY TO BOTTOM OF REAR ARRAY. 

X! 1 ) : 3 . 0 
X(2)=A.O 
X! 37=5.0 
X(4)=6 .0 
X(5)=8.0 

C NEW PRINT OPTION. 



* 



IPRINT=2 

C RUN ALL THE MONTHS. 

MONTH=0 

C SET A NEW TILT ANGLE. 

TILTC1)=2S.« 

C RESET THE NUMBER OF TILT ANGLES TO ONE. 

NA = 1 

C CHECK THE DATA FOR OPTION TWO. 

CALL DAT CK( ISWTCH ? 

C CALL SOL INS TO RUN OPTION TWO. 

CALL SOLINSC ISWTCH, MONTH, IPRINT) 

C TURN SWITCH TO OPTION THREE--THE REFLECTOR AUGMENTATION PROBLEM. 
ISWTCH=3 

C SPECIFY THE SPECTRAL REFL ECTANCEI IF REFLECTORS ARE PRESENT). 
SREFL=0.75 

C RUN THE ENTIRE YEAR. 

MCNTH=0 

C RESET THE PRINT OPTION. 

IPRINT =A 

C RESET THE ARRAY SPACING. 

X( 1 ) =3 . 

C CHECK THE DATA FOR OPTION THREE. 

CALL DATCK(ISWTCH) 

C CALL SOL INS FIVE TIMES USING OPTION THREE. 

DO 5 J=l,5 

CALL SOLINSI ISWTCH, MONTH, IPRINT) 

C INCREMENT THE ARRAY SPACING. 

X<l)=X(l)+«.0 
5 CONTINUE 

1 CONTINUE 

2 FORMAT ( 5 A4 > FI 0 . 0 ) 

A FORMAT ( 12F5 . 3 ) 

3 STOP 

END 
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SUBROUTINE SOLINS 
ALIFORNIA INSTITUTE OF TECHNOLOGY 

JET PROPULSION LABORATORY 

4800 OAK GROVE DRIVE 
PASADENA, CALIFORNIA 


JEFF H. SMITH 


THIS PROGRAM WAS WRITTEN TO SUPPORT WORK CURRENTLY BEING PERFORMED 
FOR THE LOW COST SOLAR ARRAY PROJECTCLSA) AT JPL. 

FURTHER DETAIL MAY BE FOUND IN REFERENCE I* SMITH, J.H., HANDBOOK 
OF SOLAR ENERGY DATA FOR SURFACES OF ARBITRARY ORIENTATION IN THE U.S. 
LSA DOCUMENT NUMBER 5101-91. JET PROPULSION LABORATORY; PASADENA, 
CALIFORNIA, JANUARY 1980. 


SUBROUTINE SOLINSt ISWTCH, MONTH, IPRINT ) 

THE SOLINS SUBPROGRAM CONSISTS OF SIX SUBPROGRAMS AND CALLS THE 
APPROPRIATE OPTIONS AVAILABLE DEPENDING ON THE VALUE OF ISWTCH. 

NOTE: PRIOR TO USING SOLINS, THE USER SHOULD CONSULT THE PROGRAM AND 

DATA REQUIREMENTS SECTION BELOW. 

END NOTE. 

THE THREE OPTIONS ARE: 


OPTION 1) ISWTCH=1 : GENERATES HOURLY INSOLATION VALUES FOR UP TO 9 

DIFFERENT TILT ANGLES AND ANY AZIMUTHAL ORIENTATION 
AND GROUND REFLECTANCE. 


OPTION 2) ISWTCH=2 : GENERATES MONTHLY TOTAL INSOLATION VALUES WITH AND 

WITHOUT SHADOWING FOR 5 DIFFERENT ARRAY SPACINGS AT 
ONE TILT ANGLE AND ORIENTATION PER MONTH. THIS OPTION 
ALSO CALCULATES THE TOTAL ANNUAL ENERGY WITH AND 
WITHOUT SHADOWING AND THE MONTHLY AND ANNUAL FRACTION 
OF ENERGY LOST AT EACH OF 5 ARRAY SPACINGS. 

CONTINUE 

OPTION 3) I5WTCH=3 •' GENERATES THE MONTHLY ENERGY WITH REFLECTOR AUGMEN- 
TATION AND THE FRACTION OF ENERGY GAINED. THE SHADOW 
FACTOR IS ALSO INCLUDED IN THIS OPTION AND THE ARRAY 
SPACING MUST BE SPECIFIED. 


THE MODELING APPROACH USED TO GENERATE THE INSOLATION TABLES IS DISCUSSED IN 
THE INS SUBROUTINE. 
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c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CONTINUE 

THE FOLLOWING VALUES OF IPRINT DETERMINE THE INTERMEDIATE AND FINAL OUTPUT 
OBTAINED. 


IPRINT=0 
IPRINT=1 
IPRINT =2 
IPRINT=3 
IPRINT=A 
IPRINT=5 


NO PRINTOUT 

HOURLY/DAILY (H/D) TABLE ONLY 
SHADOW TABLE ONLY 
" " AND H/D TABLES 

REFLECTOR TABLES ONLY 
" " AND H/D TABLES 


COMMON /SOL /ZKT (12), XL AT , RH0,T0WN( 5) »TILT(9),N0.NA,X(5),IUNIT, 

1 SREFL»BETA,ALNGTH,SUM(9)> SUMS ( 9 > 5) , SUMRE( 9),ZIT(2A,9),REF(2A»1)» 

2 SHD(24,1,5) 


THE FOLLOWING GO TO STATEMENT CALLS THE APPROPRIATE OPTION. 

GO TO (1,2,3), ISWTCH 

THE FOLLOWING VARIABLES ARE DEFINED FOR COMPATIBILITY PURPOSES ONLY- 
PROGRAMMING REQUIREMENT ONLY 

1 DO 10 K = 1 > 5 
10 X(K)=1.0 

BETA=1.0 

CALL INSCISWTCH, MONTH, IPRINT) 

RETURN 

2 CALL SHADOWCISWTCH, MONTH, IPRINT) 

RETURN 

3 CALL REFLCISWTCH, MONTH, IPRINT) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


REQUIREMENTS 

ASSOCIATED WITH EACH OPTION IS A SET OF RULES WHICH MUST BE ADHERED TO IN 
ORDER FOR THE PROGRAM TO EXECUTE SUCCESSFULLY. THESE REQUIREMENTS MUST BE 
SATISFIED PRIOR TO CALLING SOLINS. 

IN ADDITION TO THE OPTION REQUIREMENTS THERE ARE A NUMBER OF INPUT DATA AND 
PROGRAM REQUIREMENTS WHICH ARE ENUMERATED BELOW. 


INPUT DATA REQUIREMENTS 


CONTINUE 


SINCE THE SOLINS PROGRAM IS USED AS A SUBROUTINE, THE ONLY DATA 
REQUIREMENTS ARE THAT THE VARIABLES ZKT , XL AT , AND TOWN BE SPECIFIED 
WHEN SOLINS ( OR DAT CK IF USED) IS CALLED. ZKT SHOULD CONTAIN 12 SOLAR 
CLIMATE VALUES FOR J ANU ARY( ZKT ( 1 ) ) THROUGH DECEMBERC ZKT ( 12 ) ) . 

TYPICAL VALUES MAY BE FOUND IN THE REF. DOC. TOWN CONTAINS THE SITE 
NAMECUP TO 20 CHARACTERS) WITH FORMAT ( 5A4 ) . XL AT IS SPECIFIED IN 
DEGREES AND MINUTES SEPARATED BY A DECIMAL POINT ( DEG . MIN ) 



C 


CONTINUE 
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PROGRAM REQUIREMENTS 




THE PROGRAM DATCK DOES NOT HAVE TO BE RUN PRIOR TO CALLING SOLINS. 

THE MAIN PROGRAM MUST CONTAIN THE FOLLOWING COMMON STATEMENT! 

COMMON /$0L/ZKT(12) »XLAT »RHO. TOMN(5) .TILT(9),N0,NA.X(5).IUNIT, 

* SREFL,BETA,ALNGTH,SUM(9),SUMS(9,5>,SUMRE(9>,ZIT(2A.9),REF(2*,1>. 

« SHD(24«1>5) 

CONT IHUE 

THE UNITS OF THE TABLES MAY BE CONVERTED BY REPLACING THE VARIABLE RST 
(THE SOLAR CONSTANT) WITH ITS EQUIVALENT VALUE IN THE DESIRED UNITS. 

RST APPEARS IN THE BEGINNING OF THE INS SUBROUTINE. 

N 0 T t- THE VARIABLES USED IN THE COMMON STATEMENT ARE RESERVED FOR 
SOLINS AND ITS SUBROUTINES. 

THE ARRAY SUM(I) WILL CONTAIN THE DAILY TOTAL INSOLATION FOR ARRAY 
TILT ANGLE I FOR MONTH (NO SHADOWING OR REFLECTOR AUGMENTATION). 

THE ARRAY SUMS(I,J> WILL CONTAIN THE DAILY TOTAL INSOLATION 
(WITH SHADOWING) FOR ARRAY SPACING J AND ONE ARRAY TILT ANGLE FOR 
MONTH. 

THE VARIABLE SUMRE(l) WILL CONTAIN THE DAILY TOTAL INSOLATION WITH 
REFLECTOR AUGMENTATION AND SHADOWING FOR ONE ARRAY TILT ANGLE AND 
MONTH. 

END NOTE. 


OPTION 1 REQUIREMENTS 


CONTINUE 


UP TO BUT NOT MORE THAN 9 TILT ANGLES MAY BE RUN AT ONE TIME(STORED IN THE 

ARRAY NAMED TILT). A ZERO TILT ANGLE MAY BE USED. 

THE NUMBER OF ANGLES MUST BE SPECIFIED WITH THE INTEGER VARIABLE NA. 

THE AZIMUTHAL ORIENTATION MUST BE SPECIFIED WITH THE INTEGER VARIABLE NO: 
WHERE NO ORIENTATION NUMBER AND 

1=S0UTH 

2=S0UTHWEST 

3=WEST 

A=NORTHWEST 

5=NORTH 

6=NORTHEAST 

7=EAST 

8=S0UTHEAST 

9=H0RIZ. 

THE OUTPUT UNIT MUST BE SPECIFIED USING THE INTEGER VARIABLE IUNIT. 

THE GROUND REFLECTANCE MUST BE SPECIFIED USING THE REAL VARIABLE RHO. 

THE REAL VECTOR ZKT MUST BE SPECIFIED. 

THE TOWN VECTOR SHOULD BE SPECIFIED. 

THE LATITUDE XLAT MUST BE SPECIFIED. 

THE INTEGER VARIABLE MONTH MUST BE SPECIFIED: 

MONTH DETERMINES WHETHER ALL 12 

MONTHLY TABLES ARE COMPUTED 
(=-l), OR ONLY ONE MONTHS 
TABLES ARE COMPUTEDC E. G. 

M0NTH=3 FOR MARCH). 


ooooooo oooooooo 


c 


OPTION 2 REQUIREMENTS 


CONTINUE 


THE OPTION 1 REQUIREMENTS MUST BE SATISFIED EXCEPT THE NUMBER OF TILT 
ANGLES NA MUST EQUAL ONE. 

THE REAL VECTOR Or ARRAY SPACING VALUES MUST BE SPECIFIED. 

THE AZIMUTHAL ORIENTATION MUST BE SOUTH FACING(NO=l ) . 

THE TILT ANGLE SHOULD NOT EQUAL ZERO(HORIZONTAL) . 


OPTION 3 REQUIREMENTS 


CONTINUE 


THE OPTION TWO REQUIREMENTS MUST BE SATISFIED EXCEPT THE ARRAY SPACING 
VALUE MUST BE STORED IN XC 1 5 . 

THE REAL VARIABLE SREFL MUST CONTAIN THE SPECTRAL REFLECTANCE VALUE OF THE 
REFLECTOR. 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 1 

RETURN J 

END 


j 

t 
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SUBROUTINE DATCK 


JET PROPULSION LABORATORY— PASADENA, CALIFORNIA 
SUBROUTINE DATCKC ISWTCH) 

THIS SUBROUTINE CHECKS THE APPROPRIATE DATA REQUIRED FOR THE PARTICULAR 
ISWTCH OPTION. 

NOTE: THIS SUBROUTINE IS NOT REQUIRED TO RUN THE SOLINS PROGRAM 

AND MAY BE OMITTED AT THE USERS DISCRETION. 

COMMON /SOL/ZKTC 12) ,XLAT , RHO » TOWNC 5 ) , TILTC 9 ) ,NO,NA ,XC 5) , IUNIT, 

1 SREFL .BETA , ALNGTH, SUMC 9) , SUMS ( 9, 5) , SUMREC 9),ZITC29,9),REFC29,1), 

2 SHDC29,1,5) 

C THIS SUBROUTINE CHECKS THE INPUT DATA TO SUBROUTINE INS. 

C A A 

C THE FOLLOWING VALUES ARE ALWAYS CHECKED. 

DO 10 1=1,12 

IFCZKT ( I ) . LT . 0 .OR .ZKT(I ) . GT . 1 . 0) GO TO 100 
10 CONTINUE 

IFCXLAT.LT.0.OR.XLAT.6E.90) 60 TO 200 
IFCRHO.LT.O.OR.RHO.GT.l.O) GO TO 300 
DO 20 I S 1 »NA 

IF(TILT(I).LT.0.OR.TILTCI).GT.90.) GO TO 900 
20 CONTINUE 

IFCNO.LT.l.OR.NO.GT.S) GO TO 500 
IFCNA.LT.1.0R.NA.GT.9) GO TO 600 
GO TO Cl. 2, 3), ISWTCH 


1 RETURN 

C B B 


C THE FOLLOWING DATA CHECKS ARE PERTINENT TO THE OPTION TWO CALLS. 

2 DO 30 1=1,5 
T=X(I)/ALNGTH 

IFCT.LT.O. .OR.T.GT.20. ) GO TO 700 
30 CONTINUE 

IFCNO.NE.l) GO TO 160 
162 CONTINUE 

IFCNA.GT.l) GO TO 170 
172 CONTINUE 

IFCTILT<l).EQ.O.) GO TO 180 
RETURN 

180 WRITEC IUNIT , 181 ) TILTC1) 

181 FORMATCIOX. ’THE TILT ANGLE IS ZERO— THIS WILL YIELD THE SAME RESUL 
IT AS INS— WILL PROCEED WITH INS.') 

ISWTCH=1 

RETURN 

160 WRITEC IUNIT , 161 ) NO 

161 FORMAT C 10X, ’ORIENTATION IS NOT SOUTH FACING— ERRORS MAY BE INTRODU 
ICED’ ,19) 

GO TO 162 

170 WRITEC IUNIT , 171 ) NA 

171 FORMATCIOX, ’NA MUST EQUAL ONE— RESET BY PROGRAM’, 19) 

NA = 1 

GO TO 172 
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C B - - - B 

C 

C C - C 


C THE FOLLOWING DATA CHECKS ARE PERTINENT TO THE OPTION THREE CALLS. 
3 T=XC 1 )✓ ALNGTH 

BETA= ATANCSINCTILTC1)*3.14159/18Q.)/T) 

TB=BETA*180./3. 14159 

IFITB.LT. 0. .OR.TB.GT.90.) GO TO 800 

RETURN 


C C C 

C 

C D - D 


C THE FOLLOWING ARE THE PRINT STATEMENTS FOR THE ABOVE CHECK PROCEDURES. 


100 

WRITEdUNIT 

.101) 

C2KTCI), 

1 = 1,12) 

101 

FORMAT C 10X, 
GO TO 10 

•DATA 

ERROR 

IN 

ZKTM2F8.3) 

200 

WRITEdUNIT 

.201) 

XLAT 



201 

FORMATCIOX, 

RETURN 

•DATA 

ERROR 

IN 

XLAT* ,F8.2) 

300 

WRITEdUNIT 

,301) 

RHO 



301 

FORMATCIOX, 

RETURN 

•DATA 

ERROR 

IN 

RHO' ,F8.2) 

400 

WRITEdUNIT 

,401) 

TILTCK) 


401 

FORMATCIOX, 
GO TO 20 

•DATA 

ERROR 

IN 

TILT VALUE’, F5. 2) 

500 

WRITEdUNIT 

,501) 

NO 



501 

FORMATCIOX, 

RETURN 

•DATA 

ERROR 

IN 

NO* ,14) 

600 

WRITEdUNIT 

,601) 

NA 



601 

FORMATCIOX, 

RETURN 

'DATA 

ERROR 

IN 

NA' ,14) 

700 

WRITEdUNIT 

.701) 

XCI) 



701 

FORMATCIOX, 
GO TO 30 

'DATA 

ERROR 

IN 

X' ,F3.3) 

800 

WRITECIUNIT 

,801) 

TB 



801 

FORMATCIOX, 
GO TO 2 

•DATA 

ERROR 

IN 

BETA' , F8 . 3) 

END 


ooooo 


SUBROUTINE SHADOW 
JET PROPULSION LABORATORY— PASADENA, CALIFORNIA 


SUBROUTINE SHADOUC ISWTCH ,KMONT , IPRINT ) 

C THIS SUBROUTINE CALCULATES THE MONTHLY AND ANNUAL SOLAR INSOLATION ENERGY 
C WITH AND WITHOUT ARRAY SHADOWING. THE MODEL ASSUMES SOUTH FACING ARRAYS 
C AND UNIFORM INSOLATION LOSS WHEN SHADOWED( I . E. 50 PERCENT OF THE ARRAY 
C IN SHADOW IMPLIES 50 PERCENT AVAILABLE INSOLATION LOSS). THE EAST-WEST 
C SHADOWING IS ASSUMED NEGLIGIBLE. SEE REFERENCE DCCUMENTCCITED IN 
r SOL INS DOCUMENTATION). 

C 

COMMON /SOL/ZKTC12),XLAT,RHO.T0WNC5),TILT<9) ,N0,NA,X(5) .IUNIT, 

1 SREFL,BETA,ALNGTH,SUM(9> SUMS(9,5) ,SUMRE( 9) ,ZIT(29 , 9) ,REF(29, 1 ) , 

2 SHDC2A,1,5) 

DIMENSION LDAYC 12) ,S(5),SS(5), FRACM( 12, 5) , FRACY( 5 ) , ORIENTL 3, 9) 

REAL LAT 

DATA PI/3 . 1A15927/ 

DATA ORIENT/' SOU'.'TH ',' ’,' SOU • . ' THWE* , ' ST ',' WES', 

l'T ’,' ',' NOR' , 'THWE’ , 'ST ',' NOR'.’TH ',’ ',' NOR’, 

7 ' THEA ' , ' ST ',' EAS'.'T SOU ' . ' THEA • , ’ ST ',' HOR’. 

3'IZ. '/ 

DATA L DAY/31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 

1 30.31/ 

M0NTM=KM0NT 

XL AX=AINT (XLAT ) 

LAT=(XLAX< <XLAT-XLAX)/.6)*PI/180 . 

LAX=INT (XLAT) 

MINUTE=(XLAT-XLAX+.0001)*100 
I F( IPRINT . EQ . 0 ) GO TO 32 
WRITE(IUNIT.IO) ( TOWN (I), I =1,5), LAX, MINUTE 
10 FORMAT(1H1,30X, 'LOCATION: ', 5A9 , 1 OX, ' LATITUDE •• ',19, 

1 ' DEGREES' ,13, ' MINUTES’,///) 

WRITE(IUNIT,15) 

15 FORMATUX, 'MONTH' . IX. 'TILT', 2X, 'ORIENTATION', 2X,' GND. REFL. '. 

1 3X, 'UNOBSTRUCTED INSOL ATION’ , 10X, ' INSOLATION W/SHADOW' , 

1 9X, 'FRACTION LOST' ,/, IX, 123( • - ' ) ) 

WRITEC IUNIT ,18) 

18 F0RMATC83X, 'ARRAY SPACINGCMETERS ) • ) 

WRITE( IUNIT, 17) (XCM) ,M=1 ,5) , (X(M) ,M=1 ,5) 

17 F0RMAT(68X,2(5F5.2,5X),/,68X,56( '-'),/) 

32 CONTINUE 
TOT=0. 

DO 20 K = 1 ,5 

C THE VARIABLE SSfl) KEEPS THE RUNNING SUM FOR THE ANNUAL AMOUNTS. 

SS ( K ) = 0 . 

20 CONTINUE 

IFCMONTH.NE.O) IST0P=1 
IF(MONTH.NE.O) GO TO 13 
MONTH =0 

25 IFTMONTH. EQ. 12) GO TO 30 
M0NTH=M0NTH+1 
13 CONTINUE 

C CALCULATE THE INSOLATION VALUES FOR MONTH 
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CALL INSUSUTCH, MONTH, 1PRINT) 

C 

C CALCULATE THE MONTHLY INSOLATION 
C=SUM( 1 )»L DAY (MONTH) 

DO 3 1-1,5 
3 S<I)-0. 

C CALCULATE THE ANNUAL INSOLATION THUS FAR( IF REQUIRED) 

DO I 1=1,5 

S(I)=$(I ) +SUMS( 1 , I )*L DAY (MONTH) 

SS( I ) = S ( I ) + SS ( I ) 

C CALCULATE THE MONTHLY FRACTION OF ENERGY LOST DUE TO SHADOWING. 

1 FRA CM (MONTH, I ) = ABS( 1 .-SCI )/C) 

TOT=TOT+C 

IF(IPRINT.EQ.O) GO TO 33 
C OUTPUT THE MONTHLY RESULTS. 

WRITE(IUNIT,6)M0NTH.TILT(1),(0RIENT( J,NO) , J=1 , 3 ) , RHO , SUMC 1 ) , 
1 (SUMS(1.K),K=1,5),(FRACM(M0NTH,I),I=1,5) 

6 F0RMAT(1X»I3,2X»F5.2»AX,3AA,3X,F4.2,15X,F4.2.15X, 

1 2C5F5.2.5X),/) 

33 CONTINUE 

I F( I STOP . EQ . 1 ) GO TO 29 

2 GO TO 25 
30 CONTINUE 

IF(ISTOP.EQ.l) GO TO 29 
DO 5 1=1,5 

C CALCULATE THE ANNUAL FRACTION OF ENERGY LOST. 

5 FRACYCL)=l-5S(L)/TOT 
IF(IPRINT.EQ.O) GO TO 29 
C OUTPUT THE ANNUAL RESULTS. 

U'RITE( iUNIT , 7 ) TOT , ( SS< L ) , L = 1 , 5 ) , U RACYC L ) , L = 1 • 5 ) 

7 FORMAT (IX, 123(’ SOX, 'ANNUAL AMOUNTS s ' , FI 0 . 2 , 1 OX, 5FI 0 . 2 , 
12AX. 'ANNUAL FRACTION LOST : ' ,20X,5F10 .2) 

29 RETURN 
END 
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SUBROUTINE REFL 
JET PROPULSION LABORATORY— PASADENA, CALIFORNIA 
SUBROUTINE REFLCISWTCH, MONTH, IPRINT) 

THIS SUBROUTINE CALCULATES THE MONTHLY ENERGY AUGMENTATION WITH REFLECTORS. 
THIS MODEL IS SPECIFIC TO THE TANDEM RACK ARRAY DESIGN(SEE REF. DOC.). 

COMMON /SOL/ZKT (12),XLAT,RH0, TOWN ( 5 )>TILT!9),N0,NA,X!5),IUNIT> 

1 SREFL , BETA, ALNGTH, SUM( 9 ) , SUMS( 9,5), SUMRE( 9),ZIT!2A»9)*REF!2A,1)» 

2 SHD! 2A ,1,5) 

DIMENSION LDAYU2) 

REAL MPG.LAT 
DATA PI/3 .1415927/ 

DATA LDAY/31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31/ 

XLAX=AINT!XLAT) 

LAT=CXLAX+CXLAT-XLAX)/.6)XPI/180. 

LAX=INT(XLAX) 

MINUTE'- CXLAT-XLAX+. 0001 )«100 
IF(IPRINT.EQ.O) GO TO AQ 
C PRINT HEADINGS. 

WRITE! IUNIT, 10) ( TOWN! I ) , 1 = 1 , 5) , LAX, MINUTE, X( 1 ) 

10 FORMAT(1H1,///,30X, ’LOCATION: ', 5AA , 1 OX, ' LATITUDE: ’, 13, ' DEGREES’, 

113, IX, 'MINUTES' ,//,50X, ’ARRAY SPACING = ',F6.2,' METERS',//) 

WRITE! IUNI T ,11) 

11 FORMAT! AX, 'MONTH ', AX, ' ARR . TILT' , AX, ’REFL . TILT' , AX, 'SPEC. REFL.', 

1 AX , ' GND REFL . ’ , AX, ' INSOL . W/OUT SHADOW, REFL .' ,AX, ' INSOL . W/REFL. A 
2UG . ' , AX, ' FRAC . GAIN ’,/, 132 !'-•) ) 

AO CONTINUE 
DO 88 J=2.5 
88 X! J ) = 1 . 0 

IF(MONTH.EQ.O) GO TO 23 
ISTOP=l 
KMONTH=MONTH 
GO TO 15 
23 KMONTH=0 
TOT = 0 . 0 
TRE=0.0 

25 IF! KMONTH . EQ . 12) GO TO 30 
KMONTH=KMONTH+l 
15 CONTINUE 
C 

CALL INS!ISWTCH, KMONTH, IPRINT) 

C 

C CALCULATE THE FRACTION GAIN WITH REFLECTORS. 

MPG=SUNRE!1)/SUM! 1)-1.0 

IF(MONTH.EQ.O) TOT=TOT+SUM! 1 ) *LDAY! KMONTH ) 

IF! MONTH . EQ . 0 ) TRE=TRE+SUMRE! 1 )*LDAY! KMONTH) 

IF! IPRINT . EQ . 0 ) GO TO 50 
RT=BETA*180./PI 
C OUTPUT RESULTS. 

WRITE! IUNIT , 2 ) KMONTH , TILT ! 1 ), RT , SREFL , RHO , SUM! 1 ) , SUMRE! 1 ) , MPG 

2 FORMAT (5X,Ic , 8X, F5 . 2 ,8X, F5 . 2, 9X, FA . 2 , 10X, FA . 2, 12X, F8 . 2 , 17X, F8 . 2 , 

1 17X , F8 . 2 ) 
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50 CONTINUE 

IF(ISTOP.EQ.l) GO TO 13 
GO TO 25 
30 CONTINUE 

I F C IPRINT . EQ . 0 ) GO TO 13 
IF(MONTH.NE.O) GO TO 13 
ANF=TRE/T0T-1 .0 
WRITE! IUNIT # 62 ) TOT , TRE f ANF 

62 FORMAT! IX, 131! ' ) .//,56X, 'ANNUAL AMOUNTS'- • , F10 .2, 15X, F10 .2 , 15X, 
1 F10 . 2) 

13 RETURN 
END 
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CALIFORNIA 


SUBROUTINE INS 

JET PROPULSION LABORATORY— PASADEN . 

SUBROUTINE INS( ISUTCH .MONTH , IPRINT ) 

THIS PROGRAM WILL CALCULATE THE HOURLY SOLAR INSOLATION 
VALUES FOR A GIVEN SITE. THE USER MAY CONTROL THE TILT 
ANGLE OF THE SURFACE, THE GROUND REFLECTIVITY, THE 
ORIENTATION, AND THE MONTH. THE RESULTS REPRESENT 
AVERAGE DAILY INSOLATION VALUES AND THE METHOD OF COM- 
PUTATION IS BASED ON AN AGGREGATION OF DIRECT, DIFFUSE, AND REFLECTED 
INSOLATION COMPONENTS . 

FOR A CLEARER UNDERSTANDING OF THE MODEL USED IN THE PROGRAM, 

THE USER IS REFERRED TO THE REFERENCE DOCUMENT. 

CONTINUE 

COMMENTS INSERTED THROUGHOUT THE PROGRAM REFER TO THE ABOVE DOCUMENT 
FOR EXPLANATION. 


TO USE THE PROGRAM THE USER ISSUES THE FOLLOWING CALL: 

CALL INSCISWTCH, MONTH, IPRINT) 

WHERE: ISWTCH=1 COMPUTES AVAILABLE INSOLATION 

ISWTCH=2 " " " WITH SHADOW LOSSES 

ISWTCH=3 " " " WITH REFLECTOR 

AUGMENTATION AND SHADOWING 

MONTH DETERMINES WHICH MONTH IS RUN-IF EQUAL TO ZERO THE 
ENTIRE YEAR IS RUN 

IPRINT EQUAL TO ZERO YIELDS NO PRINTOUT 

IPRINT EQUAL TO ONE, THREE, OR FIVE YIELDS THE HOURLY AND 
DAILY TOTALS. 

CONTINUE 


DEFINE THE ARRAYS AND INITIALIZE THE DATA STATEMENTS. 

COMMON /S0L/ZKT(12) , XL AT , RHO , TOWN C 5 ) , TI LT ( 9 ) , NO , NA ,XC 5 ) , IUNIT , 

1 SREFL»BETA,ALNGTH,SUM(9), SUMS (9,5), SUMRE( 9),ZIT(24,9),REF(24, 1), 

2 SHDC 24 ,1,5) 

DIMENSION LDAYC12) .XDECC12) ,R(12) ,WAZ(12) , ORIENT ( 3 , 9 ) , 7ILDC 9 ) , 
1DNIC24) , ASH 24, 9) , RSI (24, 9) ,SUMD(9) ,SUMR(9) ,B(12) 

REAL LAT 

DATA LDAY /31 , 28 , 31 , 30 , 31 , 30 , 31 , 31 , 30 , 

1 31,30,31/ 

DATA XDEC /-19. 51, -10. 28, .2, 11. 56, 20. 14, 23. 27, 

1 20.26, 12.03, .37,-10.47,-19.58,-23.27/ 

DATA R /l. 03, 1.0207, 1.0057, .9875, .9727, .967, 

1 .9692, .9785, .9945,1.0133,1.0267,1.0327/ 

DATA WAZ / 18 n. , 22 * . ,270 . ,315. ,0. ,45. ,90. ,135./ 

DATA ORIENT/' SOU'.’TH SOU ' , ’ THWE’ , ' ST WES', 

l'T ’,' ',' NOR' , 'THWE' , 'ST ',' NOR'.'TH ',' ’,' NOR', 
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l’T ’,’ »,' NOR' , ’THWE’ , ’ST NOR'.'TH ’,’ NOR’, 

2 ' THEA ' , ' ST ’,’ EAS' , ' T ’,’ SOU ' , • THEA ’ , ’ ST HOR’, 

3' 12. ’/ 

DATA PI/3. 1415927/ 

DATA B /. 142, .144, .156, .18, .196. .205. .207, 

1 .201. .177, .16, .149, .142/ 


DEFINE THE UNITS OF THE TABLES BY SPECIFYING THE SOLAR CONSTANT. IN THIS 
CASE KW/M*«2 IS USED. 

RST=1 . 377 

DO LOOP 876 CONVERTS THE INPUT TILT ANGLES INTO RADIANS. 

DO 876 LL=1.NA 
TILDCLL )=TILT(LL5 
876 TILTCLL)=TILTCLL)XPI/180. 

C IF M=1 , THE COMPONENT TABLES OF THE TOTALS ARE PRINTED. IF M=0, THEY ARE NOT 
M=0 

C THE VARIABLE IUNIT DETERMINES WHICH OUTPUT UNIT THE OUTPUT 
C WILL BE PRINTED ON. 

C A2M CONTAINS THE CONVERTED VALUE(RADIANS) OF THE AZIMUTH 
C ANGLE. 

AZM=WAZ C NO )*P 1/180 . 

C THE FOLLOWING FOUR STATEMENTS CONVERT THE LATITUDE TO RAD- 
C IANS. 

XL AX-AINT CXLAT ) 

LAT=(XLAX+(XLAT-XLAX)/0.6)*PI/180. 

LAX=INT (XLAT) 

MINUTE=CXLAT-XLAX+ .0001)#100 
C CONSISTENCY CHECK. 

IFCNA .NE.l .AND.ISWTCH.GT.l) WRITE! IUNIT ,345) 

345 FORMAT ( 3X , ' XXXXXXXXXX' . 1 OX , ' THE NUMBER OF TILT ANGLES IS’, 

1’ GREATER THAN ONE SO MATRIX SHD MUST BE REDEFINED (24,NA,5)') 

30 FORMAT!// ,29X,' LOCATION ’,5A4,’ L ATITUDE’ , 14 , ’ DEGREES’, 

1 14,’ MINUTES') 

N-MONTH 

IFCMONTH.LT. 1) N=1 
134 CONTINUE 

RD = AINT CXDEC(N) ) 

C CONVERT THE DECLINATION TO RADIANS. 

DEC=CRD+CXDECCN)-RD)/0.6)*PI/180 . 

C CALCULATE THE SUNRISE HOUR ANGLE. 

COSWS=C-SINCLAT)/COSCLAT) )*CSINCDEC)/COSCDEC) ) 

IFCC0SWS.GT.1. .OR.COSWS.LT.-l) RETURN 
IFCIFRINT.EQ.O.OR.ISWTCH.'lT.l) GO TO 133 

IFCCN.EQ.l . OR. MONTH. GT.O) . AND . ( IFRINT . EO . 1 . OR . IPRINT . EQ . 3 . OR . IPRIN 
1T.EQ.5)) WRITECIUNIT, 41) 

I FC C N. EQ.l. OR. MONTH. GT.O) .AND. C IPRINT. EQ.l. OR. IPRINT. EQ. 3. OR. IPRIN 
1T.EQ.5)) WRITECIUNIT, 30) C T OWN C JJ ), JJ =1 , 5 ), LAX , MINUTE 
133 CONTINUE 

WS=ARCOSCCOSW$) 

TWS=US*12/PI 

SUNRIZ=12.-ABSCTWS) 

SUNSET=12.+ABSCTWS) 

COSLD=COSCLAT)xCOSCDEC) 


A-15 


rif ~ ii 








o o 


C A 


SINLD=SINCLAT)*SIN(DEC> 


A 


S = 0. 

C DO LOOP 500 CALCULATES THE VALUE OF THE TOTAL DAILY DIRECT NORMAL 
C RADIATION. 

DO 500 L=l# 39 
WW=WS*L/40. 

CZE=COSLD*COS(WW)+SINLD 

PAR=-B(H)/CZE 

APA=*ABS(PAR) 

IFCAPA.GT.80.) 60 TO 501 
ANS=EXP(PAR)*CZE 
GO TO 502 

501 ANS=0. 

502 S=ANS+S 
500 CONTINUE 

C A - A 

ANSO=EXP( -B( N)/( COSLD+SINLD) )*(C05LD+SINLD)/2. 

AI=WS/A0 . *( ANSO+S) 

C TAKE THE PROPORTION OF TOTAL DAILY RADIATION THAT REACHES THE EARTH AND 
C CALCULATE THE FRACTION THAT IS DIFFUSE RADIATION, KD. 

C THE ZKT VALUES REPRESENT THE MONTH-SPECIFIC AVERAGE DAILY FRACTION OF TOTAL 
C INSOLATION REACHING THE EARTHS SURFACE AFTER PASSING THROUGH THE ATMOSPHERE. 

ZKD=ZDCZKT(Nn 

C 

C CALCULATE THE EXTRATERRESTRIAL SOLAR RADIATION. 
H0-2A./PI*RCN)#RST*CC0SLD*SIN(WS)+WS*SINLD) 

PH^HO 

C 

THH=HO*ZKTCfO 
DHH=HQ#(ZKT (N)-ZKD) 

C 

C CALCULATE THE MONTHLY AVERAGE DAILY DIFFUSE RADIATION ON A HORIZ- 
C ONTAL SURFACE. 

RHH=HO*ZKD 

CALCULATE THE ABAR FACTOR. 

A=DHH/(24./PI*AI) 

FAC=A/ZKT(N) 

C DO LOOP 2 SETS THE ARRAYS WHICH WILL CONTAIN THE HOURLY VALUES 
C EQUAL TO ZERO. 

£ g _ ___ _ _______ ______ _______ _________ 

C ZERO ARRAYS. 

DO 2 K = 1 , 9 
DO 2 1=1 ,2^ 

DNI ( I ) = 0 . 

ASI(I,K)=0. 

RSICI,K)=0 . 

2 ZIT ( I , K ) =0 . 


C B B 

C 

C C C 


C DO LOOP 3 CALCULATES THE AVERAGE HOURLY SOLAR INSOLATION VALUES ON 
C A TILTED SURFACE. 

DO 3 1=1,24 
TIME=I-1. 
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WT=ABS(12.-TIME) 

IFCTIME-SUNRIZ) 3, 3, A 

4 IF( TIME-SUNSET ) 5,3,3 

5 W=CTIME-12.)*PI/12. 

COSZ-SIHLD+COSLD*COS(W) 

CO$k' : COS( DEC)*SIN(W) 

DO 10 K=1,NA 

25 CONTINUE 

C COSTH IS THE ANGLE OF INCIDENCE BETWEEN THE SUNS DIRECT BEAM AND A 
C NORMAL TO THE ARRAY SURFACE. 

C 

COSTH=(SIN(LAT)*COS(TILTCK))+COSCLAT)*SIN(TILT(K))* 

1 COS(AZM) )xSIN(DEC)+(COS(LAT)*COS(TILT(K) )-SIN(LAT)K 

1 SIN(TILT(K)>XCOS<AZM>)XCOS(DEC)*COSCW)- 

1 SIN(TILT(K))XSIN(AZM)XC0SW 
C 

IF( COSTH. LT.1.0E-25) COSTH=0.0 

C RR CONVERTS THE DAILY TOTAL TO AN HOURLY VALUECREFERENCE, STEP*). 
RR=PI/24.#(C0S(W)-C0S(WS) )/( SINCWS )-WS*COS(WS ) ) 

IFCRR.LT. 1.0E-25) RR=0. 

C CALCULATE THE AVERAGE HOURLY DIFFUSE RADIATION ON A HORIZONTAL 
C SURFACE. 

RHI =RHH*RR 
PAR2=-BCN)/C0SZ 
AP3=ABS(PAR2) 

IFCAP3.GT.80) GO TO 3 
DNI ( I ) = AxEXPC PAR2 ) 

IFCDNICD.LE.O.) DNICI)=0. 

DHI=DNI ( I ) XCOSZ 
IFCDHI.LE.O.) DHI =0 . 

C IF ARRAY IS HORIZONTAL WE HAVE A SPECIAL CASE-GO TO D BLOCK 
IFCTILT(K).EQ.O.) GO TO 26 

C CALCULATE THE AVERAGE HOURLY TOTAL DIFFUSE RADIATION ON A TILTED SURFACECNOT 
C REQUIRED FOR FINAL OUTPUT BUT MAY BE PRINTED BY SETTING M>1). 

RSI(I,K)=(RHI+ (RHI + DHI>XRHO)*U.-COSCTILT<K>)>/2. 

IFCTILTCKKEQ.O.) RSI(I,K)=RHI 
GO TO 444 

C d D 

C HORIZONTAL ARRAY CASE. 

26 CONTINUE 

C CALCULATE THE AVERAGE HOURLY DIRECT RADIATION. 

ASI(I,K) = DNI(I)X COSTH 

C THE FOLLOWING CHECK IS FOR ANOMALOUS NEG. VALUES NEAR SUNRISE 
C AND SUNSET. 

IFCASI (I , K) . LE. 0 . ) ASI ( I , K) =0 . 

C CALCULATE THE AVERAGE HOURLY DIFFUSE RADIATION ON A HORIZONTAL 
C SURFACECREFERENCE 1 STEP 4). 

RSICI,K)=(RHI+(RHI+DHI)xRH0)/2. 

IFCTILTCKI.EQ.O.) RSI(I,K)=RHI 
444 CONTINUE 

C D - D 

C 

C CALCULATE THE TOTAL AVERAGE HOURLY RADIATION ON A TILTED SURFACE. 

ZITU,K)=DNICI)#C0STH+.5X<(1.+C0S(TILTCK)))XRHI+ 

1 RHO*CDHI+RHI)x(l.-COS(TILT(K)))) 

IFCTILTCKI.EQ.O. . OR . ISWTCH . EQ . 1 ) GO TO 10 
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C E - - E 

C CALCULATE THE SHADOW CASE. 

DO 296 L=l,5 

296 5HDC I,K,L)=ZIT(1 , K)*5HDL0S( TILT!)!) , X! L ) , ALNGTH , COSZ) 

C E - - E 


C CALCULATE BETA, THE REFLECTOR TILT ANGLE 
BETA= ATAN!ALNGTH*SIN!TILT(K)>/X(K>) 

IF( ISWTCH . EQ . 2 ) GO TO 10 
IF(ARSIN(COSZ).LE.BETA) GO TO 392 
C CALCULATE REFLECTOR AUGMENTATION WITH SHADOW FACTOR. 
REF!I,K>=ZIT!I,K)XAUG!COSZ,TILT!IO,BETA,SREFL) 

GO TO 10 

C CALCULATE REFLECTOR CAUSED LOSS WITH SHADOW FACTOR. 
392 REF!I,K)=ZIT(I,tOxSHDLOS!TILT(K),XCl), ALNGTH»CQSZ) 
10 CONTINUE 
3 CONTINUE 


c c c 

c 

C F P 


IF(M.EQ.l) GO TO 92 
C PRINT THE APPROPRIATE HEADINGS. 

C PRINT OUTPUT 

IFCIPRINT.NE.1.AND.IPRINT.NE.3.AND.IPRINT.NE.5) GO TO 506 
IF(NONTH.GT.O) GO TO 90 
DO 90 1=9,10,3 

IFCIPRINT.EQ.O. OR. ISWTCH. GT.l) GO TO 5C6 
IF(N.EQ.I) WRITE( IUNIT ,91 ) 

IF(N.EQ.I) WRITE( IUNIT , 30 ) ( TOWN (JJ ). JJ=1 , 5 ), LAX .MINUTE 
IFC I . EQ . 9 . AND.N . NE.1.AND.N.NE.9. AND. N. NE. 7 .AND.N.NE.10) WRITE! IUNI 
IT, 101) 

90 CONTINUE 

92 IF(M.EQ.O) GO TO 93 
IF(N.EQ.l) GO TO 93 

WRITE! IUNIT ,91 ) 

IFIN.EQ.I) WRITE! IUNIT , 30 ) (TOWN! JJ ) ,JJ=1 , 5) , LAX, MINUTE 

91 FORMAT ! 1H1 ) 

101 FORMAT!/) 

IF! ISWTCH .GT.l) GO TO 506 

93 WRITE! IUNIT, 31 )N, PH, RHO.ZKTCN) ,ZKD, (ORIENT! J, NO), J=l, 3) 

31 FORMAT!/, 9X, ’ MONTH = ’, 19 , 3X, * HORIZ . HO= ', F5 . 2 , 10X, ' REFLECTIVITY ' ' 
1 , F5 . 2 > 

1 1 OX, ' KT : ’,F10.3,10X,’KD: ’ , F10 . 3//9X, 1 ORIENTATION: ', 

1 3A9.16X, 'HOURLY SOLAR RADI ATION *, 31X DAILY TOTAL'/ 


1 /21X, * 9 5 6 7 8 9 10 11’, 

« ’ 12 13 19 15 16 17 18 19 20’ ,9X, 'RADIATION' 

1 ,/10X, 'TILT *»5X, 83 ('-')) 

C p p 

C 

c G G 


C ZERO OUT DAILY TOTAL ARRAYS. 
506 DO 12 K=1 , NA 
SUMN = 0 . 

DO 297 L=l,5 
297 SUMS! K, L ) =0 . 

SUMO! K ) =0 . 

SUMR(K)=C. 


SUMREC K > — 0 . 
SUM(K)=0. 


C G-- G 

C 

C H - - H 


C SUM OVER THE HOURS IN THE DAY TO GET THE DAILY TOTALS. 

DO 14 1=1, 24 
SUMN=SUMN+DNI(I) 

SUMD(K ) =SUMD( K)+ASI ( I » K) 

SUMR(K)=SUMRCK)+RSI(I,K> 

DO 298 L = 1 , 5 

298 SUMS'K,L)=SUMS(K,L)+SHDCI,K,L) 

SUriRE(K)=SUMRE(K)+REF(I,K) 

14 SUMCK)=SUMCK)+ZIT(I,K> 

IFCM.EQ.l.AND.K.EQ.l) WRITEC IUNIT, 505 >TILD( K) ,( DNI ( I >, 1=5 , 21 > , 

1 SUMN 

505 FORMAT ( ' NORM. DIR . ' , F5 . 2 , 4X, F4 . 2, 16F5 . 2, F10 . 2) 

IFCIPRINT.EQ.O) GO TO 567 

C I - - - - 1 

C PRINT THE APPROPRIATE TABLES. 

GO TO (701,567, 703>567»705)»IPRINT 
701 WRITEC IUNIT , 32) TILDC K ) , (ZIT ( I , K ) , 1= 

1 5,21), SUMC K ) 

32 FORMAT (10X,F5.2,4X,F4.2,16F5.2,F10.2) 

35 F0RMATC/,2A4,/,3X, 'TOTAL ' , 17F5 . 0 , F10 . 0 , 4X, 4F6 . 0 ) 

GO TO 567 
703 DO 750 KL = 1 ,5 
WRITEC IUNIT , 751 ) 

751 FORMAT ( 1 OX, 'THE FOLLOWING TABLE CONTAINS SHADOW FACTOR') 

WRITEC IUNIT, 31) N , HO , RHO , ZKT ( N ) , ZKD, ( ORI ENTC J , NO) , J =1 , 3) 

WRI TEC IUNIT , 32 ) TILDC K ) , ( SHDC I , K , L ) , 1=5 , 21 ) , SUMSCK, KL ) 

750 CONTINUE 
GO TO 567 

705 WRITECIUNIT, 765) 

765 FORMATCIOX, 'THE FOLLOWING TABLE CONTAINS BOTH REFLECTOR AND SHADOW 
1 FACTORS.') 

WRITECIUNIT, 31) 

WRI TEC IUNIT, 32) TILDC K) , (REFC I , K) , I =5 , 21 ) , SUMREC K) 

567 CONTINUE 

IFCM. EQ.O) GO TO 12 

WRITEC IUNIT, 503 )TI LOCK), C ASI C I , K ) , I =5, 21 ) , SUMDC K) 

503 F0RMATC3X, 'DIRECT' ,F5.2,4X,17F5.2,F10.2) 

WRI TEC IUNIT, 504) TILDCK) , (RSI C I , K ) , 1=5 , 21 > . SUMRC K) 

504 FORMAT (3X, ' DI FFUSE ' , F5 . 2 , 4X , 17F5 . 2 , F10 . 2 ) 

c x X 

12 CONTINUE 
N = N + 1 

IFCN. GT. 12. OR. MONTH. GE. 1) GO TO 667 
GO TO 134 

C REDEFINE THE TILT VALUES IN DEGREES BEFORE RETURNING. 

667 DO 668 KP=1,NA 

668 TILT(KP)=TILD(KP) 

RETURN 

END 
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FUNCTION ZD CALLED FP.ON SUBROUTINE INS 
FUNCTION ZD(ZT) 

THIS FUNCTION INTERPOLATES AND RETURNS KD VALUES BASED ON THE 
KT VALUES PASSED AS INPUT. 

DIMENSION ZKT(6),ZKDC6> 

DATA ZKT/.3, .9, .5, .6, .7, .75/ 

DATA ZKD/. 179, .183, .188, .174, .149, .125/ 

IF(ZT-0 .3) 1,1,2 

1 ZD= . 125 
GO TO 10 

2 IFCZT-0.75) 3, 3, A 
A ZD= . 125 

GO TO 10 

3 DO 20 J=2 , 6 
T1=ZT-ZKT(J-1> 

T2=ZT-ZKT(J) 

TEST=T1«T2 
IF( TEST ) 5,6,20 

5 Y1=ZKDCJ-1> 

Y2=ZKD(J) 

ZD=Y1+CY2-Y1)K(ZT-ZKTCJ-1))/(ZKT(J)-ZKT(J-1)) 

GO TO 20 

6 IFCT1) 8,9,8 
9 ZD=ZKD(J-1) 

GO TO 20 
8 ZD=ZKD(J) 

20 CONTINUE 
10 RETURN 
END 
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FUNCTION SHDLOS CALLED FROM SUBROUTINE INS 
FUNCTION SHDLUSC ALPHA *XA>ALNGTH > SINE) 

THIS SUBPROGRAM COMPUTES AND RETURNS ONE MINUS THE SHADOW LOSS FRACTIONCTHE 
PERCENT RECEIVING INSOLATION) DURING HOUR SINE. 

XX=XA/ALNGTH 

E=AR5IN(SINE) 

$HDLOS=(COS( ALPHA )+XX)*SIN( E)/SIN(ALPHA+E) 

I F( SHDLOS . LT . 0 . ) SHDLOS = 0. 

IF(SHDLOS . GT . 1 , ) SHDLOS=1.0 

RETURN 

END 


A-21 










o oo o ooo 


FUNCTION AUG CALLED FROM SUBROUTINE INS 
FUNCTION AUG(G,T,B,SREFL> 

THIS SUBPROGRAM COMPUTES AND RETURNS ONE FLUS THE INCREASE IN HOURLY ENERGY 
DUE TO REFLECTOR AUGMENTATION. 

E=ARSIN(G) 

IFCE.LT.2.*B> GO TO 1 

AUG=SREFLNSIN(T)*SIN(F-B>*(1 ./SIN(B)-(COS(B)/SIN(B) 
1+C0SCT)/SIN(T))/(SJN(B)#C0S(E-2.*B)/SIN(E-2.#B) 
l + COSC B ) ) ) +1.0 

RETURN 

1 AUG=SREFL*SIN(E-B>*SIN(T)/SIN(B> +1.0 
RETURN 
END 
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ORIENTATION: SOUTHEAST HOURLY SOLAR RADIATION DAILY TOTAL 
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LOCATION: TUCSON 
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THIS PROGRAM GENERATES THE M FACTOR TABLE. 

DIMENSION DDEC<12),B<12).DLAT(51>.FACC12),Rtl2) 

DATA DDEC /-I 9 . 51 , -1 0 . 28 , . 2 , 

1 11. 56, 20.1 A. 23. 27, 20. 26, 12. 03, .37, -10. 97. -19. 58, -23. 27/ 
DATA B /.192, .199, .156, .18, .196, .205, .207, .201, 

1 .177, .16, .199,. 192/ 

DATA R /l. 03, 1.0207, 1.0057, .9875, .9727, .967, 

1 .9692, .9785, .9995, 1.0133, 1.0267, 1.0327/ 

PI =3 . 1916 
C SET THE OUTPUT UNIT 
IUNIT=6 

C SET THE SOLAR CONSTANTCAND THUS THE UNITS). 

SC=1.577 
WRITEC IUNIT , 60 ) 

60 FORMAT (/ , 20X, ' THE VALUE OF THE M FACTOR ’) 
WRITECIUNIT.70) 

7 0 FORMAT </,5X, 'EAT.', 27X , 'MONTH* ) 

WRITECIUNIT ,80) 

80 F0RMATC/5X, 'DEG. 12395678' 

1 ,' 9 10 11 12') 

WRITEC IUNIT ,S1 ) 

81 FORMAT (/ ) 

DO 10 1=1,51 
DLAT(I)=I-1. 

XL AT = INT CDIATCI ) ) 

RLAT = CXLAT+CDLATU)-XLAT>/.6)*PI/180. 

DO 20 J = 1 , 12 
XDEC r INT DDECC J ) ) 

RDEC=(XDEC+( DDECC J ) -XDF C )/ . 6 ) *PI/180 . 

COSWS=(-SINCRLAT)/COSC' L AT ) ) * ( SIN C RDEC )/COS ( RDEC ) ) 
WS=ARCOSCCOSWS) 

DWS=US*18C ./PI 
COSLD=COS(RLAT)*COSCRDEC) 

SINLD=SINCElAT)*SINCRDEC) 

S = 0 . 

DO 30 N=1 ,39 
W=iJS»«N/90 . 

COSZ = COSLD>iCOS(W)+SINLD 
PAR=-B( JJ/C0S2 
APA=ABS(PAR) 

IF(APA.GT.SO) GO TO 90 
AN5 = EXPCP.AR)*C0SZ 
GO TO 50 
90 ANS =0 . 

50 S=ANS+5 
30 CONTINUE 

ANS0=EXPC-BCJ)/CC0SLD+SINLD))*CC0SlD+$INLD>/2. 

AINT=WS/90 .*CANSO+5) 

THG=R( J)*<SC*(COSLD*SINCWS)+WS*SINlD) 

FAC( J ) =TH0/AINT 
20 CONTINUE 

WRITE! I UNIT, 90) DLAT(I) . ( FAC(MM) ,MM=1 , 12 ) 

90 FCRMAT(5X,F3.0,2X.12F5.2) 

DC 1 K-10,50,10 
IF(I.EQ.K) WRITEC IUNIT ,81 ) 

1 CONTINUE 
10 CONTINUE 
STOP 
END 


APPENDIX B 


SOLAR CLIMATE PARAMETER DATA FOR 235 LOCALITIES 


The following pages provide solar climate parameter data for 235 
localities in the United States (see Reference 10). Latitudes are in 
degree-minute form. 


K t * average daily total radiation for month M 

1 average daily total extraterrestrial radiation for month M 

M - 1,2, ...» 12 

LAT. 


LOCATION 

(DEG.MM) JAN 

FEB 

MAR 

APR 

MAY 

Jim 

JUL 

AUQ 

SEPT 

OCT 

NOV 

DEC 

A 04K 

jm 

*1.51 

0. 301 

0. 316 

0.142 

0.155 

0.119 

0.119 

0.122 

0.329 

-0.363_ 

0.397 

0.406 

0.112 

"nPUTri 

AK 

55.02 

0.29$ 

0.315 

0.368 

0.409 

0.428 

0.397 

0*418 

0.416 

0.417 

0.364 

0.368 

0*310 

0=Th£l 

AK 

60.47 

0. 303 

0.361 

0.444 

0.458 

0.432 

0.414 

0.303 

0.353 

0.420 

0.43* 

0.412 

0.111 

JUNr AJ 

Ak 

68.22 

0.267 

0.2M 

0.141 

0.1,7 

0 . 1,0 

0 . 1,5 

0.17. 

.9.166 

0-.137 

0-327 

0.340 

0.242 

kint. ' 

AK 

58.41 

0.349 

0.381 

0.451 

0.440 

0.437 

0.419 

0.407 

0.390 

0.4)9 

0.492 

0.494 

0.176 

K 30 1 A< 

AK 

57.45 

0.321 

0.341 

0.430 

0.444 

0.404 

0.416 

0.413 

0.4)0 

0.417 

0.411 

0.451 

0.344 


AK 

64.50 

0. 187 

0, 132 

3.*#, 

0.47. 

0.474 

0.477 

0.425 

0.400 

0.453 

0.472 

0*421 

0.001 

SUM*! T 

AK 

6 l.ifl 

0.270 

0.340 

0.4?5 

0*439'- 

0.490 

0* 445 

0.422 

0. 414 

0.457 

0.484 

0.511 

0.241 

VAKUTAT 

AK 

69.51 

0. 264 

0.281 

0.361 

0.394 

0.375 

0.366 

0.357 

0.356 

0.3*6 

0.376 

0.164 

0.247 


_iL 

55.34 

0.402 

0.675 

0.461 

0.506 

0.517 

0.521 

0.505 

0.526 

0.S20 

0*540 

0.492 

0.42* 

MOBIL: 

AL 

50.41 

0.434 

0.456 

0.465 

0.515 

0.523 

0.511 

0.400 

0.495 

0.303 

0.549 

0.504 

0. 444 

MCNriiJMcRV 

4L 

52.16 

0.412 

0.4 35 

0.470 

0.521 

0. 52 9 

0.537 

0.514 

0.530 

0*518 

0,549 

0.505 

0*441 

F OR T SN1 T H 

AR 

35.20 

3.447 

0.455 

0.476 

0.493 

0.533 

0.565 

0.576 

0.377 

0.54, 

0.556 

0.515 

H.41I 

TTT^PT^TT^ 

“TT 


TW 

0. 4*1 

d.4?3 

0*4^0 

0 V*37 

3.573 

TTTT7 - 

TTTTT 

0.550 

0.S‘" 

0,501 


PHOENIX 

AZ 

55.26 

0.579 

0.603 

0.644 

0.712 

0.746 

0.744 

0.694 

0.699 

0.720 

0.702 

0.657 

0.597 

PR:SC)TT 

A l 

14. 19 

0. 599 

0.600 

0.640 

0.692 

0.732 

0 . 74 , 

0.644 

0.641 

0.70, 

0.704 

0,675 

0.6241 

TUCiOi 

”a T 

3 i .5 7 

57 $99 

3.617 

TTTTT" 

TTT1T" 

0.745 

TTTTT" 

3.654 

0.4.2 

0.697 

0.695 

0.66) 

0.610 

MINSL3N 

AZ 

35.01 

0.56b 

0.601 

0.644 

0.695 

0.723 

0, 734 

0.655 

0.657 

0.701 

0.696 

0.671 

0*607 

VOL. 

_JuL 

32.40 

0. ,07 

0.621 

0.675 

0.727 

0.761 


3_,6A5 

Q.JQ0_ 

(3 « 72fc 

Q.712 

0-670 

0.424 

ARClTA 

CA 

40.59 

J. 389 

<3.412 

0. 4*»4 

0.499 

0. 516 

0.528 

0.506 

0*501 

0*528 

0.495 

0*440 

0.410 

BAKE R*>f I ElD 

CA 

15.25 

3.462 

0.503 

0.5 7*1 

0.639 

0.699 

0. 744 

0.749 

0.744 

0,726 

0.676 

0.572 

0*467 

China lake 

CA 

55.41 

0. 563 

0.565 

0.632 

0.68? 

0.710 

0.741 

0 *f 29 

0.805 

0*726 

0.6,7 

0.611 

0* 505 

6*4oetT ‘ 

TT 

16.57 

0. 567 

"TTTTT" 

0.640 

0.692 

6.622 

0. 749 

0. 726 

0. 7)1 

0.729 0.695 

0.647 

9*591 

EL TORO 

CA 

33.40 

0.540 

0.545 

0.573 

0.584 

0.577 

0.595 

0.660 

0.657 

0.622 

0.60? 

0.590 

0. 5*2 

MEJH3 

_Si_ 

56.4. 

0. 414 

0.4 76 

0.5 7, 

0.641 

0.692 

0. 71, 

Q. 749 

Q.T5Q 

.0,717 

0.601 

0.564 

0-417 

LCNT, BEACH 

CA 

33.44 

3.532 

3.537 

0.574 

0. 5a 7 

0.5 75 

0. 580 

0.642 

0.641 

0.610 

0.595 

0.579 

0.550 

LOS A <GR cS 

CA 

33.56 

0. 533 

0.538 

0.578 

0.591 

0.574 

0.575 

0.644 

0.635 

0.604 

0.592 

0.502 

0. 554 

MOONY SHASTA 

?» 

*UI9 

0« 410 

0*450 

0. 452 

0.554 

0.612 

0.655 

0.722 

0.701 

0.606 

0.614 

0.496 

0.44ft 

TITTiT! 

CA 

34.46 

3.561 

0.610 

0.658 

0.705 

0. 739 

0. 756 

0.709 

0.696 

0. 730 

0.703 

0.660 

0.614 

PA ALAND 

CA 

37,44 

0, 461 

0.489 

0.545 

0.593 

0,617 

0.634 

0.640 

0,636 

0.640 

0.592 

0.540 

0.409 

KOINT *UGU 


34 . 0 T 

0.53* 

0.542 

0.5,5 

0.592 

0.542 

0. 557 

0.591 

0*391 

0.579 

0.505 

0.516 

0.562 

RED BLUFF 

CA 

40.09 

0.406 

0.454 

0.524 

0.598 

0.664 

0.700 

0.74? 

0.729 

O.fl? 

0.636 

0.507 

0. 429 

S ACMM£NrO 

CA 

38.31 

0.400 

0.460 

0.55 1 

0.621 

0.679 

0. 723 

0.751 

0.740 

0.725 

0.654 

0.528 

0.420 

SAN DIEGO 

CA 

32.44 

0.542 

0.548 

0.575 

0.584 

0. 558 

0. 561 

0.610 

0,625 

0.609 

0.603 

0.594 

0.565 

TiT FUTiCTSM 


' 57. 37 

- 3 : 45 ?' 

0.484 

0.544 

TTTTT" 

tttTt 

0.641 

0 .*3?" 

TTTTT 

0..54 

3. 5.5 

TTTTT 

TTSJT 

SANTA MARIA 

CA 

34.54 

0. 506 

0. 515 

0.571 

0. 585 

0.596 

0.6)6 

O.f .53 

0.646 

0.628 

0.621 

0.501 

0.544 

SUNNYVALE 

_e» 

37.25 

0.475 

0.495 

0.553 

0.599 

0.615 

0.662 

0.681 

0*673 

0.659 

0,605 

0.547 

0*492 

C4ljAidJ 5I»o mtt 

CfT 

30.49 

0 . to5 

"TTTTT" 

TTTTT" 

0.6 00 

1.194 

0.61, 

0.618 

3.614 

0.671 

3.6,1 

0.645 

0.610 

DENVER 

CO 

39.45 

0.589 

0.566 

0.589 

0.597 

0. 596 

0. 633 

0.635 

0.643 

0,667 

0.667 

0.625 

0.601 

<4016 

-ta. 

39. 39 

0. 52 7 

0.542 

0.577 

0.6J1 

0.630 

0.675 

0*667 

0.655 

0.6,2 

0.66, 

0.612 

0.567 

GRAND JJNCMDN 

CT 

39.07 

TTTTT" 

TTTTT" 

0.492 

6.61, 

0. 664 

0.700 

3.6,9 

0.684 

o.toi 

0.679 

0.634 

0. 506 

PUEBL3 

CC 

38.17 

0.594 

0.570 

0.509 

0.606 

0.603 

0.656 

0.645 

0.656 

0.674 

0.6 73 

0.6)9 

0.605 

MAR T F 3RD 

CT 

41.56 

0.365 

0.381 

0,389 

0.417 

0.439 

0.453 

0.462 

0.454 

0.461 

0.463 

0.304 

0.162 

MASHI 4G7QN- STEF L I NX 

1C!) 

D. 3B^* 

TTTTT 

3.43, 

TTTTT - 

0.410 

0.512 

O.~50 8 

o. sot 

0* 512 

6. 555 

0.446 

T.nr 

# ILM ING T 0N 

OE 

39.40 

0. 399 

0.416 

0.44? 

0.462 

0.470 

0,507 

0,510’ 

0.500 

0.509 

0.503 

0*454 

0. 401 

APALACHICOLA 

n 

29.44 

0.436 

0.4fc| 

0.503 

0.561 

U. 304 

0. 54 7 

0. 508 

0,508 

0.528 

0.569 

0*5)5 

0*464 

D»yTo<U S?Ach 


29.ll 

0. 4Pi 

0.492 

TTTTT - 

6.561 

o.iSo 

0.501“ 

TTTTT" 

TTTTT - 

0 . soi 

TTTTT 

TTTJT" 

0 . 6,6 

JACKSONVILLE 

fl 

30.30 

0.469 

0.463 

0. 523 

0. 555 

0. 546 
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APPENDIX C 


AVERAGE DAILY SOLAR INSOLaIION DATA BY MONTH AND ANNUAL TOTALS 

FOR 235 LOCALITIES 


The following Cablet contain average daily insolation estimates 
and annual totals for south-facing surfaces at tilt angles from 0 to 
90 degrees for 235 localities in the United States. The intermediate 
angles are provided to facilitate interpolation. The average ground 
reflectivity assumed is 0.20. Each array tilt angle has 12 estimates 
associated with it — one for each month — plus the total annual 
amount. The average daily amount for the entire year is computed by 
dividing the annual total by 365. The values within the tables are in 
units of kilowatt hours per square meter. Localities are listed in 
alphabetical order by state and then by city within the state. 

The tables were printed using an IBM 3800 laser printer. 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLA'ION ESTIMATES (KWH/SO. 

SITE: ANNETTE AK 

LATITUDE: 55 DEGREES 2 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 



AVERAbE DAILY TOTAL TERRESTRIAL INSOLATIPH ESTIMATES (KMH/S9. M) 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/SQ. M) 

SITE: KING SALMON AK 

LATITUDE: 58 DEGREES A1 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/ 

SITE: KODIAK AK 

LATITUDE: 57 DEGREES 45 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: SUMMIT AK 

LATITUDE: 63 DEGREES 20 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: YAKUTAT AK 

LATITUDE: 59 DEGREES 31 MINUTES 






INSOLATION ESTIMATES (KWH/SQ. 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 

SITE: LITTLE ROCK AR 

LATITUDE: 34 DEGREES 44 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/SQ. M) 

SITE: PHOENIX AZ 

LATITUDE'- 33 DEGREES 26 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. 

SITE: PRESCOTT A Z 

LATITUDE: 3* DEGREES 39 MINUTES 



average daily ,om terrestrial isolation estimates irmrasr 

SITE: TUCSON *7 

LATI.UDt': 32 DEGREES 7 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: WINSLOW AZ 

LATITUDE: 35 DEGREES 1 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ES 1MATES (KWH/SQ. M) 




AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. 

SITE: BAKERSFIELD CA 

LATITUDE: 35 DEGREES 25 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: CHINA LAKE CA 

LATITUDE: 35 DEGREES 41 MINUTES 





average DAUT TOIAl TERRESTRIAl INSOLATION ESTIMATES (KUHXSO. M, 

SITE: DAGGETT CA 

LATITUDE: 34 DEGREES 52 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. 

SITE: LONG BEACH CA 

LATITUDE: 33 DEGREES 49 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. M) 

SITE: LOS ANGELES CA 

LATITUDE: 33 DEGREES 56 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (Kk’H/SQ. 

SITE: MOUNT SHASTA CA 

LATITUDE: 41 DEGREES 19 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KMH/SQ. M) 

SITE: NEEDLES CA 

LATITUDE: 3* DEGREES A6 MINUTES 



^ 4 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/SO. Ml 

SITE: OAKLAND CA 

LATITUDE- 1 37 DEGREES A« MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SO. M) 
SITE: POINT M'JGU CA 



«o 


total terrestrial insolation estimates (KWH/SO. M) 

, SITE: RED BLUFF n 

LATITUDE: AO DEGREES 9 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 
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SITE: SAN DIEGO CA 

LATITUDE: 32 DEGREES «<♦ MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: MIAMI FL 

LATITUDE: 25 DEGREES 48 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 

SITE: ORLANDO FL 

LATITUDE: 28 DEGREES 33 MINUTES 
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average oaky total terrestrial insolation estimates 
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LATITUDE: 43 DEGREES 34 MINUTES 
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daily total terrestrial insolation estimates (KWH/SQ. M) 

SITE-- INDIANAPOLIS IN 

LATITUDE: 39 DEGREES 49 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KKH/SQ. 

SITE: SOUTH BEND IN 

LATITUDE'- 41 DEGREES 42 MINUTES 



average daily total terrestrial insolation estimates OCWH/SQ. M) 

SITE: DODGE CITY KS 

LATITUDE: 37 DEGREES 46 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 
SITE: GOODLAND KS 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/S3. 

SITE: TOPEKA KS 

LATITUDE: 39 DECREES A MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/S9. 

SITE : WICHITA KS 

LATITUDE: 37 DEGREES 39 MINUTES 





average daily TOTAL terrestrial insolation estimates irwh/so. M) 

SITE: LEXINGTON KY 

LATITUDE: 33 DEGREES 2 MINUTES 





AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KU'H/SQ. 

SITE: LOUISVILLE KY 

LATITUDE: 3S DEGREES 11 MINUTES 



.30 2.28 2. AO 2.82 3.29 3.59 2.87 2.27 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: LAKE CHARLES LA 

LATITUDE: 30 DEGREES 7 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. M) 

SITE: NEW ORLEANS LA 

LATITUDE: 29 DEGREES 59 MINUTES 





AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. 

SITE: SHREVEPORT LA 

LATITUDE: 32 DEGREES 28 MINUTES 
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AVERAGE DA 1 1 if TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: PATUXENT RIVER MD 
LATITUDE: 38 DEGREES 17 MINUTES 
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VERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/S 

SITE: ALPENA MI 

LATITUDE: A 5 DEGREES A MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: DETROIT MI 

LATITUDE: 42 DEGREES 25 MINUTES 



.SI 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: GRAND RAPIDS MI 

LATITUDE: 42 DEGREES 53 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. 

SITE: HOUGHTON MI 

LATITUDE-- A 7 DEGREES I 1 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. MJ 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (JCWH/’SQ. 

SITE: TRAVERSE CITY MT 
LATITUDE: 44 DEGREES 44 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SO. PI) 

SITE: DULUTH MN 

LATITUDE: 46 DEGREES 50 MINUTES 



A\ AGE DAILY TCTAL TERRESTRIAL INSOLATION ESTIMATES CICMH/SO. Ml 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIHATES (KWH/SO. Ml 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/SQ. M) 

SITE: ROCHESTER MN 

LATITUDE: A3 DEGREES 55 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: COLUMBIA MO 

LATITUDE: 38 DEGREES 49 M T NUTES 




AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKUH/SQ. 

SITE: KANSAS CITY MO 

LATITUDE'- 39 DEGREES IS MINUTES 



AVERAGE PAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 
SITE: SPRINGFIELD MO 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTI,-,a.cS iKWH/SQ. M) 

SITE: ST. LOUIS MO 

LATITUDE: 38 DEGREES 45 MIKoTES 





AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKUH/SQ. M) 

SITE: JACKSON MS 

LATITUDE: 32 DEGREES 19 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 

SITE: RALEIGH-DURHAM NC 
LATITUDE: 35 DEGREES 52 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: BISMARCK HD 

LA.irUDE: A6 DEGREES A6 MINUTES 



•93 2.93 3.37 3.87 A. 15 A. 12 3.13 2.59 1203 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M» 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/SQ. 

SITE: RENO NV 

LATITUDE: 39 DEGREES 30 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ . M) 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. 

SITE: WINNEMUCCA NV 

LATITUDE: AO DEGREES 5A MINUTES 





AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/SQ. 

SITE : YUCCA FLATS NV 

LATITUDE: 36 DEGREES 57 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: BUFFALO NY 

LATITUDE: 42 DEGREES 56 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: TOLEDO OH 

LATITUDE: 41 DEGREES 36 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 



UJ 

o 















<>“ 

HO 

m 

m 


s© 

s© 

« 

S3 

in 

S3* 

N 

© 

n. 

fO 

u < 

• 

» 

* 

• 

* 

. 

• 

. 

. 

* 

* 

. 

• 

« 

UJO 

> 

< 

K) 

HO 

HO 

fO 

tO 

rO 

HO 

ho 

HO 

HO 

HO 

rO 

<SJ 

CSJ 


-»h 


s3* 

•o 

CM 


« 

O' 

«o 

O 

*6 

H 

CM 

O' 


<z 

• 

* 


* 

• 

* 

• 


* 

• 

* 

♦ 

• 

• 


*A 

o 

CM 

•o 

s© 


H 

•0 

o 

«3* 

•o 

O' 

O 

<0 

zo 

o 

r-* 

O' 

O 

rH 

H 

H 

O' 

M 

in 

00 

« 

P- 

»3* 

ZL 

CM 

<M 

CM 

K> 

HO 

HO 

HO 

CM 

CM 

CM 

H 

o 

O' 

•0 

< 

H 

H 

H 

H 

rH 

rH 

H 

pH 

»H 

H 

H 

rH 






1 

1 

o 

CM 

HO 

CM 

rt 

O' 

sO 

CM 

rv 

rH 

n 

Mf 

O' 

O' 



1 

o 

CM 

HO 

sf 

n 

in 


r-* 

h' 

*0 

00 

00 


M> 



Q 1 

* 

• 

• 

» 

• 

• 


« 

* 




• 




1 

1 

H 

rH 

rH 

H 


»H 

rH 

rH 


rH 

rH 

rH 

rH 

rH 



1 

1 

s3* 

•o 

O' 

00 

h- 

in 

H 

s© 

0 

ro 

-3* 

O 

rH 

00 



1 

o* 

SO 


•0 

O' 

o 

rH 

rH 

CM 

CM 

CM 

CM 

rH 

Os 



x i 

. 

* 

• 

• 

• 

• 

• 

• 

• 

• 

• 

» 

* 

* 



i 

1 

rH 

rH 

rH 

rH 

H 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 



1 

1 

O' 

m 

O 

'd* 

SO 

sO 

»o 

O' 

CM 

ro 

O' 

s© 

in 

s© 



1 

S» 

© 

CM 

H0 

S3* 

n 

SO 

s© 

r*. 

r>* 

s© 

n 

ro 

O 



O 1 

* 

• 

• 

♦ 

• 

« 

• 

• 

* 

* 

• 

. 

• 

* 



1 

1 

CM 

HO 

HO 

K> 

HO 

HO 

HO 

HO 

H0 

ro 

to 

H0 

HO 

ro 



1 

1 

sD 

SD 

h» 


s3* 

00 

o 

O' 

«© 

rH 

CM 

n 

O' 

so 



1 

h* 

O 

r-t 

CM 

HO 

H0 


HO 

HO 

HO 

rH 

00 

o* 

o 



to 1 

. 

. 

• 

* 

• 

* 

* 

» 

• 

* 

* 

» 

« 

• 



1 

HO 

'O 

s3* 

«o 

s3* 

'tf* 

«3- 


«sT 

s3* 

•O 

rO 

HO 

ro 

tn 


1 















XUJ 

X 

1 















Oh 

h- 

1 















3 

2 

1 

m 

CM 

s© 

00 

n- 

sr 

00 

O' 

00 


o 

00 

oo 

CM 

X 

O 

1 

rs. 

O' 

O' 

O' 

O' 

O' 

«0 

h. 

s© 

in 

CM 

r- 

CM 

h- 


E 

< 1 

» 

* 

* 

» 

• 

♦ 


• 

» 

» 

» 

• 


* 

E 


t 

S3* 

s3* 

vT 

'J* 

<r 

<r 

*3* 

s3* 

*3“ 

s3* 

S3* 

ro 

HO 

CM 


>- 

1 















sO 

CO 

1 















H 


I 
















t/3 

1 

Hv 

HO 

r*H 

h* 

O' 

O' 

sD 

O 

HO 

ro 

Hv 

S3* 

M3 

h. 

xm 

H* 

1 

>3* 

in 

m 

s3* 

H0 

CM 

rH 

0 

00 

M3 

rH 

s© 

o 

»3* 

3UJ 

2 

■o 1 

» 

. 

• 

» 

• 

• 

• 

* 

» 


• 

• 

• 

* 

oiu 

3 

1 

in 

m 

in 

in 

n 

in 

in 

in 

s3* 

s3* 

S3* 

ro 

ro 

CM 

H- 

O 

» 















to o 

E 

f 















O uj 

<E 

I 















2 Q 


1 

m 

m 

CM 

in 

in 

ro 

00 

o 

»H 

O 

CM 


o 

rH 

3 

>- 

1 

in 

in 

in 

S3* 

h0 

CM 

o 

O' 

Hh 

in 

0 

sr 

O' 

ro 

O f* 

-J 

“0 I 

. 

• 

* 

• 

• 

• 

. 

* 

. 

* 

• 

• 

. 

* 

>- nJ- 

M 

1 

m 

in 

in 

in 

n 

n 

in 

»3- 

v3* 

s3“ 

s3* 

H0 

CM 

CM 


< 

1 















•» 

Q 

1 















UJ LU 


1 















h-Q 

LU 

1 

o 

in 

s3* 

o 

H0 

sr 

CM 

00 

CM 

n3* 

HO 

in 

HO 

© 

>H Z> 

O 

1 

o 

o 

O 

o 

O' 

«o 

n- 

in 

>3“ 

CM 

00 

HO 

00 

ro 

i/) h* 

< 

E I 

• 

* 

• 

• 

» 

* 

. 

♦ 

* 

« 

« 

» 

* 

• 

HH 


1 

m 

in 

m 

in 

•sT 

-j* 

s3* 

s3“ 

sJ- 

>3- 

rO 

ro 

CM 

CM 

h* 

UJ 

1 















< 

> 

1 
















< 

1 

















1 

ro 

sO 

O' 

o 

O' 

n 

o 

CM 

rO 

CM 

HO 

00 

oo 

ro 



1 

o 

rH 

H 

CM 

K 

rH 

rH 

o 

O' 

00 

in 

rH 

hv 

ro 



< 1 

• 

• 

» 

♦ 

• 

• 

• 

« 

• 

* 

» 

• 

• 

* 



1 

1 

s3* 

s3- 

s3" 

s3* 

stf* 

s3* 

'J* 

s3* 

HO 

HO 

ro 

rO 

CM 

CM 



1 

1 

O 

00 

stf - 

O' 

rO 

ST 

in 

HO 

O 

in 

rH 

CM 

fs. 

r*v 



1 

«0 

O' 

O 

0 

H 

H 

rH 

rH 

rH 

o 

O' 


s3- 

rH 



E l 


* 

. 

» 

• 

. 

♦ 

» 

• 

t 

. 

. 

» 

» 



1 

1 

CM 

CM 

HO 

HO 

HO 

HO 

ro 

ro 

HO 

HO 

CM 

CM 

CM 

CM 



1 

1 

in 

CM 

O' 

in 

0 

s3* 

<o 

00 

00 

Pv 

CM 

CM 

00 

o 



1 

oo 

o 

o 

rH 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

O' 

00 



U, 1 


• 

* 

« 

• 

» 

• 

* 

• 

» 

• 

. 

* 




1 

1 

rH 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

rH 



1 

1 

CM 


h* 

h* 

SO 

HO 

o 

s© 

o 

s3* 

Mi 

S3* 

H* 

Mf* 



1 

CM 

s3* 

m 

%© 

H'. 

•0 

O' 

O' 

o 

o 

o 

O 

O' 

00 



-o 1 

• 

• 

. 

• 

• 

* 

* 

* 

• 

• 

• 

* 

. 




1 

1 

»h 

rH 

H 

H 

H 

rH 

rH 

rH 

CM 

CM 

CM 

CM 

rH 

rH 


> 

< ♦— oooooooooooooo 

(X-i ... * * 

QC »-h oomoinoinomooooo 

< h- HHNNWKltf^invDNIO^ 




AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KI4H/SQ. M) 

SITE: OKLAHOMA CITY OK 
LATITUDE: 35 DEGREES 24 MINUTES 





AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. M) 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKUH/S9. M) 

SITE: ASTORIA OR 

LATITUDE: 46 DEGREES 9 MINUTES 



.56 2.40 2.75 3.01 3.44 2.89 2.03 1.37 884.3 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. 

SITE: BURNS OR 

LATITUDE: A3 DEGREES 35 MINUTES 
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average daily total terrestrial insolation estimates ckmh/sq. mj 

SITE: PORTLAND OR 

LATITUDE: 45 DEGREES 36 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. M) 

SITE : SALEM OR 

LATITUDE- 44 DEGREES 55 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. K) 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/SQ 

SITE: ERIE p. 

LATITUDE: A2 DEGREES 5 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL ISOLATION ESTIMATES LRUHySR 

, . ti SITE: HARRISBURG PA 

LATITUDE: 40 DEGREES 13 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: PHILADELPHIA PA 

LATITUDE: 39 DEGREES 53 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: PITTSBURGH PA 

LATITUDE: 40 DEGREES 30 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. 

SITE: WILKES-BARRE-SCRANPA 
LATITUDE: 41 DEGREES 20 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKUH/SO. M) 

SITE! RAPID CITY SD 

LATITUDE: 44 DEGREES 3 MINUTES 



82 2.82 2.19 3.78 4.41 4.57 3.80 3.04 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CRNHySR. M> 

SITE: SIOUX FALLS SD 

LATITUDE: A3 DEGREES 3A MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ . M) 

SITE: MEMPHIS TN 

LATITUDE: 35 DEGREES 3 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 

SITE: NASHVILLE TN 

LATITUDE: 36 DEGREES 7 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M> 

SITE-’ ABILENE TX 

LATITUDE: 32 DEGREES 26 MINUTES 





AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKWH/SQ. 

SITE: AMARILLO TX 

LATITUDE: 35 DEGREES l<t MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 

SITE: BROWNSVILLE TX 

LATITUDE: 25 DEGREES 56 MINUTES 
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AVERAGE DAILY tqt.l TERRESTRIAL ISOLATION ESTIMATES (KUH/S9. M> 

SITE: SAM ANGELO TX 

LATITUDE: 31 DEGREES 22 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. 

SITE: SAN ANTONIO TX 

LATITUDE: 29 DEGREES 32 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KMH/SQ. 

SITE: SHERMAN TX 

LATITUDE: 33 DEGREES A3 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. H> 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KMH/SQ. M) 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. 

SITE: BRYCE CANYON UT 

LATITUDE: 37 DEGREES 42 MINUTES 
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AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. 

SITE: SALT LAKE CITY UT 
LATITUDE: AO DEGREES A6 MINUTES 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 




<>- 

CM <S- 

m 

>c 



* 

M 

m 

*■ 

90 

** 

O' 

Qt < 

tuo 

• * 

-4* 

st 

• 

*■ 

• 

* 

• 

ST 

• 


* « 

tO 

M 


< 


«3h 

© 

m 

O 

CM 

e 

tO 

CM 


>n 

CM 

#n 

pH 

*• 

Pm 

«z 

• 

• 

• 

• 

* 

* 

* 

♦ 

• 

• 

* 

• 

• 

• 

S3 

<6 

•n 

** 

in 


© 

CM 


CM 


<£ 

pH 

pH 

© 

zo 

CM 

H 


M) 

K 

is. 

N 

m 

tO 

© 

e 

•0 

to 

m 

Z£ 

in 

« 



M3 


*6 


M> 

in 

IT) 

to 

CM 

e 

< «t 

H 

H 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 




1 

pm 

ps 

in 

CM 

ps 

o 

CM 

pH 

O' 


O' 

M3 

to 

to 



1 

© 

to 

in 

n. 

oo 

o 

pH 

cm 

CM 

to 

to 

to 

CM 

o 



Q 1 

• 

• 

* 

• 


• 

• 


• 

• 

• 

* 


• 



1 

» 

pH 

CM 

CM 

CM 

CM 

to 

to 

to 

to 

to 

to 

to 

to 

to 



1 

1 

m 

K) 

m 

m 

CM 

Ps 

© 

pH 

O' 


rs. 

© 

CM 




1 

m 

© 

CM 

o- 

«C 

Ps 

O' 

O 

o 

PH 

pH 

pH 

© 

so 



2 1 

* 

• 

* 

• 

♦ 

• 

* 

♦ 

• 

• 

• 

• 

* 

• 



1 

1 

CM 

to 

to 

to 

to 

to 

to 


O’ 

'tf- 

*0* 

'd* 

to 

to 



1 

1 

M 

o* 

CM 

•0 

pH 

CM 

o 

sO 

(0 

•0 

O 

pH 

pH 

CM 



i 

O' 

•0 

o 

pH 

to 

O’ 

in 

m 

m 

in 

m 

to 

O 




O 1 

« 


• 

• 

* 

• 

* 

• 

• 

• 

• 

• 

• 

• 



1 

1 

to 

to 

** 

O' 

O’ 

O' 

o* 

O' 

o- 


»o- 


Atf* 

to 



1 

1 

© 

o 

CM 

pH 

Ps 

O 

© 

ps 

CM 

o* 

O' 

o- 

O 

•a 



1 

to 

r% 

•0 

© 

© 

© 

o 

O' 

© 

to 

m 

CM 

00 

CM 



to 1 

• 

• 


* 

• 

* 

• 

• 

• 

* 

• 

• 

» 

• 



1 

o* 

o- 

O’ 

St 

Iff* 

m 

in 

o- 

O" 

«o 

H9- 

o- 

to 

tO 

40 


1 















C Ui 

X 

I 















> 1“ 

K 

1 















=D 

2 

1 

© 

m 

to 

oo 

in 

c r 

o 

O' 

m 

•0 


in 

r% 

Atf* 

2 

O 

1 

to 

o- 

St 

o 

o* 

to 

to 

PH 

o 

00 

•O 

© 

to 


*-H 

£ 

< 1 

* 

• 

* 

• 

• 

• 

• 

« 

♦ 

• 

• 

« 

• 

• 

£ 


1 

m 

in 

in 

m 

m 

in 

m 

m 

m 

'tf* 

•o 

to 

to 

CM 


> 

I 















*0- 

A 

1 















in 


1 
















m 

1 

o> 

M3 

CM 

in 

in 

pH 

in 

rs 

ps 


pH 

pH 

© 

to 

tfi 

*- 

1 

•0 

00 

90 

rs 

M3 

m 

rO 

pH 

O' 

n* 

CM 

M3 

© 

to 

UJ 

2 

“1 I 



* 

» 

* 

* 

• 

* 

* 

• 

« 

• 

» 

• 

UJ 

2 

1 

m 

m 

m 

in 

m 

m 

m 

in 

o- 

'O* 

& 

to 

CM 

CM 

^ Ckf 

O 

» 















~» O 

£ 

1 















O UJ 

< 

1 















u. a 


1 

pH 

pm 

pH 

© 

rs 

© 

© 

®© 

O’ 

«0 

&0 

CM 

'tf* 

rv 

c* 

>- 

1 

to 

CM 

CM 


© 

•0 

M3 

to 

pH 

00 

CM 

<M3 

© 

CM 

O M3 

_l 

—) 1 

• 


* 

• 

* 

♦ 

• 

» 

• 

• 

• 

* 

• 

» 

2 to 

PH 

1 


-c 

M3 

HO 

m 

m 

in 

m 

in 

O' 

<r 

to 

CM 

CM 


< 

1 















•• •• 

a 

1 















iuuj 


1 















t- Q 

UJ 

1 

m 

rs. 

to 

M> 

M3 

CM 

sO 

ps 

Ps. 

to 

o* 

•0 

Mf 

n* 

»-h Z5 

o 

1 

© 

© 

© 

00 

rs 

MO 

o- 

CM 

o 

00 

CM 

M3 

O 

to 

in t- 

< 

£ 1 

• 

* 

* 

• 

• 

« 

• 

• 

• 

♦ 

* 

* 

• 

• 

M 

K 

» 

m 

in 

m 

m 

m 

m 

m 

in 

in 

o* 


to 

to 

CM 

t— 

UJ 

1 















< 

> 

1 















— J 

< 

1 

















1 

© 

o 

ps 

Ps 

O’ 

00 

© 

rs 

to 

nO 

m 

m 

rs. 

•d* 



1 

CM 

o* 

o* 

o* 

O’ 

to 

CM 

pH 

© 

00 

V3- 

© 

to 

rs 



< 1 

• 

* 


• 

• 

* 

• 

* 

• 

« 

• 

* 

• 

• 



1 

in 

m 

ut 

in 

in 

m 

m 

in 

m 



to 

to 

CM 


1 


3 


I 

o 


» 

to 

© 

© 

90 

to 

m3 

sD 

rO 

00 

o 

00 

o 

rs 

o 

1 

© 

to 

© 

'tf- 

m 

in 

m 

m 

O' 

o 

pH 

oo 

O' 

o 

£ 1 

. 

* 


« 

• 

» 

• 

* 

• 


• 


♦ 

* 

1 

1 

o 





■d* 

sS* 

%** 

O' 

O' 

O' 

to 

to 

to 

1 

1 

O’ 

M> 

O 

CM 

pH 

© 

in 

00 

© 

to 

© 

to 

00 

M> 

t 

© 

CM 

vr 

in 

M3 



rs 

fs. 

r s 

r* 

m 

CM 

© 

U. 1 

» 

• 

« 

• 

• 

♦ 

• 

• 

• 

• 

• 

• 

• 

. 

1 

1 

CM 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

CM 

1 

1 

to 

© 

in 

© 

pH 

CM 

pH 

© 

O' 

Pm 

00 

pH 

M3 

o- 

1 

pH 


M3 

© 

© 

© 

pH 

pH 

CM 

CM 

CM 

CM 

© 

00 

-3 1 

« 

• 

• 

• 

« 

» 

• 

• 

• 

• 

• 

* 

• 

• 

1 

CM 

CM 

CM 

CM 

CM 

to 

to 

to 

to 

to 

to 

to 

to 

CM 


© 

o 

O 

O 

O 

o 

© 

o 

o 

o 

© 

O 

o 

O 

© 

© 

m 

O 

in 

© 

m 

© 

in 

© 

© 

© 

© 

© 


pH 

pH 

CM 

CM 

to 

to 

O' 

o 

in 

M> 

rs 

•0 

© 




217 


K> 


CM 

* 


O' 

m 

S 

X 


w 

%n 


X 

H 

t- 

tO 


o 

►H 

t- 

< 


ifi 
C LU 
> t- 


O 

in 


X 

»H 

X 


< 

t-H 

0* 

H- 

tO 


to 

«/> 


QUJ 
X O' 

o o 


ft: 

ft’ 


X Q 
U 

*-• r* 
ft: tn 


_i uj lu 

< ►- Q 

H- H3 
O tO H- 


*< 

Q 


UJ 

O 

< 

ft: 


> 

< 


UJ 

o 

2$ 

LUO 

> 

< 


<M m 


t*> 

» 

** 


fO 


CM «H 


ft* 

*n 


ft 

m 


CM 

* 

tO 


Jh 

<X 

CM 

* 

40 

r* 

*9 

m 

rs 

• 

CM 

rs. 

ft 

ft 

m 

ft 

33 


•0 

4ft 

•ft 


o» 

tn 

•0 

ft 

m 

O' 

»o 

X o 


H 

49 

•ft 


r-» 

fN 

in 

tO 

o 

ph 

o 

z x 


m 

m 

in 

m 

m 

m 

m 

m 

in 

49 

to 

■< HK 

rt 


H 

H 

ph 


H 

*h 

H 

H 

fH 

H 


ft 

49* 

ft 


1 

ft 

ft 

to 

*0 

PH 

pH 

O' 

O 

r- 

Ps 

O' 



1 

O 1 

CM 

ft 

*0 

O' 

H 

CM 

CM 

to 

to 

to 

Cii 

P-4 

•0 

1 

1 

to 

tn 

to 

to 

ft 

ft 

49* 

ft 

49- 

ft 

40 

49- 

to 


X 

o 

X 


>- 

« 

to 


X 

3 

o 

X 

< 


< 

o 


LU 

o 

< 

ft 

UJ 

> 

< 


fO 

CM 


m 

pH 

ft 

O' 

pH 

H 

o 


pH 

ft 

CM 

H 

to 

49 

in 

n. 

to 

O' 

O' 

O 

© 

o 

* 

* 

» 

• 

• 

• 

• 

• 

, 



CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

to 

to 

to 

to 

CM 

O' 

ft 

to 

O' 

to 

m 

© 

CM 

ft 

N 

O' 

O 

CM 

to 

ft 

in 

ft 


n. 

ft 

CM 

CM 

to 

fO 

to 

to 

to 

to 

to 

to 

to 


O' 

CM 


in 

to 


49 

CM 


%9 

m 


in 

ft 


«0 tO 


fO 

•0 


•0 

49 


ft 

f4- 


•0 

ft 


in 

•9* 


CM 

pH 

49 


O' 

ft 


to to 


i 

in 

O' 

CM 

to 

© 

•9* 

ft 

in 

CM 


IS. 

o 



ft | 


pH 

CM 

CM 

CM 

pH 

o 

O' 

•0 

ft 

CM 

•0 

CM 

ft 

1 

1 

I 

in 

m 

in 

in 

m 

in 

in 

•9* 

49* 

49* 

49“ 

to 

to 

CM 

1 

1 

O' 

pH 

to 

pH 

CM 

o 

m 

rs 

to 

fN 

to 


CM 


1 

"» 1 

in 

ft 

m 

in 

49* 

to 

pH 

O' 

n. 

m 

o 

m 

O' 

to 

1 

1 

| 

in 

in 

in 

m 

m 

in 

in 

49* 

49* 

ft 

49* 

to 

CM 

CM 

1 

1 

o 

to 

CM 

to 

pH 

in 

f«s 

fs. 

ft 

CM 


ft 

in 

to 

1 

I 

O' 

to 

to 

t>. 

ft 

49* 

CM 

O 

to 

ft 

O 

49* 

to 

CM 

» 

1 

I 

in 

m 

m 

tn 

m 

m 

m 

in 

•9* 

49* 

ft 

to 

CM 

CM 

1 

1 

in 

n. 

49* 

•0 

O' 

ft 

CM 

m 

ft 

in 

ft 

O' 



1 

X 1 

m 

in 

m 

49* 

to 

CM 

pH 

O' 


in 

o 

49* 

O' 

to 

1 

l 

1 

m 

m 

m 

m 

m 

in 

in 

ft 

ft 

49* 

49* 

to 

CM 

CM 

1 

1 

ft 

<0 

pH 

pH 

O' 

to 

in 

*9* 

pH 

ft 

to 

CM 

© 

CM 

ft 1 

O' 

O 

PH 

pH 

o 

© 

O' 

to 

P4. 

in 

pH 

f". 

pH 

ft 

1 

1 

1 

ft 

m 

m 

m 

in 

m 

49* 

49* 

Vf 

49* 

ft 

to 

to 

CM 

1 

t 

pH 

ft 

m 

to 

rv 

o 

© 

to 

to 

ft 

m 

ft 

O' 


X i 

•0 

o 


CM 

CM 

to 

to 

CM 

CM 

pH 

O' 

ft 

CM 

to 

» 

i 

1 

to 

49* 

ft 


•9* 

49* 

4tf* 

49* 

49* 

49* 

to 

to 

to 

CM 

1 

1 

rv 

fN 

o 

© 

© 

n. 

CM 

m 

ft 

ft 

to 

CM 


O' 

U_ | 


© 

CM 

to 

49* 

49* 

in 

tn 

in 

in 

49* 

to 

o 

r-. 

i 

i 

i 

CM 

tO 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

CM 

i 

i 

O' 

CM 

ft 

04 

pH 

pH 

O' 

ft 

pH 

49* 

m 

to 

tn 

49* 

1 

O' 

to 

ft 

in 

n* 

to 

to 

O' 

© 

© 

© 

O' 

to 

ft 

1 

1 

pH 

CM 

CM 

CM 

pj 

CM 

CM 

CM 

tn 

to 

to 

PJ 

CM 

CM 

>- 

ft t- 
QC -J 

© 

O 

o 

o 

© 

o 

o 

© 

© 

o 

o 

© 

o 

o 

0C PH 

O 

o 

m 

o 

in 

o 

m 

O 

in 

o 

© 

© 





pH 

pH 

CM 

CM 

to 

to 

49* 

•o* 

m 

ft 


to 

O' 




.20 3. *7 3.25 2.74 1005.2 2.8 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/S9. M) 



-Jh 


m 

N 

(M 

tO 

tO 

tO 

rvi 

© 


e 

rg 

m 

rs 

*x 


• 

* 

4 

* 

• 

• 


* 

« 

■ « 

4 



ss 

ft* 


p* 

to 

© 

to 

p* 

IO 


e 

to 

in 

to 

o 

zo 

*£} 


p% 


RH 

M 

© 

cp 

r» 

♦ 

M 

to 


to 

Z£ 

3 

in 

in 

in 


3 

s© 

m 

m 

in 

♦ 

to 


o 

<< 

•H 

H 

p*4 

^4 

^4 

H 

«-4 


r-4 

H 

H 







» 

ph 

in 

CM 

to 

CM 

s3* 

m 

3 

H 

fNi 

*4 

CO 

fs. 

CO 



1 

•0 

CM 

3 

m 

p* 

to 

O 1 

© 


fH 

CM 

H 

© 

CO 



» 

1 

*-4 

CM 

ISJ 

CM 

N 

CM 

CM 

to 

to 

to 

to 

to 

to 

CM 



i 

i 

*"4 

M> 

3 

3 

© 

3 

sO 

3 

to 

co 

f* 

s3 

CO 

CM 



i 

3 

•0 

© 

CM 

s3- 

in 

sO 

p* 

co 

C0 

OS 

C0 

sO 

3 



Z 1 

t 

* 

• 

♦ 

• 

• 

* 

• 

• 

* 

• 


i 

• 



1 

1 

CM 

CM 

to 

to 

to 

»o 

to 

to 

to 

to 

to 

to 

to 

to 



» 

1 


3 

to 

O' 

to 

3 

CM 

to 


p-4 

to 

a3- 

m 




1 


CO 

o 

r-4 

to 

A* 

m 

m 

S0 

sO 

m 

to 

© 

© 



O I 

♦ 


• 

* 

* 

♦ 

• 

• 

• 

• 


• 

* 




1 

1 

tO 

to 

3 

3 

3 

3 

3 

s3- 

AO 

a3 

3 

A3- 

3 

to 



1 

1 

•0 

to 

O' 

to 

3 

r% 

p* 

in 

Os 

CM 

CO 

•3* 

CM 

^4 



1 

CM 

in 

© 

p* 

•0 

to 

c§ 

•0 

ph 



H 

r*- 

CM 



in t 

* 


• 

* 


• 



• 

• 

• 

* 

* 




* 


3 

3 

3 

3 

3 

s3* 

A3- 

s3- 

a3* 

* 

a3 

to 

to 

<n 


i 















< Ui 

X 

i 















> i- 

t— 

i 















=> 

Z 

i 

pH 

3 

CP 

O' 

p. 

H 

CM 

H 

CO 

CM 

CM 

S3* 

O' 

CO 

Z 

O 

I 

rH 

CM 

CM 

CM 

CM 

CM 

F-l 

O 

to 

n. 

to 

CO 

CM 

S0 

t-4 

e 

< i 

« 

» 

• 

* 

• 

« 

• 

* 


* 

• 




r 


I 

in 

in 

in 

m 

m 

m 

in 

in 

•3- 

•3* 

A3* 

to 

to 

CM 


V 

i 















CP 

CO 

i 















f*H 


i 
















in 

i 

3 

to 


to 

•0 

sO 

© 

to 

to 

CM 

p-4 

to 

3 


in 

i- 

i 

3 

3 

3 

in 

sJ- 

to 

CM 

o 

to 

sO 

H 

m 

O' 

to 

UJ 

z 

-» i 

» 

• 

• 

* 

» 

• 

• 

• 


* 

* 

• 

• 

* 

LU 

r> 

» 

m 

in 

m 

m 

in 

in 

m 

in 

s3- 

a3 

a3* 

to 

CM 

fg 

LU C* 

c 

i 















o 

£ 

i 















o UJ 

< 

» 















z o 


t 


t*4 

in 

M3 

to 


© 

o 

C0 

to 

to 

tv 

m 

to 

c 

V 

i 


O' 

to 

rv 

3 

s3- 

to 

H 

to 

sO 

© 

A3- 

to 

CM 

o p* 


"» i 

* 

* 

» 

* 

♦ 

* 

♦ 

♦ 


* 

c 

» 

• 

* 

sx »o 

►H 

i 

in 

in 

m 

in 

in 

in 

in 

m 

Atf- 

s3- 


to 

CM 

CM 


< 

i 















H •• 

Q 

i 















UJ UJ 


i 















b- Q 

UJ 

» 

3 

•0 

m 

O' 

o 

p* 

CM 

m 

AO 

m 

AO 

Cts 

H 

O' 

►H 3 

o 

i 

m 

m 

m 

' T 

to 

CM 

p-4 

Os 

p* 

in 

o 

S3- 

3* 

CM 

ini- 


E 1 

• 

• 

• 

■ 

♦ 

* 

• 

• 

• 

• 

c 

• 

• 

« 

*-i 

te 

1 

m 

m 

m 

m 

m 

m 

*.n 

a3- 

s3- 

s3 

S3* 

to 

CM 

CM 

H- 

UJ 

1 















< 

> 

1 















-i 

< 

1 

















1 


»o 

3 

3 

to 

p. 

Os 

<0 

m 

© 

H 

in 

^4 

to 



1 

cp 

H 

H 

t—i 


o 

O' 

to 

P. 

AO 

CM 


CM 

A0 



< 1 

« 

• 

• 

• 

* 

• 

* 


• 

• 


• 

* 




1 

1 

3 

in 

m 

m 

m 

m 

s3- 

3 

3 


a3- 

to 

to 

CM 



f 

1 

O' 

in 

m 

CM 

n. 

o 

o 

to 

to 

in 

s3- 

s3* 

sO 

CM 



1 

«0 

H 

CM 

to 

to 

3 

3 

to 

to 

CM 

O 


to 

O' 



E l 


• 

* 

♦ 

• 

• 

• 

• 

• 


• 


• 

• 



1 

1 

to 

3 

3 

3 

S3* 

a3- 

A3- 

3 

stf- 



to 

to 

CM 



1 

» 

3 

m 

CO 

O' 

O' 

aO 

CM 

in 

sO 

© 

C0 


co 

h* 



1 

•0 


CM 

to 

S3- 

m 

sO 

sO 

sO 

sO 

in 

3 

H 

CO 



1 

CM 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

Oil 



1 

1 

to 

to 

to 

CM 

in 

so 

in 

CM 

p- 

^•4 

CM 

m 

H 

O' 



1 

© 

3 

in 

p. 

to 

O' 

o 

r* 

F«4 

CM 

IM 

H 

© 

PH 



"> 1 

• 

* 

• 

♦ 

• 

* 

* 

# 

* 

c 


c 

* 

* 



1 

1 

CM 

CM 

CM 

CM 

CM 

CM 

to 

to 

to 

to 

to 

to 

to 

CM 


>- 

o 

© 

o 

© 

G 

© 

o 

© 

0 

© 

© 

o 

© 

© 

Qt -i 

c 

• 

* 

» 

* 

* 

« 

* 

* 

• 

« 

• 

• 

• 

OCM 

© 

o 

m 

o 


© 

in 

© 

in 

O 

© 

© 

© 

o 

< t- 


H 

c—4 

CM 

CM 

to 

to 

3 

3 

in 

sO 

fH. 

to 

O' 


O' 

CM 

I 

O 


I 




AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 




N 

to 

9* 

* 

O 

m 

fv 

to 

to 

9 


N 


♦ 

< Z 

» 


* 

• 

t 

• 

* 



• 

* 

• 

• 

• 

S3 

in 

•0 

Ifl 

s© 

o 


in 

fv 

N 


H 

m 

s© 

CM 

x© 

m 

<M 

in 

f*s 

9* 



to 

CV 

m 

9* 

e 

O' 

fN 

ZC 

pH 

CM 

cm 

CM 

N 

CM 

CM 

<M 

m 

CM 


PH 

9 

to 

< < 

H 

PH 

H 

H 

PH 

PH 

H 

PH 

pH 

H 

pH 

PH 






i 

t 

© 

o 

a 

9 

FV 

m 

CM 

to 

to 

tv 

CM 


CM 

<0 



l 

fv 

9 

o 

© 

pH 

CM 

to 

to 

*■ 

s3* 

in 

in 

m 

s3 



a i 

t 

• 

♦ 

• 

t 

• 

• 

♦ 

• 

• 

• 

t 

t 

• 



1 

1 

O 

© 

pH 


pH 

H 

pH 

pH 

pH 

pH 

pH 

H 

pH 

PH 



1 

1 

fv 

9 

9 

9 

IV 

m 

pH 

rv 

pH 

m 

9 

00 

CM 

to 



i 

0 

CM 

to 

s3 

in 


fv 

rv 

to 

•0 

m 

•0 

00 

rv 



X t 

* 

t 

t 

* 

t 

• 

• 

• 




‘ 

t 

• 



1 

I 

H 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

PH 

pH 

pH 



1 

1 

pH 

pH 

in 

rv 

oo 

tv 

m 

pH 

m 

to 

00 

© 

SO 

s0 



1 

O 

tO 

*■ 

m 


rv 

oo 

9 

9 

9 

9 

9 

tv 

m 



O 1 

• 

t 

• 

• 

t 

• 


* 

* 

• 

• 

• 

• 

* 



1 

1 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 



1 

i 

in 

CM 

rv 

o 

pH 

9 

S3* 

rv 

tv 

s3* 

pH 

rv 

Mf 

CM 



1 

S0 

O 

pH 

to 

*3* 

*3* 

m 

m 

m 

in 

s3 

pH 

•0 

S3 



in i 

t 

• 

* 

* 

t 

. 

• 

• 

t 

• 

* 

t 

• 

* 



i 

K» 

s3 


S3* 

S3- 

s3- 

s3* 

s3* 

<3* 

s3" 

s3 

9 

to 

to 

in 


» 















< UJ 

X 

i 















3i- 

►“ 

i 















=> 

X 

» 

00 

m 

CM 

rv 

9 

9 

m 

9 

O 

00 

O 

in 

m 

to 

X 

o 

i 

00 

pH 

CM 

cm 

CM 

CM 

CM 

pH 

pH 

9 

sO 

CM 

rv 

pH 

►H 

E 

< i 

• 

t 

* 

• 

• 

t 

« 

* 

« 

• 

• 

• 

# 

* 

E 


i 

•9 

in 

in 

in 

m 

in 

in 

m 

m 

s3- 

s3 

s3 

to 

to 


>“ 

i 















*0 

0 

i 















m 


i 
















in 

i 

%3 

rv 

9 

tv 

CM 

■3- 

CM 

tv 

9 

9 

CM 

s3 

9 

CM 

in 

p- 

i 

o 

iH 

pH 

pH 

rH 

© 

9 

tv 

m 

to 

9 

to 

sO 

o 

UJ 

X 

~9 1 

• 

* 

t 

• 

* 

* 

* 

• 

♦ 

♦ 

* 


# 

* 

UJ 

D 

1 

so 

SO 

sO 

>0 

sO 

sO 

in 

in 

m 

in 

s3* 

s3 

to 

to 

< <* 

O 

I 















►H O 

E 

1 















O. UJ 

< 

1 















r q 


1 

S3 

9 

rv 

to 

rv 

tv 

in 

o 

to 

S3 1 

CM 

CM 

fv 

CM 

>- 

>- 

1 

tO 

to 

to 

to 

CM 

pH 

o 

9 

r* 

in 

pH 

•*> 

o 

m 

—1 «0 

mJ 

—) 1 

• 

♦ 

* 

t 

• 

t 

• 

• 

• 

• 

• 

• 

t 

« 

O *3 

PH 

» 

m 

»n 

in 

m 

in 

m 

in 

S3* 

S3* 

s3- 

S3 

to 

to 

CM 


< 

1 















»• *. 

Q 

1 















UJ UJ 


1 















t- a 

UJ 

1 

*3 

m 

sO 

s3 

o 

CM 

to 

O 

in 

00 

9 

CM 

9 

s3 


O 

t 

H 

CM 

CM 

CM 

CM 

pH 

o 

9 

tv 

in 

pH 

IV 

pH 

sO 

m 

< 

E i 

♦ 

t 


# 


* 

• 

• 

* 

* 

* 

• 

• 

t 

►H 

or 

1 

m 

in 

m 

m 

in 

m 

in 

S3- 

S3* 

S3* 

s3 

to 

to 

CM 

t— 

UJ 

1 















< 

> 

1 















-ml 

< 

1 

















1 

s© 

to 

9 

CM 

to 

CM 

•0 

to 

m 

m 

© 

00 

9 

m 



I 

9 

pH 

pH 

CM 

CM 


•H 

pH 

o 

9 

tv 

to 

9 

in 



< 1 

* 

• 

* 

* 

» 

* 

• 

• 

* 

• 

• 

» 

t 

t 



1 

» 

tO 

s3 

s3* 

»3* 

s3“ 

S3* 

S3 

S3" 

s3 

to 

to 

to 

CM 

CM 



1 

1 

sO 

rv 

sO 

to 

9 

to 

in 

SO 

in 

CM 

pH 

s3 

pH 

s3 



1 

sO 

00 

9 

o 

© 

pH 

pH 

pH 


pH 

© 

oo 

sO 

to 



E i 

• 


♦ 

« 

• 

t 

• 

• 

. 

« 

• 



t 



1 

1 

CM 

CM 

CM 

to 

to 

to 

to 

to 

to 

to 

to 

CM 

CM 

CM 



1 

1 

00 

fv 

in 

CM 

00 

to 

fv 

9 

pH 

CM 

9 

CM 

pH 

sO 



1 

in 

fv 

00 

9 

9 

o 

© 

O 

pH 

pH 

© 

© 

9 

fv 



U. ! 

• 

t 


• 

• 

* 

* 

* 

• 

• 


• 

• 

• 



1 

1 

rt 

pH 

pH 

pH 

pH 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

pH 

•H 



t 

1 

m 

m 

s3 

CM 

© 

tv 

to 

«0 

to 

sO 

© 

O 

sO 

00 



1 

•0 

o 

pH 

CM 

to 

to 

s3 

<3- 

in 

m 

Srf) 

sO 

in 

s3 



“> 1 

• 

» 

• 

♦ 

• 

t 

* 


• 

• 

• 

• 

* 

• 



1 

1 

o 

pH 

pH 


pH 

rH 

pH 

pH 

pH 



pH 

pH 



>* 

►- oooooooooooooo 

CKlmJ 

qc *-* ooinoinoinoin©oooo 
<h ph ph cm fsi »o to »3 s3 m «o rv «o 9 


o 


H 


I 'i 

i 1 

I 1 




4 


122 - 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/S9. M) 



average daily total terrestrial insolation ESTIMATES (RUHsSQ 

. .SITE: SPOKANE WA 

LATITUDE: 47 DEGREES 38 MINUTES 




AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 



< ui 
2»- 


E 

O 

< (St 

-j 

mm 
•HlU 
Ui 
> t*c 
tu o 
to UJ 

a a 
►-» 

z«o 

3 


UJ w 
f-O 


Ui 

o 

ac < 
UJfi 
> 

< 


<z 

33 

20 

ZE 


2 

O 

£ < 

>- 

© 

4/> 


3 

O 

E 


> 

PM 

< 


UJ 

o 


> 

< 



**> 

* 

fH 


tO 

p 

« 

« 

« 

fN 

m 

to 

© 

© 


» 

* 

• 

• 

• 

» 

• 

• 

* 

• 

* 

* 

t 

• 


tO 

m 

tO 

m 

m 

m 

m 

in 

in 

to 

Kl 

to 

to 

fy 


r- 

(SI 

. 

f> 

p 

o 

o 

n* 

40 

o 

in 

to 

© 

(SI 


m 

M 

in 

• 

« 

o 

(SI 


to 

in 

m 

© 

© 

to 


#■* 

O 

m 

© 

f* 

« 

o 

• 

K 

m 

o 

© 

© 

m 



»n 

m 

tn 

in 

»n 

in 

m 

in 

to 

N 

ry 

© 

© 



H 

pm 

pm 


p< 


PM 

PM 

PM 

PM 

PM 

PHI 


1 

1 

m 

N 

o 

•0 

40 

m 

o 

o 

«0 

(M 

fN 

to 

m 

© 

1 

f* 

O' 

° 


P4 

cy 

M 

m 

tn 

O 

O 

© 

© 

in 

1 

1 

e 

o 

pM 


PM 

PM 

PM 


pm 

pM 

PM 

PM 

PM 

PM 

1 

1 

fy 

© 

m 

m 

40 

•o 

m 

m 

O 

m 

pM 

pH 

fN 

© 

1 

H 

© 

m 

40 

fV 

«0 

o 

e 

° 

PM 

(SI 

fy 

PM 

© 

I 

1 

PM 

H 

pm 

PM 

pM 

p* 


<M 

rj 

fJ 

<v 

ry 

fy 

fy 

1 

© 

(SI 

to 

fy 

m 

40 

in 

in 

«0 

(SI 

o 

f^ 

fn 

fy 

* 

© 

♦ 

in 

k. 

tO 

o 

o 

pM 

H 

fV 

fy 

pM 

© 

to 

1 

1 

(SI 

fy 

<Si 

fy 

IM 

N 

in 

K) 

in 

to 

m 

90 

m 

fy 

* 

1 

o 


W 

m 

m 

o 

i"i 

m 

40 

o 

fsi 

© 

© 

to 

1 

n* 

PM 

fsi 

O' 

m 

40 


n. 

h> 

n* 

© 

90 

© 

© 

1 

1 


© 

© 

o 

o 

o 

o 

o 

O 

o 

o 

© 

© 

m 

1 

1 

fy 

H 

c» 

fN 

P* 

p* 

o 

in 

O 

fy 

pM 

© 


© 

1 

o 

K> 

© 

o 

in 

m 

o 

o 

m 

CJ 

O 

© 

© 

m 

1 

1 

in 

in 

m 

m 

in 

m 

m 

m 

in 

m 

O 

© 

to 

m 

I 

1 

in 

O 


K> 

p 


o 

o 

40 

in 

fN 

© 

PM 

pM 

1 

Ol 

© 

© 

o 

in 

in 

pM 

G 

*0 

40 

p4 

m 

© 

fy 

1 

1 

40 

© 

40 

40 

40 

o 

40 

40 

m 

m 

in 

© 

to 

to 

1 

1 

tO 

© 

© 

in 

© 

o 


pM 

K> 

to 

to 

in 

in 

m 

1 

fN. 

•0 

fN 

fN 

40 

m 

o 

k> 

pM 

o 

o 

© 

tn 

fN 

I 

1 

m 

in 

m 

m 

in 

m 

m 

m 

in 

o 

o 

90 

in 

fsi 

1 

1 

in 

fN 

© 

© 

in 

•0 

« 

in 

o 

pM 

o 

© 

© 

0 

1 

m 

© 

40 

40 

40 

m 

o 

k» 

rt 

© 

m 

© 

m 

© 

t 

1 

IA 


in 

m 

in 

m 

m 

in 

m 

in 

o 

© 

m 

fy 

1 

J 

© 

m 

cs t 

n* 

© 

o 

40 

o 

in 

90 

© 

fy 

© 

fj 

t 

N 

© 

in 

in 

m 

m 

in 

m 

o 

to 

© 

fN 

in 

to 

1 

1 


© 

O' 

o 

o 

o 

p 

o 

o 

p 

© 

w 

in 

fy 

t 

1 

© 

© 

* 

n. 

in 

PM 

o 

40 

40 

o 

© 

© 

pM 

© 

1 

tO 

pm 

»• 

K> 

o 

m 

in 

m 

m 

m 

© 

(SI 

© 


I 

* 

ISJ 

m 

m 

in 

K) 

in 

in 

m 

in 

in 

in 

in 

90 

fy 

1 

•0 

pM 

pm 

pm 

© 

40 

^4 

40 

o 

o 

© 

in 

fy 

in 

1 

© 

O' 

o 


PM 

rsi 

m 

in 

iO 

o 

m 

in 

ry 

© 

1 

1 


pM 

ISI 

N 

rj 

(SI 

<M 

m 

(si 

m 

ry 

ry 

ry 

fy 

1 

I 

O' 

© 

m 

O 

40 

m 

in 

e 

40 

pM 

fN 

«o 

© 

© 

t 

«o 

*-H 

INI 

in 

o 

in 

o 



«0 

to 

to 

•0 

fs 

1 

1 

e 

PM 

pm 

p< 

*M 

pM 

pM 

H 

pM 

*+ 

“ 

pM 

pH 

“ 



© 

o 

o 


o 

O 

e 

o 

© 

© 

© 

o 

O 

-J 

• 

• 

« 

* 

» 

• 

• 

• 

* 

* 

* 

t 

t 

• 

PM 

o 

o 

in 

o 

m 

O 

in 

e 

in 

O 

© 


© 

© 

f» 


pM 

pm 

fy 

M 

in 

in 

o 

o 

m 

© 

i‘ 

to 

© 


o 


i 


■fr -ri.A 


-224 


AVERAGE DA II V TOTAL TERRESTRIAL INSCIATIU>< ESTIMATES (KUH/SQ. M) 



Ui 


< >- 

• 

© 

tn 

* 

© 


rs> 

© 

© 

m 

K> 

© 

© 

N 

at< 

* 

• 

* 

• 

• 

t 

• 

> 

• 

• 

* 

• 


• 

UiQ 


© 

♦ 

© 

© 

♦ 

© 

© 

© 

© 

© 

♦ 

r» 




• 

O 

e 

e 

tn 

© 

©» 

ph 

K» 

© 

© 

m 

N 

<0 

<Z 

• 

* 

• 

• 

• 

* 

* 

* 

* 

* 

* 

* 


• 

33 

» 

« 

ph 

aT 

f** 

• 

90 

f*. 

A 

m 

* 

m 

© 

© 

zo 

f-N 

© 

tn 

© 

•e 

O 

© 

; O' 

<e 

m 

n* 

© 

m 

*A 

zz 

© 

m 

© 

© 

© 

fH 

n* 

© 

© 

© 

m 

© 

to 

ph 

< < 

ph 

pn 

ph 

PH 

pH 


P< 

R"H 


PH 

pH 

*4 


ph 




I 

* 

HI 

n* 

to 

•0 

« 

<0 

«r 



RH 

K 

« 

<r 




i 

p t 

* 

PH 

ti 

RO 

o 

m 

«o 


n» 

•0 

•0 

« 

« 

(N 



i 

i 

o 

pH 

pH 

pH 


pH 

PH 



pH 

PH 

PH 

pH 

pH 



i 

i 

19 

m 

(SI 

rs. 

o 

to 

O’ 

nr 

M 

« 

m 

nO 

O' 

© 



i 

* 

c* 

O' 

o 

(SI 

to 

nr 

m 

nO 

© 

n* 

fH 

* 

tn 



i 

i 

PM 

pH 

PI 

N 

(SI 

(Si 

fw 

(SI 

(M 

cvl 

rJ 

(SI 

(SI 

CM 



i 

i 

PH 

© 

•0 

© 

■0 

m 

O' 

p^ 

O' 

in 

to 

* 

O' 

«» 



i 

« 

K> 

m 

K 

O' 

pH 

(SI 

< 

nr 

m 

m 

o 

(SI 

O' 



C 1 

* 

♦ 

• 

• 

• 

• 

* 



• 

• 



* 



1 

1 

(SI 

in 


to 


O 

o 

nr 

nr 

nr 

nr 

nr 

nr 

tn 



1 

1 

to 

© 

•0 

to 

m 

•5 

fn. 


«r 

(SI 

© 

nO 

pH 

© 



1 

© 

(SI 

HT 

HO 

•0 

O' 

«v 

pH 

H 

H 

O' 

nO 

(SI 

n» 



90 t 

• 

• 

• 

• 


• 

* 

• 

* 

• 

• 






i 

*r 

m 

in 

m 

IT, 

aO 

HB 

o 

nO 

nO 

m 

m 

tn 

© 

V* 


i 















<* U4 

X 

i 















at»~- 

t- 

i 















3 

z 

i 

K> 

pH 


Ml 

n* 

to 

in 


tn 

o 

O' 

»n 

m 



o 

i 

N 

© 


•0 

to 

to 

to 


nO 

m 

o 

in 

•0 

© 

►M 

c 

< i 

• 

* 

• 

* 

t 

t 

* 

• 

• 

• 

• 

* 


• 

z 


i 

© 


s0 

HD 

HD 

nO 

A0 

nO 

o 

nO 

nO 

in 

«r 

© 


S- 

i 















© 

U> 

i 















to 


i 
















© 

» 


N 

in 

O 

•s 

«0 

to 

O’ 

© 

«n 

nr 

© 

nr 

n. 

tn 

p* 

i 

* 

<0 

<0 

H& 

in 

O 

to 

pH 

o* 

nO 

o 

(Si 

nr 

us 

Ui 

z 

*» i 

* 

• 

« 


« 


t 

* 

• 


• 




Ui 

3 

i 

fH 

n* 


n*. 

fH. 

n» 

fH. 

r*. 

«o 

•o 

nO 

tn 


RO 

tt 

o 

i 















< o 

z 

i 















Zuj 

«c 

i 















Nfi 


i 

to 


© 

O 1 

m 

pn 

© 

nO 

Rn 

to 

pH 

RO 

O' 

© 

*£ 

► 

i 

to 

e* 

p* 


K 

© 

nr 


o 

fs. 

(Si 

m 


o 

Pf W> 

-J 

1 i 

• 

♦ 

* 

• 

t 


• 

• 

« 

* 


» 

« 

* 

>• ^ 

+4 

1 

© 

© 

<0 

<0 

<0 

nO 

nO 

<o 

nO 

m 

m 

o 

tn 



< 

l 















M «• 

a 

1 















Ui Ui 


l 















►- o 

Ui 

l 

*n 


© 

•5 

RO 

O 

pH 

in 

nO 

in 

tn 

nr 

< 

© 

«3 

o 

1 

to 

© 

*■ 


o 

•n 

(Si 

© 

to 

nO 

pH 

m 

© 

pH 

tn »- 

< 

Z i 

» 

* 

* 

• 

* 

* 

• 

• 


* 

• 

• 


« 

*-h 

oc 

t 


© 

«0 

*0 

• P 

o 

© 

nO 

m 

m 

m 

nr 

to 

RO 

H 1 

Ui 

I 















<Z 

> 

1 















•i 

«r 

1 

















I 

© 

o 

•0 


«0 

m 

(SI 

m 

m 

to 

o 

n. 

m 

© 



1 

O 

K> 

•n 

O’ 

o 

•o 

■r 

RO 

nj 

pH 

© 

»n 

© 

(Si 



© \ 

* 

♦ 

« 

* 

* 

• 

• 

• 

• 

« 

* 

• 

« 

t 



1 

1 

US 

m 

m 

m 

m 

in 

in 

m 

in 

«n 

o 

o 

RO 

to 



1 

i 

© 

o 


nO 

m 

K> 

M 

© 

© 

fs. 

Hf 

pH 

o 

o 



i 

in 

•0 

© 

pH 

w 

to 

K) 

nr 

«r 

to 

(ni 

O 

© 

(Si 



I 

1 

MS 

RO 


o 

o 

< 

O’ 

nr 

nr 

nr 

*r 

* 

© 




i 

1 

© 

tt 

*& 

to 

© 

to 

nO 

(M 

nO 

O' 

•0 

© 

in 

R"S 



l 

pH 

© 

in 

N& 

n. 

•0 

O' 

© 

© 

© 

© 

o 

© 

© 



U. 1 

♦ 

• 

♦ 

* 

t 


• 

* 

• 

• 

* 

* 





l 

1 

(M 

M 

M 

(SI 

w 

(Si 

(SI 

RO 

to 

to 

in 

to 

(si 

«v 



i 

l 

in 

© 

pH 

pH 

© 

«0 

nO 

(si 

r» 

pH 

m 

nr 

« 




l 

pH 

<r 

in 


K 

n. 

to 

O' 

O' 

o 

© 

o 

O' 

© 



-"> l 

• 

♦ 

• 

* 

« 

• 

• 

• 

* 

« 

« 

■ 

* 




1 

I 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 


(Si 

N 

pH 

pH 



C-225 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES CKMH/SQ. M) 



-226 


average daily total terrestrial insolation estimated (KWH/SQ. m> 

SITE: GREEN BAY WI 

LATITUDE: AA DEGREES 29 MINUTES 



i- - w .v^MMKnCUWikf^ 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. H) 



average daily total terrestrial insolation estimates ocuh/sq. M) 

SITE: MADISON UII 

LATITUDE: 43 DEGREES 8 MINUTES 



AVERAGE DAILY TOTAL TERREST 



AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. 

SITE: CHARLESTON WV 

LATITUDE: 38 DEGREES 22 MINUTES 



> 

< e 
fie -j 

fi^H O 


o O 
© in 




m 

C\J 


© 

K) 


in 


in 


m 



«o & 




n ^ e e e 

^ ^ 


UJ 

O 

<>■ 
DC < 
UJO 
> 

< 


^'JtflOWKtKIrtN 


«z 

zo 
ZE 
« «* 


60 

sa- 

h* 

CM 

m 

60 

O 

m 

ps. 

pH 

CO 

sO 

sO 

Os 

m 

CM 

in 

* o 

r*v 

fs. 

sO 

ro 

pH 

a 

pH 

ST 

pH 

a 

pH 

a 

pH 

a 

pH 

a 

pH 



O' 

in 

\ 

X 

2 


in 

LU 

K 

< 


in 

LU 


o 

M 

>- 

< 

mJ 

o 

in 

z 


m 

> LU 

3h 

X 


CM 

CM 


in 

Ui 

a: 

ft: 

LU 


ZO 
hh LU 
H-O 

z 

©oo 
Z K) 


— > LULU 

< H- O 
>- WD 

o in k- 

»“ M 

►- 

>* < 


< 

© 

LU 

O 

< 

ft£ 

LU 

> 

< 


I 

I 

I 

© I 
I 
I 
I 
I 
i 
t 

z > 

i 

t 

i 

i 

i 


sa 


sa 

r*. 


sD 

00 


O' 


o 

CM 


.16 

a 

CM 

.30 

.35 

.38 

pH 

a 

60 

ro 

29 

15 

CM 

CM 

CM 

CM 

CM 

CM 

eg 

CM 

CM 


10 • 

.36 

CM 

in 

.67 

o 

60 

.91 

.00 

so 

a 

PH 

60 

pH 

O 

CM 

sa- 

i-4 

01 

o 

60 

CM 

CM 

CM 

CM 

CM 

CM 

to 

to 

to 

to 

ro 

ro 

rO 

CM 


O 

z 

>- 

© 


t 

i 

i 

sa- 

a 

sa- 

a 

a 

•a- 

a 

sa- 

sa- 

a 

a 

a 

ro 

ro 

i 

i 

i 

■< j 

.98 

a 

pH 

.17 

.18 

.15 

.10 

.03 

92 

79 

a 

o 

sD 

CM 

o 

60 

LZ 

68 

i 

i 

> 

a 

in 

in 

m 

m 

m 

m 

sa- 

sa- 

s3- 

a* 

ro 

ro 

CM 


r 

< 

>- 


< 

© 


o 

<£ 

CZ 

LU 

> 


in 


K) *•> 


•a* 

r>- 


a* 

eo 


cm ro i-t 

o *-• CM 


pH rj 

sa m 


p*. 

so 


O pH 

r-. 


SO 

CM 


Os 

SO 


Os 

CM 


sO 


Os 

CM 


f-. 

m 


rO 

CM 


in 

tO 


Os 

r«. 


a- 

a* 


fO 


ro k) 


to 

so 


-J 

z m 

m 

.58 

pH 

60 

CM 

ro 

pH 

S0 

Cs 

< 

©LU 

z - 

s© 

in 

m 

a 

to 

pH 

OS 

►H 

ft: 

H» LU 
o ft: 

© 

o 

1 in 

i 

in 

m 

m 

m 

in 

m 

sa- 


oo 

stf- 


O pH 
SO pH 


in 

m 


so 

Os 


m 

ro 


CM CM 


CM 

60 

.80 

a 

.66 

sa- 

in 

o 

sa 

23 

03 

CM 

oo 

59 

m 

in 

in 

m 

m 

in 

IT) 

m 

sa 

a- 


SO 

•a 


so 

00 


m 

CM 


n M cm 


Os 

K> 


CM 

a 


Os 

rO 


ro 

ro 


m m in in m 


sa 

o 

ro 

sO 

m 

o 

in 

60 

D 

pH 

c» 

60 

HD 

a- 

OS 

sa 

60 

to 

in 

m 

sa 

sa 

>a 

ro 

to 

CM 

CM 


si 


s a 


ro 


ro 

rv 


n- 


sO 

sO 


Os 

m 


m 

a 


co 

rO 


stf- 

CM 


Os 

ro 


60 

stf* 


60 

a 


m ro cm 


I 

I 

C f 

l 

I 

l 

i 

I 

i 

u. t 

i 

i 

i 

i 

i 

i 

-> i 
t 
i 


.36 

.58 

.65 

.71 

.75 

.77 

.77 

75 

71 

65 

PH. 

a 

CM 

CM 

06 

53 

ro 

ro 

rO 

ro 

ro 

rO 

ro 

ro 

ro 

to 

ro 

to 

CM 

CM 

.39 

.63 

.73 

CM 

60 

.89 

.95 

. 99 

.01 

02 

01 

»a 

O' 

© 

60 

61 

to 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

ro 

to 

to 

CM 

CM 

CM 

CM 

.65 

.90 

.01 

.11 

.19 

.27 

.33 

60 

ro 

pH 

sa 

ro 

a 

ro 

sa 

37 

sO 

CM 

10 

pH 

pH 

CM 

CM 

CM 

CM 

CM 

CM 

cm 

<M 

CM 

CM 

CM 

CM 


CM 

ro 

CM 

l 

o 


>- 

<h- 
0 C-J 
ft:»-H 
< h- 


m 


eg 


m 

CM 


in 

rO 


o 

a 


m 

a 


in 




AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SQ. M) 


l 

4 


Ul 

o 















<>- 

on 

4 

1 n 

N 

•0 

10 

to 

«o 

to 

p^ 

m 

#H 

N* 

e 

oe< 

• 

• 

• 

* 

» 

• 

t 

• 

* 

• 

• 

• 

t 

• 

uJ£j 

> 

4 

♦ 

tf\ 

m 


m 

in 

m 

m 

m 

in 

m 

m 

«■ 

♦ 


_J »- tONNONtttOrN — OtvfOrs 


<z 

t 

t 

• 

• 

• 

« 

• 

• 

• 

* 

• 

t 

« 

* 

S3 

tO 

•0 

•0 

nO 

rH 

m 

W 

•0 

tO 

nO 


to 

M 

n* 

ZO 

o 

in 

rH 

n© 

© 

CM 

tO 

CM 

rH 

•0 

On 

m 

fN. 

in 

zr 

•0 

On 

O 

O 

rH 

rH 


rH 

H 

o 

O' 

oo 

NO 

nt 

4< 

fH 

H 

IV 

CM 

CM 

CM 

CM 

CM 

IM 

cm 

H 

rH 

H 

H 




1 

1 

n*. 

in 

CM 

n. 

o 

N 

0 

hN 

rH 

CM 

Cs 

rH 

to 

*r 



1 

to 

nT 

r* 


CM 

<T 

sO 

Cn» 

Os 

o 

rH 

M 

rH 

O' 



O 1 


• 

• 

• 

t 

* 

• 

• 

t 

• 

• 

• 

* 

* 



1 

1 

H 

CM 

CM 

CM 

to 

to 

to 

to 

to 

*3“ 


nT 

NT 

to 



1 

» 

N 

to 

O' 

n- 

to 

NO 

h* 

m 

O' 

rH 

m 

in 

<\l 

so 



1 

nt 

O 

»o 

SO 

O' 

rH 

to 

m 

so 

00 

O' 

O' 

to 

in 



Z 1 

« 

t 

• 


• 

* 

• 

» 

* 

• 

* 

* 

• 

• 



1 

1 

«SI 

tO 

to 

to 

to 

NT 


S& 

St 

St 


NT 

NT 

nT 



1 

1 

NT 

sO 

to 


to 

m 

nJ* 

to 

Os 

NO 

n* 

CM 

CM 

nO 



t 

to 

in 

to 

rH 

O' 

sO 

to 

O' 

o 

rH 

rH 

O 

n- 

CM 



O 1 


* 


t 

• 

t 


• 

* 

♦ 

• 

* 

• 

* 



1 

1 

tO 


Ai* 

m 

m 

m 

m 

m 

so 

nO 

NO 

NO 

m 

m 



1 

1 

CM 

*T 

O' 

rH 

to 

rH 

o 

st 

to 

00 

in 

m 

On 

to 



1 

m 

rH 

to 

NO 


On 

o 

o 

o 

O' 


to 

Ps. 

e 



in i 

• 

t 

* 

• 

• 

* 

• 

• 

• 

* 


• 

• 

• 



i 

m 

sO 

nO 

sO 

sO 

nO 

rs. 

n- 

n* 

sO 

nO 

NO 

in 

m 

*n 


i 















>- UJ 

X 

i 















31- 

1- 

i 















=0 

Z 

i 

CM 

© 

CM 

© 

to 

rH 

Nrf- 

to 

00 

00 

P's 

ON 

00 

NO 

z 

o 

i 

o 

nT 

in 

>e 

NO 

NO 

in 

Nfl- 

CM 

o 

in 

to 

o 

rH 

H 

E 

< i 

• 

♦ 

• 

t 

t 

* 

♦ 

• 

• 

• 

* 

* 

• 

• 

E 


t 


n- 

r-v 

IS. 

p^ 

Ps. 

r** 

N 

rv 


NO 

m 

in 



>- 

i 















m 

CO 

i 















in 


i 
















<n 

i 

O' 

in 

m 

rH 

rH 

rv 

to 

NT 

n* 

r*v 

sO 

rH 

PH 

CM 

in 

»— 

t 

o* 

rH 

rH 

rH 

O 

oo 

sO 


rH 

to 

rH 

to 

to 

>T 

LU 

z 

-> i 

t 

• 

t 

* 

• 


• 

• 

• 

• 

t 

* 

• 

• 

UJ 

r> 

» 

P** 

00 

00 

oo 

OO 

r*N 

hv 


n* 

NO 

NO 

in 

NT 

to 


o 

i 















o: o 

E 

i 















LU UJ 

< 

i 















Q.Q 


i 

«0 

tO 

Ch 

O 

n. 

o 

to 

H 

CM 

rH 

Os 

nT 

NO 

nT 

m 

>- 

i 

00 

O' 

00 

00 

sO 

m 

CM 

O 

r>» 

n3* 

sO 

oo 

O 

o 

C CM 

-J 

-> i 

• 

• 

• 


• 

♦ 

* 

* 

* 

• 

• 

• 

* 

• 

O <T 

t-H 

i 

n* 

r% 

r*» 

h* 

r>. 

Ch 

rN. 

h- 

sO 

SO 

in 

NT 

to 

to 



i 















•• •• 

Q 

i 















UJ UJ 


i 















t-Q 

UJ 

i 


>o 

sO 

CM 


CM 

NO 

NO 

to 

90 

00 

«0 

O' 

On 

tM ID 

O 

i 

O' 

© 

© 

o 

o 

00 

NO 

Nj- 

CM 

O' 

to 

NO 

oo 

O 

in 

< 

E i 

• 

• 

* 

• 

* 

• 

t 

• 

• 

* 

• 

• 

• 

• 

t-H 

0£ 

1 

*0 


P's 

r-. 

sO 

NO 

NO 

N© 

NO 

in 

m 

NT 

to 

to 

K 

UJ 

1 















< 

> 

1 















-J 

< 

t 

















1 

to 

© 

00 

to 

t 

rH 

in 

in 

to 


in 

CM 

00 

NO 



1 

00 

rH 

rH 

CM 

CM 

CM 

rH 

o 

O' 

n. 

to 

00 

rH 

nT 



< I 

♦ 

* 


# 

• 

t 

• 

• 

* 

♦ 

* 

• 

• 

* 



1 

1 

m 

nD 


SO 

NO 

sO 

NO 

NO 

in 

in 

in 

St 


to 



1 

1 

•t 

O' 

rv 

CM 

Ntf- 

to 

00 

o 

O' 

Nd- 

st* 

CM 

00 

to 



1 

in 

O' 

H 

to 

nT 

in 

m 

nO 

m 

m 

to 

o 

in 

o 



E 1 

• 

* 

t 

♦ 

t 

» 

• 

• 

• 

• 

* 

* 

• 

• 



1 

1 


nT 

m 

in 

in 

m 

m 

in 

in 

m 

in 

in 

•T 

t 



1 

1 

O' 

rH 

nT 

in 

to 

00 

o 

o 

fv 

o 

O' 

in 

e 

NT 



1 

rH 

n* 

O' 

rH 

K> 

nT 

NO 


n* 

00 

C*. 

NO 

nT 

© 



U. 1 

• 

# 

• 

* 

* 

♦ 

• 

• 

* 

• 

* 

t 

t 

• 



J 

1 

►o 

IO 

to 


NT 

st 

'tf' 

st 

nJ* 


NJ* 

‘T 

*T 

nT 



1 

1 

sO 

o 

in 

«0 

o 

00 

nJ* 

O' 

o 

O' 

O' 

O' 

PN* 

nT 



1 

rH 

n* 

O' 

rH 

to 

in 

n* 

OO 

o 

e 

rH 

rH 

o 

10 



“> * 

* 

• 

• 

t 

t 

• 

» 

• 

» 

• 

• 

♦ 

• 

• 



1 

1 

CM 

CM 

CM 

»o 

to 

to 

to 

rO 

St 


NT 

NT 

nT 

to 


>- 

<i- oooooooooooooo 

a 

ai w eomomotneinoeooo 

< ►— <Hrlt'JfS|IOK5s}-C , LO»Or-'«0O > 


o 




-233 


AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KUH/SQ. M) 



UJ 

o 

<>- 
a t< 

K 

M 

CSI 

d 

1A 

IA 

IA 

m 

IA 

d 

CM 

•0 

to 

•0 

UJO 

> 

d 

IA 

in 

lA 

IA 

IA 

IA 

IA 

IA 

IA 

IA 


d 

to 

< 
















in 
>UJ 
3 •— 
s 


E 

o' 


in 

UJ 
LU ill 
2 CkC 

z o 

UJ UJ 
>- Q 


Xph 

O d 


UJ UJ 
KO 
M3 

in m 


-J M 

< z 


fV 

•0 

d 

to 

lA 

CM 

O' 

rO 

CM 

O' 

O' 

pH 

00 

3D 

rv 

O' 

d 

00 

o 

o 

•0 

90 

O' 

K) 

lA 

♦A 

IA 

to 

Z O 

M 

IA 

M 

IA 

O' 


H 

N 

O' 

fv 

00 

IA 

00 

00 

z e 

A* 

•0 

O' 

O' 

O' 

o 

O 

O 

O' 

O' 

00 

rv 

IA 

K) 

< 

pH 

PH 

M 

M 

M 

CM 

CM 

CM 

M 

pH 

M 

H 

pH 

pH 



i 

i 

N 

d 

CM 

O' 


nO 

IA 

90 

fv 

o* 

90 

in 

m 

IA 


i 

O 1 

M 

fv 

O 

CM 

Ul 

fv 

O' 

pH 

CM 

IA 

m 

m 

d 

CM 


1 

1 

1 

CM 

CM 

IA 

IA 

IA 

IA 

IA 

d 

d 

sJ* 

d 

d 

d 

d 


1 

1 

e 

sO 

»0 

d 

O 

CM 

f\J 

O' 

IA 

d 

nD 

d 

00 

© 


1 

Z 1 

<0 

CM 

IA 

00 

pH 

IA 

IA 

NO 

00 

O' 

o 

o 

<0 

*o 


1 

1 

1 

CM 

IA 

IA 

IA 

d 

M- 

d 

NT 

d 

Nd 

m 

in 

d 

d 


1 

1 

H 

O' 

O' 

M> 

O' 

O' 

NO 

O' 

fv 

CM 

pH 

IA 




1 

O 1 

O* 

IA 

00 

M 

IA 

lA 

fv 

00 

O' 

O 

O 

00 

m 

o 


1 

1 

1 

K1 

d 

d 

in 

IA 

lA 

lA 

IA 

m 

NO 

'O 

m 

in 

m 


1 

1 

no 

•0 

O' 

|v 

O 

© 

NO 

•0 

m 

O' 

d 

d 

pH 

in 


1 

(A 1 

CM 

rv 

O' 

H 

lA 

d 

d 

Ntf* 

d 

IA 

pH 

rv 

CM 

m 


1 

in 

IA 

IA 


NO 

NO 

NO 

NO 

nO 

NO 

'O 

in 

m 

d 

X 

1 















M 

1 















z 

1 

o 

00 

'O 

© 

© 

M> 

O' 

fv 

IA 

in 

00 

O' 

© 


o 

E 

1 

< 1 

CM 

d 

lA 

<sD 

nO 

IA 

d 

IA 

CM 

o 

m 

O' 

IA 

m 

V 

1 

nO 

*o 

O 

nO 

nO 

sO 

nO 

NO 

NO 

nO 

in 

d 

d 

90 

00 

1 















in 

1 

CM 

CM 

pH 

IA 

IA 

CM 

d 

IA 

O' 

IA 

o 

rv 

rv 


z 

» 

*A 1 

O 

pH 

pH 

O 

O' 

00 

NO 

d 

pH 

O' 

IA 

m 

rv 

O' 

X 

o 

1 

rv 


fv 

rv 

no 

NO 

NO 

NO 

<o 

in 

m 

d 

IA 

fj 

r 

1 















< 

I 
















1 

CM 

IA 

00 

O 

rv 

© 

O' 

in 

O' 

H 

in 

•d 

CM 


-J 

1 

~> 1 

pH 

PH 

o 

o 

00 

rv 

d 

CM 

O' 

rv 

o 

CM 

in 

rv 

►H 

1 

K 

IV 

rv 

rv 

«o 

nO 

-o 

<o 

in 

in 

in 

d 

IA 

CM 

Q 

1 















UJ 

1 

O' 

y£> 

IA 

o 

CM 

o 

d 

in 

in 

CM 

00 

d 

in 

d 

< 

1 

E 1 

CM 

90 

IA 

!A 

CM 

pH 

O' 

rv 

in 

IA 

rv 

pH 

d 

fv 

a: 

UJ 

l 

»0 

V? 

l* 

nO 


NO 

IA 

m 

m 

m 

vd 

d 

IA 

CM 

> 

l 















< 

1 
















1 

O' 

pH 

00 

pH 

o 

rv 

o 

O' 

'O 

o 

O' 

fv 




1 

< 1 

lA 

00 

00 

O' 

O' 

00 

00 

o 

m 

'd 

O' 

d 

00 

pH 


l 

1 

1 

IA 

lA 

IA 

IA 

lA 

m 

IA 

in 

m 

in 

d 

d 

IA 

IA 


1 

1 

CM 

90 

O' 

CM 

CM 

o* 

IA 

d 


m 

d 


*c 

CM 


1 

E 1 

IA 

O* 

o 

CM 

K) 

90 

d 

NO 

d 

IA 

pH 

•0 

90 

00 


1 

i 

1 

d 

d 

IA 

IA 

IA 

m 

IA 

m 

m 

in 

in 

d 

d 

IA 


1 

1 

O 

fv 

e 

O 

rv 

pH 

IA 

CM 

00 

© 

rv 





1 

U. l 

K> 

00 

PH 

IA 

d 

NO 

rv 

00 

00 

o* 

00 

rv 

d 

© 


1 

1 

1 

m 

lA 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 


t 

1 

CM 

o 


© 

IA 

CM 

o 

d 

nO 

m 

d 

CM 

rv 

CM 


* 

“> f 

d 

o 

CM 

in 

rv 

O' 

PH 

CM 

90 

Nd- 

in 

m 

IA 

pH 


1 

1 

CM 

lA 

IA 

IA 

lA 

IA 

d 

d 

d 

d 

d 

d 

d 

d 


>- 
















< M 
CK -J 

© 

o 

O 

o 

© 

© 

© 

o 

o 

© 

O 

© 

o 

o 


£* M 

o 

e 

IA 

o 

IA 

© 

IA 

o 

in 

o 

O 

© 




< M 


M 

pH 

CM 

CM 

IA 

IA 

d 

Nd- 

in 

NO 

rv 

00 

O' 


o 


i 



a 


-23A 


N •« N o. 

■n m m in 


iu 

o 

« v 

K< 

UQ 

> 


0 0 


n n M 
<« <« m' in m' «.* 


<z 
=>s 
z o 
ZE 
<< 



4ft 

s 

X 

2 

* 


4ft 

UJ 


4ft 

UJ 


o 

4ft 


4ft 
>- UJ 

3 h 

Z> 

X 


so 
ift rft 
O 

Z eft 


oe 

K 

ift 

UJ 

Of 

DC 

UJ 


Of UJ 

a. of 
4 ft O 
UJ 
KO 
o 

ocs* 




I 

< I 

I 

J 

I 

I 

I 

I 

r i 

t 

» 

I 

I 

I 

I 

U. I 
I 
I 


•ft 




<NJ 

•0 


•ft 

stf- 

•ft 


o 

10 

•ft 

Pft 

o 

stf- 

stf- 

4ft 

Ift 

Ift 

sO 

SO 

sO 

sC 

sO 


rj 


•o 

& 


•ft in & 
r% to 


ift ift it* it, 


ift 


CM 


CM nH 
MT <0 


*T 


Pft 

SO 


•0 

4ft 


•0 

sJ- 


CM 1% 
•ft PH 

sO *0 




•ft 

•0 


in 


"> i 


cm «o 
•ft to 

CM IM 


•ft ft. 
»H K) 

•ft •ft 


•0 ft * 

4ft ft. 

•ft Pft 


*■ «0 



*0 

PH fg 

s*' stf- 


Ift 

ft* 

<ft 

CM 

4ft 


PM 

4ft 

Ift 

sJ* 

s* 

•ft 


^ Ift 
*© CM 

4ft Ift 


ft* 

ft. 

gr- 



ft* 

Ift 

CM 

ft* 

•ft 

•ft 

CM 



* 

* 

•ft 


ift 

•o 

CM 

I 

o 




CM CM 


4ft 

Pft 


O 




AVERAGE DAILY TOTAL TERRESTRIAL INSOLATION ESTIMATES (KWH/SO. M) 


ui, tmmm 


%n 

>ui 



E 

sO 

c- 

lA 
UJ 
ZUJ 
< Of 

a o 

HlU 

a: Q 

UJ 

x <r 

tn >t 


UJ UJ 
A- O 



i 

r 



UJ 

O 


<>- 

cm 

in 

a% 

*0 

•0 

O' 



© 

•0 

A© 

ro 

© 

At 

Or < 

* 

» 

• 



• 

» 

* 

* 

♦ 

• 

• 

* 

* 

wo 

st 

*- 




*■ 

At 


At 

st 


At 

ro 

rO 


> 


_J 

A* 

e 

sO 

© 

SO 

o 

© 

ro 

o 

CM 

ro 

A* 

ro 

A* 

rs 

< 

z 

• 

• 

• 


« 

• 

• 

• 

• 

• 

• 

p 

* 

* 

33 

At 

© 

© 

pH 

O' 

m 

CM 

•0 

ro 

e 

ro 

P0 

90 

ro 

30 

ro 

in 

o 

O* 

*o 

•0 

O' 

•0 

A. 

m 

As 

in 

O 

ro 

ZC 

an 

© 

A. 

A* 

As 

As 

As 


As 

As 

A0 

m 

At 

CM 

< 

< 

pH 

pH 

PH 

pH 

#H 

pH 

H 

pH 

pH 

pH 

PH 

pH 

IH 

pH 


1 

1 

a* 

© 

O' 

A* 

ro 

no 

pH 

ro 

ro 

pH 

CM 

At 

o 

A O 


1 

ro 

•0 

O' 

H 

ro 

At 

aO 

As 

90 

O 

O 

o 

o 

90 


t 

1 

pH 

pH 

pH 

CM 

Cm 

CM 

CM 

CM 

CM 

CM 

ro 

ro 

CM 

CM 


1 

1 

© 

in 

«0 

O' 

•0 

© 


•t 

m 

St 

O 

in 

m 

sO 


1 

«0 

ro 

m 

A. 

O' 

pH 

ro 

Aj- 

m 

SO 

As 

A- 

A© 

At 


Z 1 


• 

« 

* 

• 

♦ 

• 

• 

p 

t 

p 

p 

* 

P 


1 

1 

pH 

rsj 

<v 

CM 

CM 

rO 

ro 

rO 

ro 

ro 

ro 

ro 

ro 

rO 


1 

1 

a* 

m 

pH 

At 

in 

ro 

•0 

O 

O' 

in 

A© 

m 

H 

m 


1 

pH 

p-* 

CT 

eg 

sr 

© 

A» 

Os 

O' 

o 

o 

o 

As 

rO 


O 1 

» 

• 

. 

* 

• 

• 

p 

* 

p 

p 

p 

t 

t 

* 


1 

t 

K> 

ro 

*■ 

gr 

O' 

sS- 

O’ 

© 

AS- 

m 

m 

At 


At 


1 

1 

rO 

so 

«0 

AO 

pH 

CM 

o 

st 

St 

o 

CM 

o 

ro 

in 


1 

An 

(Si 

At 

s0 

*0 

O' 

o 

O 

c? 

o 

•0 

At 

e 

st 


4A 1 

• 

* 

• 

p 


♦ 

p 

» 



» 



• 


I 

•t 

in 

m 

m 

in 

in 

sO 

aO 

aO 

A© 

m 

in 

in 

At 

X 

1 















>- 

1 















z 

1 

ro 

m 

o 

An 

pH 

© 

in 

SO 

pO 


fO 

m 

•t 

At 

o 

1 

ro 

sO 

to 

«0 

O' 

O' 

to 

As 

aO 

O* 

o 

At 

A* 

O' 

E 

< t 

• 

* 

p 

• 

* 

♦ 

t 

p 

p 

P 

* 

♦ 

p 

P 


1 

sO 

© 

© 

SO 

© 

sO 

© 

aO 

AO 

aO 

AO 

m 

At 

ro 

V 

1 















0Q 

» 















in 

1 

* 

© 

pH 

© 

pH 

® 

HA 

•o 

O' 

aS- 

CM 

As 

CM 

As 

H* 

1 

rO 

m 

in 

<t 

At 

ro 

pH 

O' 

s0 

AT 

•0 

O 

CM 

ro 

z 

-> 1 

• 

• 

• 

• 

» 

• 

• 

p 

p 

p 



p 

• 

3 

1 

A 

An 

r-N 

An 

An 

As 

As 

AO 

A© 

a£> 

m 

m 

At 

rO 

o 

1 















E 

1 















< 

1 
















t 

o 

sO 

ro 

A. 

A- 

ro 

© 

in 

pH 

aO 

«o 

o* 

in 

© 

>- 

f 

•0 

•0 

•0 

A* 

SO 

in 

ro 

pH 

O' 

A© 

o 

ro 

A© 

O 

pj 

“> 1 

• 

f 

* 

• 

t 

* 

• 

p 

• 

• 

• 

• 

♦ 

• 

AH 

I 

sC 

© 

A© 

sO 

© 

sO 

sO 

aO 

in 

in 

m 

At 

rO 

CM 


1 















o 

1 















UJ 

1 

ro 

At 

in 

CM 

sO 

As 

•o 

P0 

o 

o 

CM 

in 

© 

m 

o 

1 

© 

o 

© 

© 

O' 

•0 

As 

m 

st 

CM 

As 

pH 

in 

•0 

< 

E 1 

• 

• 

• 

• 

• 


* 

p 

• 

• 

• 

• 

* 


Of 

1 

in 

sO 

s© 

AO 

m 

•O 

in 

in 

in 

in 

At 

At 

ro 

CM 

UJ 

1 















> 

I 















■< 

1 
















1 

m 

An 

ro 

A. 

•0 

SO 

CM 

St 

A?- 

p^ 

O' 

As 

As 

0 


1 

•0 

© 

H 

pH 

pH 


pH 

O 

O' 

•0 

At 

O 

m 

o 


< 1 

♦ 

♦ 

» 

• 

p 

* 

• 

p 

p 

p 

P 

p 

p 

p 


1 

1 

•t 

m 

m 

m 

m 

in 

in 

m 

AS- 

aS- 

At 

At 

ro 

ro 


1 

1 

a* 

m 

o 

CM 

CM 

O' 

ro 

in 

*o 

o 

in 

O 


H 


1 

A* 

H 

ro 

At 

m 

in 

aO 

«o 

aO 

'O 

<t 

pH 

•0 

At 


E 1 

• 

♦ 

# 

* 

« 

• 

• 





* 


P 


t 

l 

rO 

At 

© 

O' 

O’ 

st 

© 


Atf* 

St 

At 

At 

ro 

ro 


1 

1 

© 

sO 

ro 

•0 

pH 

CM 

pH 

P0 

ro 

m 

At 

At 

in 

O 


1 

At 

•0 

© 


ro 

At 

m 

m 

a£> 

sO 

A© 

in 

»o 

© 


U. 1 

• 

♦ 

• 

t 

* 

• 

• 

p 

• 

t 

P 

• 

p 



1 

! 

CM 

CM 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 


1 

1 

rO 

ro 

pH 

A. 

CM 

sO 

« 

to 

As 

ro 

pH 

o 

CM 

© 


1 

© 

o 

(si 

ro 

in 

SO 

A. 

m 

O' 

o 

pH 

pH 

© 

to 


•*> 1 

* 

♦ 

* 

• 

• 

* 

# 

p 

p 

p 

p 

■ 

p 



1 

1 

pH 

CM 

(si 

CM 

CM 

CM 

CM 

CM 

CM 

ro 

ro 

ro 

rO 

CM 


>- 
















HC H 

e 

e 

© 

o 

C» 

© 

o 

O 

O 

O 

O 


O 

© 


Of -J 

• 

* 

* 

• 

* 

• 

• 

• 

• 

* 

• 

• 

• 

* 


oc m © o in .>iAoiAoin©ooo© 




C-236 


